Thermal fluctuations in pinned elastic systems: field theory of rare events and droplets 
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Using the functional renormalization group (FRG) we study the thermal fluctuations of elastic 
objects (displacement field u, internal dimension d) pinned by a random potential at low temperature 
T, as prototypes for glasses. A challenge is how the field theory can describe both typical (minimum 
energy T = 0) configurations, as well as thermal averages which, at any non-zero T as in the 
phenomenological droplet picture, are dominated by rare degeneracies between low lying minima. 
We show that this occurs through an essentially non-perturbative thermal boundary layer (TBL) 
in the (running) effective action r[u] at T > for which we find a consistent scaling ansatz to 
all orders. The TBL describes how temperature smoothes the singularities of the T = theory 
and contains the physics of rare thermal excitations (droplets). The formal structure of this TBL, 
which involves all cumulants of the coarse grained disorder, is first explored around d = 4 using 
a one loop Wilson RG. Next, a more systematic Exact RG (ERG) method is employed, and first 
tested on d = models where where it can be pushed quite far. There we obtain precise relations 
between TBL quantities and droplet probabilities (those are constrained by exact identities which 
are then checked against recent exact results). Our analysis is then extended to higher d, where we 
illustrate how the TBL scaling remains consistent to all orders in the ERG and how droplet picture 
results can be retrieved. Since correlations are determined deep in the TBL (by derivatives of r[u\ 
at u = 0), it remains to be understood (in any d) how they can be retrieved (as m = 0"*" limits in the 
non-analytic T — effective action), i.e. how to recover a T = critical theory. This formidable 
"matching problem" is solved in detail for d = 0, A'' = 1 by studying the (partial) TBL structure of 
higher cumulants when points are brought together. We thereby obtain the /3-function at T = 0, all 
ambiguities removed, displayed here up to four loops. A discussion of the d > 4 case and an exact 
solution at large d are also provided. 



I. INTRODUCTION 

Complex and disordered materials are often domi- 
nated, in their static properties, by rare but large fluctu- 
ations, and in their dynamics by ultra-slow relaxations. 
Such behavior occurs in a variety of systems involving 
many interacting degrees of freedom, from supercooled 
liquids to spin glasses^. Understanding the nature of 
the low lying excitations in these highly complex sys- 
tems and even more so, obtaining an analytic theory of 
the low temperature behaviour in realistic dimensions, is 
a formidable task. A class of problems which shares some 
of the basic physics and seems more amenable to analytic 
approaches are elastic manifolds pinned by quenched dis- 
order. They are of interest for numerous experimental 
systems, such as interfaces in magnetsSi^, contact lines of 
fluids wetting a rough substrate^, charge density waves^ 
(CDW), Coulomb glasses^ and Wigner crystals^iS in het- 
erojunctions and on the surface of heliumiS^, or vortex 
lattices in superconductor o^^i^^i^'^ . It seems likely that 
a broader class of problems, including random fleld sys- 
tems, can be understood by similar methods. 

In this paper we will thus focus on the low temperature 
equilibrium behaviour of elastic manifolds pinned by a 
random potential. Such systems are usually described 
by a displacement fleld u{r) which measures deformations 
away from a reference ordered state in the absence of any 
thermal, disorder induced, or quantum fluctuations. In 



their classical version they are typically described by a 
model Hamiltonian 



H[u]= J d'^T 



]^\Vu\^ + W{u{v),v) 



(1) 



where the random potential W{u, r) can be chosen, 
for concreteness, Gaussian with zero-mean and variance 
W{u,r)W{u' ,r') — R{u—u')5^'^\r—r'). We indicate here 
and in the following a disorder (ensemble) average by an 
overline. At equilibrium and low temperature the elastic 
manifold is usually pinned by disorder into non trivial 
rough conflgurations. In these so called glass phases, de- 
formations scale with size L asu ^ L'^ with roughness ex- 
ponent Q, and are dominated by the competition between 
elastic and disorder energies. Different physical cases are 
described by varying longitudinal (i.e. internal) dimen- 
sionality (c?), transverse (i.e. embedding) dimensionality 
{N — the number of vector components of m - we focus 
on = 1 for simplicity here), and R(u). For instance 
R{u) is a periodic function in the case of CDW or lat- 
tices, while in the case of interfaces in a ferromagnet it 
is short range for random bond disorder, and long range 
{R{u) ^ \u\) for random field disorder. 

As in other energy dominated glass phases in related 
systems, in the pinned elastic medium temperature is for- 
mally irrelevant in the RG sense as the ratio of typical 
temperature to energy scales flows to zero as Tl ^ L^^ . 
Here, the exponent 6* = d— 2-f 2C (> 0) as a consequence 
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of the so called statistical tilt symmetrjsi^ii^ (STS) , which 
expresses that the probability distribution of the disor- 
der is translationally invariant. However, because of the 
multiplicity of low energy configurations, thermal effects 
are rather subtle and do in fact dominate the behaviour 
of some observables at small but non zero temperature. 
The basic mechanism is that even if local quasidegen- 
eracies of ground states may occur only rarely, they can 
induce large deformations. 

A minimal physical description of such thermal ex- 
citations in disordered glasses is provided by the phe- 
nomenological droplet picture^S*^. In its simplest form, 
it supposes the existence, at each length scale L, of a 
small number of excitations of size Su ^ L'^ above a 
ground state, drawn from an energy distribution of width 
5E ^ with constant weight near 6E = 0. This means 
that with probability ^ T / this excitation can be ther- 
mally active with 5E ~ T. The consequence is that 
static thermal fluctuations at a given scale are dominated 
by such rare samples/regions with two nearly degenerate 
minima. For instance, the (2n)*'' moment of u fluctua- 
tions is expected to behave as 

- ~ c„(T/L«)L2<. (2) 

as a straightforward consequence (angular brackets de- 
note the thermal average with respect to the canoni- 
cal distribution). This simple picture seems to describe 
models such as Eq. ^ relatively well, at least in low 
dimensions r^iiSiSS. The statistics of low energy excita- 
tions has been much studied numerically recently in other 
systems such as low dimensional spin glasses^^. The 
droplet picture has been refined in many ways to account 
for more complex systems, and it is still controversial if 
this scenario fully describes systems such as spin glasses. 
Issues such as the fractal dimension d/ of these droplet 
excitations ("fat" droplets) and the possible global mul- 
tiplicity of true ground states are still debated. In any 
case, it is clear that the basic mechanism of stochastic lo- 
cal ground state quasi-degeneracies must be part of any 
realistic analytical theory of such systems. 

Similar physics occurs in the dynamics of these sys- 
tems. Although we will not address it in details here, it 
is useful to study, in parallel to the statics, equilibrium 
dynamics." The Langevin equation of motion is 

r^dtUrt = Vlurt + F{urt,r) + C(r, t), (3) 

with friction 77, thermal noise C(^: t) with 
< (:ir,t)C{r',t')) = 2r^TS'^{r - r')6{t - t'), and ran- 
dom force F{u,r) — —dV{u,r)/du (noting u{r,t) = Urt) 
with correlator denoted A(u) = ~R"{u). Equilibrium 
dynamics at non-zero temperature T > provides 
an equivalent way, through the fluctuation dissipation 
relations, to study the statics. Extension of the droplet 
picture to the dynamics supposes that equilibrium 
dynamics in the glass phase is dominated by thermal 
activation between the thermally active low lying 
quasidegenerate minima, with typical barriers scaling as 



Ub ^ ■ Little is known about the distribution of these 
barriers, but there is some evidence^^ that ip k, 0. 

Despite the wide applicability and extensive theoreti- 
cal studies of the elastic model, Eqs. (|ll3f) . there are few 
analytical results for the ground state properties, and 
even less for thermal excitations. There are exact so- 
lutions in the mean-field (A^ ^ 00) Hmife^iiSSiSSiSl and 
for fully connected modcls^'*'^^ (i.e. the d ^ 00 limit). 
These limits are interesting and probably capture part 
of the physics but they also have peculiarities which are 
still not fully understood. In both, thermal activation 
over divergent barriers {Ul ~ L^) is not included, and 
the infinite N limit (or Cayley trea^S*^) leads to global 
state degeneracies (with 5E ~ T) - replica symmetry 
breaking (RSB)"^^ - which may, or may not, be lifted at 
finite A^. For finite N oi d exact results exist only in two 
cases. There are some exact results for the roughness 
and the free energy distributions of the 1-1-1 dimensional 
directed polymer problem (d = 1 A^ = 1 case )"'^--" --^-'-^-- 
and the related Burgers turbulence, and KPZ growth or 
exclusion problem322i^2ii£. The d = limit of a particle 
in a random potential is simpler but also interesting as 
in some cases the full structure of the zero temperature 
fixed point can be solved exactly and many results ob- 
tained. This is the case for the A^ = 1 Sinai or so called 
toy model (i.e. RF) case {0 = 2/3): there the droplet pic- 
ture is "exact" and the droplet probabilities are exactly 
known^ii^. Other interesting results have been obtained 
for any A^ in the marginal case [9 = 0) of a log correlated 
R{u), where RSB seems to occur 

A potentially promising method is the field-theoretic 
Functional Renormalization Group (FRG). It was orig- 
inally developed at T = to lowest order in an ex- 
pansion in e = 4 — d ^ 1 (one loop) for the statics^ 
and later for depinning'*^ (around the upper critical di- 
mension due =4). It is called "functional" since one 
must follow the RG flow of the full correlator R{u) which 
becomes marginal in d = 4, as u becomes dimension- 
less. It is shown that R{u) becomes non-analytic be- 
yond the Larkin length as the random force correlator 
A(u) — —R"{u) develops a cusp singularity at u = 0. 
Physically this cusp is associated to the existence of many 
metastable states and its meaning has been discussed 
jp_27, 44,46^ At the depinning transition it yields a non-zero 
threshold force as fc ^ |A'(0''")|. Most importantly, the 
non-analyticity of the effective action seems to allow to 
evade the so called dimensional reduction phenomenon 
which renders perturbation theory in an analytic R{u) 
trivial and yields the incorrect result C, = e/2. Instead, 
in the space of non analytic functions, well behaved fixed 
points do exist, yielding reasonable values for the expo- 
nents (C = e/3 for random field disorder, C, = 0.283e for 
random bonds^ and a logarithmic growth of displace- 
ment for CDW and Bragg glass^^). 

An outstanding but crucial question is the inher- 
ent consistency of this FRG method given this non- 
analyticity, a highly anomalous occurence in conven- 
tional field theory. In particular, the occurence of non- 
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analyticity in the effective action at finite scale usually 
does not occur, because fluctuations average out any sin- 
gularities associated with multiple local minima. Here, 
the zero temperature field theory contains only disorder 
(sample to sample) and not thermal fluctuations, and 
these considerations do not apply. Instead it is tempting 
to compare the integration over fast modes to some (iter- 
ative) minimization procedure, which can then produce 
non-analytic energy landscape for fast modes ("shocks"). 
This picture can be made precise in simple single mode 
d — "toy RG" models^. The emergence of such shocks 
is well known in the equivalent Burgcrs-KPZ system in 
d — 1. In any case, the usual justification of pertur- 
bative renormalization group calculations (i.e. e expan- 
sions) must be re-examined. We will comment further on 
the physics associated with this structure below. 

A sign of the difficulties inherent in justifying the ap- 
proach occurs already in attempts to extend this FRG 
method beyond the one loop approximation, and also at 
non-zero temperature. Formidable difficulties appear al- 
ready at two loop order for the beta function at T = 
(and even in the one-loop corrections to the two point 
correlation function). Because the effective action is 
non-analytic, the renormalizcd perturbation theory is 
amhiguouiM'. An exact solution of the FRG at large N 
to any number of loops shows that special care has to 
be given to how any specific correlation is precisely de- 
fined, i.e. how the system is preparedSi. For N — 1 
quasistatic depinning the ambiguity can be fully lifted, 
as at any steady state velocity w > 0, the elastic mani- 
fold always moves forward. The resulting field theory is 
found rcnormalizable to two loop**^, and the predictions 
fully supported by numerical simulations. A solution of 
the ambiguity for the statics, based on a requested (and 
hoped for) renormalizability property, was proposed in*^. 
An alternative route, technically more complex but more 
transparent physically, which we will follow here, is to 
study the system at non-zero temperature T > where 
the action remains analytic and no such ambiguities ap- 
pear. 

Extensions of the FRG to non-zero temperature have 
been studied previously to one loopiSiSSiSiiSi. Since the 
renormalized temperature Tl ~ tends to zero, one 
expects that the asymptotic fixed point is still the T — {) 
one. This convergence was found to be highly non- 
uniform. From the FRG flow equation truncated to follow 
the second cumulant alone ignoring feedback from higher 
cumulants, it was found that A{u) = —R"{u) indeed re- 
mains analytic and differs from the T = flxed point 
function only in a thermal boundary layer u ^ where 
the cusp singularity is rounded. A highly desirable con- 
sequence of this early one loop analysis^Si^S*^ was that 
the correct scaling for creep barriers Ub ^ was recov- 
ered as a consequence of the boundary layer, from the 
divergence of — i?""(0) ~ I/^l, the quantity which en- 
ters the friction (ie the time scale) renormalization. This 
very same divergence however, is also very troublesome 
in considering any extension to higher number of loops 



(as was found inS) since one easily realizes that adding 
one loop extracts one power of I/Tl, thus adding one 
new divergence. Until now the status of the boundary 
layer beyond one loop thus seemed rather unhealthy. 

The occurrence of these technical issues is not surpris- 
ing when one reflects upon the non-trivial physics that 
the fleld theory must describe. On the one hand, typi- 
cal samples and their ground states - which control static 
correlations at T = - are critical in a fairly conventional 
sense. For instance, they are described by scaling expo- 
nents (C, 0, ■ ■ .) and amplitudes, which apparently (hope- 
fully!) can be computed systematically in an e expansion. 
On the other hand, thermal fluctuations of the manifold 
are dominated by rare low-lying metastable conflgura- 
tions, which play no role in the T = theory. The distri- 
bution (over disorder) of the thermally fluctuating part 
of the correlations (e.g. {{u — (u))^)) is thus extremely 
broad, and must be encoded somehow also in the field 
theory at T > 0. Both types of quantities are formally 
defined extremely similarly as derivatives of the effective 
action (e.g of the renormalized R{u)) at u = 0. Never- 
theless, they scale very differently and represent rather 
different physics. The coexistence of these two features 
strongly argues for an unconventional structure of the 
field theory. 

In this paper we solve these issues in d = and pro- 
pose a consistent form for the effective action of the field 
theory in any dimension to any number of loops. We es- 
tablish the clear connection between the formal thermal 
boundary layer and the physical droplet picture in any d 
(see^ for an early calculation). Thus we obtain a "field 
theory of droplets", i.e. we show that the present ex- 
tension of the field theoretic FRG captures correctly the 
basic feature of glasses, rare events due to quasi-ground- 
state degeneracies. The mechanism is quite non-trivial. 
A short account of these results has been given in^'^ , to- 
gether with a discussion of the thermal boundary layer 
structure which also arises in the dynamics and yields 
there to fluctuations of barriers, as discussed in detail 
iiM. 

First, in Section^ within a Wilson formulation, we 
derive and carefully analyze the full one loop FRG equa- 
tions. We discover that the thermal boundary layer ex- 
tends to all cumulants of the disorder. Thus the trunca- 
tion to the second cumulant performed in previous stud- 
ies is not strictly consistent at T > 0. It does however 
give the correct overall scaling, since we find that a con- 
sistent ansatz can be found for all cumulants. This results 
in a hierarchy of FRG equations for the scaling forms of 
each cumulant or the random potential, with good prop- 
erties. Unfortunately, there is no longer a small parame- 
ter (e has disappeared as a small parameter!) which pre- 
vents a perturbative solution. We must thus assume that 
there exists a global solution, which we call the Thermal 
Boundary Layer Ansatz (TBLA). 

To further explore this novel theory we turn to the 
Exact RG method (ERG), explained in Section ITTll The 
ERG follows the full (non local) effective action func- 
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tional under coarse graining. To better understand the 
physics, we start in Section ll VI with a detailed analysis 
of the d = case, where the effective action becomes 
a (replica) function rather than a functional, and can 
be studied without any truncation. There, the droplet 
picture is exact, and a variety of exact results are avail- 
able for comparison with the methods developed here. 
From the TBLA and the ERG, we derive exact relations 
for droplet probabilities. These are tested on the exact 
solution of the toy model, and are found to pass this 
non-trivial check. In Section we introduce and dis- 
cuss the problem of matching the zero temperature and 
boundary layer regimes. In particular, there are a variety 
of correlation functions describing the critical (typical) 
correlations of the ground state, and as such should be 
determined physically by the T — non-analytic field 
theory. Formally, however, they are obtained as indi- 
cated above from derivatives of the effective action at 
M = 0, deep within the TBL. Thus, some information 
should be transmitted through the TBL. We show how 
these correlations are indeed determined from particu- 
lar limits of the (non-analytic) zero temperature "outer" 
solution. This is accomplished through a careful consid- 
eration of the Partial Boundary Layer (PBL) structure 
that obtains when any two relative replica displacements 
{ua — Ub) become of order T^. Thus we will describe 
the full solution of the ambiguity problem (for d — and 
= 1) in the zero temperature beta function. As an 
example of the power of this prescription, we calculate 
the beta function and various correlations as a series in 
R up to four loop order. 

With this physical understanding of the FRG equa- 
tions in the thermal boundary layer regime in hand, we 
then turn to d > in Section IVII We show that a 
fully consistent boundary layer ansatz can be found for 
the whole non-local effective action, and is self-consistent 
when applied to the ERG equations. This result provides 
a non-trivial "solution" to the above-mentioned disease 
of proliferation of l/T^ divergences by resumming them 
order by order into a finite result. We show explicitly 
that the results from the droplet picture are recovered. 
Perturbative control (i.e the e expansion) is naively ex- 
pected in the outer non-analytic solution. The necessary 
extension of the matching understood in d = to obtain 
physical quantities (even zero temperature ones) remains 
a challenge. We discuss this issue and some perspectives 
on future work at the end of this section. Whether the 
method could account for other scenarios, more complex 
than the initial droplet picture, is also addressed. 

A panoply of appendices give various technical details 
and a description of related models. In particular the case 
d > 4 is examined in details: it is shown, through anal- 
ysis of the one loop FRG, exact solution at large d and 
qualitative arguments, that a phase with non- analytic 
effective action exists at T = and strong enough disor- 
der, but that the TBL remains of finite width, at variance 
with d < 4. 



II. BOUNDARY LAYER WITHIN THE WILSON 

FRG 

In this Section we establish the Wilson one loop FRG 
equations, which are a convenient starting point for our 
investigation of the thermal boundary layer. We include 
from the start all cumulants of the random potential since 
we do find, a posteriori, that they must all be considered. 



A. model, properties and notations 

From now on we consider the equilibrium statics at 
T > 0, defined through the equilibrium Hamiltonian in 
(Q. Formally, all static quantities of interest can be ob- 
tained from the replicated generating (partition) function 



Z[ha] = fl /i^^..exp[-^+ / 

a=l"' ^ 



ha{r)ua{r) 



(4) 



The number of replicas n can be taken to zero, or kept 
arbitrary as a bookeeping device. Since some of our cal- 
culations are independent of n, we will explicitly indicate 
when the limit ti = is considered. 

The connected cumulants of the random potential 
W{u,r) in Eq.jTlare denoted [k > 2): 



W{ui,ri)..W{uk,rk) 

= (-)'=5W(wi,---,7.fe)5'^(ri,...,rfc) 



(5) 



with S**^^^ (m, u') = R{u — u'). The cumulants higher than 
second are generated by coarse graining, and are thus 
included here from the start. We ignore for the moment 
correlations between different points in internal space: 
this choice and its extensions are discussed below. The 
probability distribution of the random potential is chosen 
translationally invariant, so the cumulants satisfy 



(mi A, • • • , Wfe + A) - ^t'^) (ui , • • • , ufe). 



(6) 



They are also even so that S''^'^'^ (— mi, • • • , — Ufe) = 
5^*^^ (mi, • • • , Ufc), the 5*'-'^^ being by construction fully 
symmetric functions of their arguments. 

Performing the disorder average one obtains: 



Z\h 



u-\^^V[u]), 



(7) 

(8) 



where we denote vectors in replica space as w = 
(mi, • • • , Un), and we define the characteristic function of 
the pinning disorder 



1 



k=2 ai---Ofc 



(9) 



which is a fully symmetric, translationally-invariant func- 
tion of the replica vector u. The property © results 
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in STS, e.g. in Fourier space \nZ[h = h] = In Z[h = 
0] — Tnh{q)h{—q) / {2q^) resulting in infinite sets of rela- 
tions between correlation functions and other important 
consequencesi^ii^ , to which we will return later. One 
consequence is that the "one replica" part of the action 
in JHJ is not corrected by disorder under coarse graining 
(which is also why one can set the elastic coefficient to 
unity) . 

Let us give some useful notations and properties. For 
notational simplicity we will often write: 

^('■H«i2...fc) = ^<'H«i,---,wfc) (10) 

for any function, and for a function of one argument 
R{ui2) = R{ui — M2). Since any function antisymmetric 
in any two replicas gives zero contribution to the action 
when summed over replicas, such functions may be added 
freely to any cumulant without physical consequences. 
This freedom will be used to simplify some expressions. 
Usually the cumulants S'^'^'-' will be taken fully symmet- 
ric in their arguments. It is then useful to define the 
symmetrization with respect to a given set of arguments: 

symi...fc0(wi...fc) = X! '^('"^(i)---'^(fc)) (11) 

Pfc being the set of permutations of k objects. Next we 
note two useful properties of the above decomposition 
of the action in cumulants, i.e. in number of replica 
sums Translational invariance © implies relations 
between derivatives, e.g.: 

k 

(^5,) =0. (12) 

This implies properties such as: 
symS'|o^o["i---fc] = 0, 

symS'^o^ o[ui...fc] = -(fc - l)sym5ji^ o[Mi...fc], (13) 
where here sym = sym^...;., and we denote everywhere 

(k) 

partial derivatives S'n/„2..„j. for derivatives in the i-th 
argument. 

Finally, there is a "gauge invariance" property: the 
change 

5W[ui...fe] -> 5W[ui...fc] +symi...,0(wi...fc_i)(14) 

yields an extra term in ((SJ proportional to an explicit 
n. Thus, if one considers n — 0, one can add freely to 
each k replica term any function of p < fc replicas. Phys- 
ically, this gauge freedom corresponds to the fact that 
the random potential may be shifted by a w-independent 
random function of x, changing only the free energy in a 
trivial way but leaving all displacement correlation func- 
tions unchanged. Explicitly, this shift is of the form 

W{u,x) ^W{u,x) +w{x), (15) 



where w{x) can have arbitrary cross-correlations with W, 
i.e. 

W{ui,xi) ■ ■ ■ W{uk, Xk)w{xk+i) ■ ■ ■ w{xp) 
^5,,...,^<i>(P'''\u,,---,Uk), (16) 

with arbitrary <i)(P^'^) . 

B. conventional scaling and one loop FRG 
equations 

Let us first summarize the conventional view on the 
behaviour of the above action ^ under coarse graining 
and RG as it has emerged from previous worka^^i^2i5iiS^. 
In the usual zero-temperature power counting, namely 
T b^T, X — > bx, and u — > 6^u, all the cumulants 
S^*^) for k > 2 are irrelevant, and the second cumulant 
R{u — u') is marginally relevant just below d = 4. Simi- 
larly, terms describing correlations between different in- 
ternal space points are strongly irrelevant, which justifies 
to restrict to ((SJ) with fc = 2 (in the field theory language 
the only needed countcrterms are for the local part of the 
second cumulant of the disorder). 

More precisely, within the Wilson approach with a run- 
ning UV momentum cut-off A; = Ae^', one makes the 
above action dimcnsionlcss S[u\ = AdS[uK^i] as: 

S[u] = / d'^fiY, -^IV.kI^ - V[v\l (17) 

a 

with f = rAi . We have introduced the ratio of thermal 
energy T to pinning energy which defines the running 
temperature: 

fi = AdTAl, (18) 
= d-2 + 2C>O. (19) 

The rescaled temperature T/ flows to zero, indicating the 
dominance of energy over entropy as appropriate for a 
pinned glass phase. For later convenience we also in- 
troduced Ad = Sd/i2nf = l/{2'^~^n'^^^r{d/2)). The 
rescaled generating function reads: 

V[u] ^ AdAfV[uA'l], (20) 

^["l^Efl^ E ^''■^["a,r--,"aj, (21) 
fc=2 I ' ai---afc 

and corresponds to the rescaled dimensionless local cu- 
mulants: 

Riu) = A-^Af-'^-^'^R{uA'i), (22) 

S^'^^ [Ua, , • • • , Ua,] - Al^-'^Af-'^'S^'^^ Af , • • • , Ua.Aj] • 

These equations embody the conventional scaling. 
Within the Wilson formulation one finds that for e = 
4 — c? <C 1 there are well defined fixed point functions 
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R*{u) ^ e <^ 1 associated to C also of 0(e). One also 
expects that for k > 2 the fixed point value of the func- 
tions S'('=)*(u) < R*{u) (more specifically S^'''>* e*^ < e 
for k > 2). This smallness has been argued to justify 
the truncation S'^'^^ ~ for k > 2 at one-loop, to lead- 
ing order in e at T = 0. Prior investigations of thermal 
effects have presumed this truncation remains valid at 
T > O^SiSiiS also to obtain the leading order in e. We 
now check this assumption. 

To capture this leading-order behavior, we therefore 
employ a one loop FRG analysis. A general one loop 
FRG scheme contains in fact much more than the mini- 
mal physics needed - in principle - for a consistent 0{e) 
analysis by the usual power-counting. However, we will 
pursue it in the general form in order to have a single 
calculation which encompasses both this simple regime 
and the thermal boundary layer. We emphasizes that at 
this stage it has mainly a heuristic value, since, as we find 
later, the restriction to one loop is not necessarily justi- 
fied within the boundary layer. Similarly, the considera- 
tion of non-local terms in internal space for the disorder 
cumulants will also, in the end, become necessary. We 
ignore for now all of these further complications, which 
will be captured by the more rigorous (but heavier) ex- 
act RG scheme studied below. The present much simpler 



one loop Wilson RG will give us a first good idea of what 
is happening in this theory. 

Expanding the action (jHJ to quadratic order in the 
"fast" modes with Fourier component non-zero only in 
the shell A/ < g < Aj+^i and performing the resulting 
Gaussian functional integral one obtains the differential 
flow equation for the function V{u) as: 



a,y(") = -iAXTrln 



rah 



T 
A 



p>i 



T 



Tr{{V")P) 



(23) 
(24) 



where the Tr denotes a trace in replica space only. For 
convenience we have introduced the parameters Cp = 
l/(2p) since an almost identical equation will be derived 
below for the exact RG in d = (with Cp = 1). 

We can now use the expansion in ^ and collect terms 
with a given number of replica sums on either side of 
(|23|l . The general form of this equation is derived in 
Appendix^ We give here only the result for the projec- 
tion of the FRG equation H23() onto two and three replica 
components: 



diR{u) ^{e~AC + Cudu)R{u) + fiR"{u) + 5^i^(0, 0, u) + -R"{uf - R"{0)R"{u) (25) 
diS<^^\u,23) = (2e - 2 - 6C + CuA^)S^^Hu,23) + sym S^^„{u,23) + ^R"(«i3)R"(«23)^ 

+ 2Sym^noo("ii23) (26) 



where here sym = sym]^23j and 7=1. Since i?"(0) does 
not feed into higher cumulants (for n — 0) we define: 

R{u) = R{u) - R"{0)u^/2 (27) 

Note that since additive (u- independent) constants in 
(^5)1 are pure gauge (for n = in the sense discussed 
above) we will in general ignore them below. Equations, 
with similar structure exist for higher cumulants (see Ap- 
pendix ^ for details). 



C. analysis within conventional scaling 

One first notes that H25I I26|l does not close, since there 
is a feedback of cumulant fc + 1 in cumulant fc, though 
only in a single term. At zero temperature, in the con- 
ventional scaling discussed above this should not be a 



problem. Let us first discuss that case and formally set 
fi =0 in (|25I26|I . We will call the solutions of the result- 
ing equation the "outer" solution. This terminology is 
motivated by the thermal boundary layer analysis which 
is performed below. 

1. zero temperature 

For T = the solution Ri{u) of (^5)1 (neglecting the 
third cumulant term S*^^^) converges to a zero temper- 
ature fixed point function^ R*{u), with a small set of 
universality classesiSiiS. For, e.g. W{u, x) periodic on 
the interval < u < 1, one finds C = and the force- 
force correlator A*{u) = -R*"{u) = Min„g2|[(u - n- 
1/2)^ — 1/12]. Its non-analytic behavior at small u, 

A*iu) = A*{0)-x\u\ + O{u^), (28) 
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is, however, super-universal, i.e. the same for all disor- 
dered elastic models. We denote x = I ^'(0^)1 — ^Xi 
where x is a 0(1) constant whose numerical value de- 
pends on the model. The physical significance of this 
cusp has been discussed by several authorai2Li4ii&. It is 
related to the existence of multiple metastable minima 
in the effective potential, since these are implied by the 
divergence of the mean squared curvature of the pinning 
potential, -A^^o(O) = {dlV{u)y = -l-oo. We will see 
further consequences below. Note that the second deriva- 
tive of (|25|l at u = yields the fixed point constraint 

-i?*"(0) = A*(0) = xV(e-2C). (29) 

Note that this will be the precise definition of the pa- 
rameter X in the remainder of this paper. The relation 
between x and the first derivative given above then holds 
only within the one-loop Wilson approximation, and is 
corrected at higher order. 

The neglect of iS^^^ then appears consistent. Given 
that the second cumulant flows to a fixed point R* ~ e 
then H26|) shows that the third cumulant also flows to the 
fixed point: 

-lR*"iu,2)R*"iu23)R*"M). (30) 

Thus, at least formally, s!^^q ^ e^. The other terms, in- 
cluding the fourth cumulant iS^lg ^ (and higher which 
have fixed points^^ similar to Ij^UI) ) can be neglected self 
consistently in that picture, as well S'^"^^ in the equation 
for R(u), as they yield only higher-order corrections in e. 

Although this picture is consistent at T = to the 
lowest order in e, closer examination shows possible am- 
biguities in the feeding term ^^^q(0,0,m) in Eq. if 
one uses (|30|) which make the next order in e more prob- 
lematic. In particular, differentiating H3U|) and taking 
V = ui — U2 small yields: 

S[%{v, 0, u) « ^R*"{u){R*"'{ur - R*"'{vr) 
-^R*"{v)R*'"{uf - jR*"{v)R*""{v)R*"{u). (31) 

It turns out that in this expression has well defined and 
coinciding limits as v — *■ 0^. However, in the formal 
FRG equation ^ s[%{v,0, u) is evaluated exactly at 
?; = and generally the replacement of that expression 
by either limit w — *■ 0^ or u — > 0~ may not be valid 
in particular renormalisation scheme. More general am- 
biguities arise at higher orders and cannot be so easily 
resolved regardless of the scheme. Even in the present 
Wilsonian formulation, treatment of this term alone is 
not consistent since other two-loop contributions arise at 
order e^, thus it is not clear that the smoothness of this 
limit is sufficient to guarantee a well behaved e-expansion 
even to this order. Studies beyond Wilson^ are discussed 
elsewhere, so we will not address them here. Instead, 



we follow a different route and focus on T > 0, where no 
ambiguity appears, and first summarize the naive T > 
analysis still to lowest order in e. A systematic study 
of the resolution of these ambiguities by connecting zero 
and non-zero temperature quantities will be performed 
for the case of d = in Section HVI 

2. non-zero temperature 

Reinstating the T/ terms in (|25I26() . the naive T > 
analysis still assumes that S''-^-' <C ei? and converts l|25|) 
into a decoupled fiow equation for Ri(u). Since T; = 
e~^' flows rapidly to zero, Ri{u) does converge to the 
T = fixed point. However, since at T > the function 
remains analytic around u = the convergence to the non 
analytic T = fixed point function is highly non-uniform. 
One finds^Si^ that for u ^ T;/e the cusp singularity is 
rounded in a thermal boundary-layer (TBL). Its detailed 
solution yields two regimes at large I (precisely Ti <C e^): 

7?/(u) -> i?*(u) for u-0(l)>f//e, (32) 

Ri{u) = ^R*"{0)u'+f^{ex)-Mu), 

for u=^^0(l), (33) 
i; 

where we introduced the TBL scaling variable u and scal- 
ing function r{u). The matching of the thermal boundary 
layer H33() to the zero temperature behaviour outside l|28() 
requires that: 

r"(u)-|u|, for|u|>l. (34) 

Substituting the TBL form (fS^ in if^ one finds that 
(up to additive constants) the leading terms are 0{Tf) 
yielding the constraint that 

(e - 2C)^^y + \{r"{u) - r"{Q)f + r" {u) 

(35) 

is a constant, which must thus equals r"(0). Note that 
the diR term is 0{Tf /e^) and the rescaling term (apart 
from the one acting on the purely quadratic part of R) is 
even smaller as 0{Tj^/e) so both can be safely neglected 
at large I such that f/ < e^iS. Using ^ yields^ 

r"(u) -r"(0) = VTT^- 1, (36) 

which is smooth for \u\ of 0(1) and has the desired match- 
ing property (|34|l . Note that the flow of R'/iO) is not 
fully specified by the solution. However the existence of 
the above single TBL together with the matching implies 
that the value to which it converges must be equal to the 
T = fixed point value, as can be seen from (|35|l (a 
different asymptotic value would require a more complex 
TBL scenario, not supported here). 
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Thus the problem of matching is successfully solved 
within this approach. The convergence to was stud- 
ied in^, where it was found that an expansion in higher 
powers of T/ could be solved order by order. In^^ it was 
found that the form H36I) was particular to Wilson, as 
other one loop schemes (with the same truncation to 
second cumulant) give other analytical forms. Although 
r{u) was found non universal, its large u H34|l and small 
u behaviour must be. For instance, the divergence of the 
mean curvature is constrained to be 

fiRnO) ^ (37) 

(i.e. r""(0) = 1). This simple consequence of the TBL 
was found to yield for susceptibility fluctuations the same 
scaling as predicted by droplets^^. Also, since a similar 
analysis (and truncation) in the dynamics^^ yields that 
the friction grows as dilnrji = i?""(0) the TBL nicely 
recovers the scaling of barriers Ub ^ expected in the 
creep motion. Thus there must be some truth in this 
TBL scaling, although this simple analysis is far from 
being sufficient, as we now show. 

Before doing so, we note that above four dimensions, 
the situation is rather different, and discussed in Ap- 
pendix^ where this case is analyzed by FRG and large 
d methods. For a weak random potential with smooth 
correlations, no cusp arises even at T = 0, and there is 
therefore no TBL. For stronger disorder or a sufficiently 
non-smooth potential, the zero temperature cusp does 
arise. In this case, however, the TBL has a width ^ T 
instead of T; , and so does not scale to zero with the infra- 
red cutoff. The physics of this difference ~ which we be- 
lieve is associated to the short-scale nature of metastable 
states for d > 4 - is discussed in Appendix IhI 

D. thermal boundary layer: the systematic analysis 

We need now to estimate the magnitude of the term 
s['^Q in the second cumulant equation (|25l) . in the TBL 
region u Ti/e, which has not been done before. For 
that we must carefully reexamine the equation for the 



third cumulant 1)26(1 . First one sees that S'-^\ui23) will 
also have a TBL form for m ~ Ti/e since the 0[{R")^] 
terms in (|26|l feed their boundary-layer scaling forms into 
S^^\ Naively balancing these terms with the rescaling 
part (« -2S^^^) suggests S^^^ ~ f^. This is, however, in- 
consistent, since the remaining terms linear in S*'^' would 
then be 0{Tf). Instead, the only consistent TBL ansatz 
for the third cumulant is: 

S^^\ui,U2,U5) = {ex)-^fts^^Hui,U2,us) (38) 

The second derivative term feeding into Eq. H25() then 
reads: 

sf^^{0,Q,u) = ffs%{Q,Q,u) for u^O{l) (39) 

and is thus of the same order 0{T^) than the other sur- 
viving terms in the boundary layer for R{u). Thus the 
above described result for r(u) appears as an (uncon- 
trolled) approximation, since we now have two coupled 
equations for the TBL scaling functions r and s^'^h The 
equation H25|) now leads to: 

0=-y + i(r"(u)-r"(0))2 
+r"(i7)-r"(0) + ,sfil,(0,0,u), (40) 

while s^^-* (wi, {t2, M3) satisfies a TBL equation simply ob- 
tained from (|26|l by setting S^^'^ s'^^\ R ^ r, Ui ^ Ui, 
T; — > 1, suppressing the linear rescaling term and setting 
the l.h.s. to zero. 

Of course we must now worry about the fourth cumu- 
lant term 5'{ioq in the equation for S^^\ so in fact there is 
an infinite set of coupled equations and the scaling forms 
in I(33I38|I must be extended to all cuniulants. We will 
now establish the consistent form of all cumulants and 
derive the TBL equation for all of them. There are a few 
subtleties in that construction. 

Let us first write the one loop FRG equation for the 
fc-th cumulant in a very schematic form: 



k 



J2 s^"' 



■ s 



fe-1 



^'11 ^('=i)"...5('='")"(41) 



where O denotes the rescaling linear operator, the two 
next terms are tadpoles, and each of the last two terms 
is a loop with m > 2 disorder vertices. The counting 
of T factors results from there being m propagators per 
loop, hence a factor T™, and a factor T"*^' associated 
to each /ci-replica vertex S'^'^^K One easily sees that only 



these terms can appear. This property and their detailed 
structure is worked out in Appendix but is irrelevant 
for the present discussion. Among them there is a feeding 
term of 0[(R")'=]. 

A look at (|41|l suggests that the natural generalization 
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of H33I38|I . namely 

• ■,Uk)^ (xe)-2rj'-+^s«(t2i, ■■■Uk) (42) 

should produce a consistent TBL equation with all terms 
in H41() of order 0{Tj^), except the rescaling and left hand 
side which are OiT^'^^ / e^) and thus disappear from the 
TBL equation as T/ <C . This is essentially correct, 
apart from one simple, but important, modification. In 
the TBL equation for r (|35l I4U|) one term survives in 
the rescaling part, coming from the quadratic part of 
R(u), and this is because this term scales naturally as 
eR*'"{0)u^ ^ (while all other even monomials scale as 
larger, hence subdominant, powers of T/). One can then 
ask whether one should consider similar terms in the TBL 
equation for s'^^\ which is 0{Ti ). The construction of 
the /c-replica terms is detailed in Appendix IbI One easily 
sees that for k = 2p even, the monomial of lowest degree 
one can construct which is a genuine fc-replica term is 
indeed of order u'^ and thus scales as Tl^. It reads: 

q2pUiU2-U2p, (43) 

where q2p = .§^^''^(0, • • • , 0), and is fully symmetric, even 
as it must be, and translationally invariant (see Ap- 
pendix . A /c-replica monomial must contain at least 



one of each replica, i.e. the product (|43|l . otherwise it is 
a pure gauge. No such term exist for k — 2p+l odd, and 
one shows that the term of smallest order in that case is 
of order u'''^^. For instance for the third cumulant, one 
checks that the term of order u"^, namely 

symi23W?W2M3 (44) 

is forbidden by translational invariance (see Ap- 
pendix . Thus the third cumulant starts at higher 
order S'^'^^ ^ u^. The conclusion is thus that the TBL 
form (|42|l is correct for odd cumulants but for even cu- 
mulants it should be replaced by the correct TBL form: 

S^'^p\ui- ■ -Uk) = q2pUiU2-.U2p (45) 
+ m-'Tf'^'P^\u,---U2p). 

Having determined the scaling for all cumulants within 
the TBL, we briefly return to the third cumulant TBL 
equation. Using (|46|l . we can express the feeding terms 
from the fourth cumulant in their appropriate TBL form. 
Remarkably, as for the second cumulant, a formal an- 
alytic solution for s^'^'> is possible (and this appears to 
generalize to all cumulants). Using the gauge freedom 
and STS relations (|13|l one directly rearranges the third 
cumulant equation into the form 



Sno('"i23) - Sno("ii3) = — ],,- s 1^(^123) + a(Mi3) + a(u23) + symi23 \]-r" {uis,)r" {u2s.) 
+lr"{u^2)r"{u^:,f - |r"(fii2)r"(S23)r"(ui3) + ^4t'oo("ii23)] }, 



(46) 



where G'(t2i23 ) ^ G(i)(mi23) + G(2)(y^23)^ ^here 
are translationally invariant functions symmet- 
ric in 1 <-> 2 and respectively antisymmetric in 1 3, 
2 <-> 3, and a{u) is an even function. These are con- 
strained by the requirement that s'^', or equivalently 

(3) . 

s\-{^ is a symmetric function, and in the latter case, fully 
gauge invariant. Although this has been written as a 
formal solution for s'^^\ it clearly involves the fourth cu- 
mulant quantity Sxioo(''^ii23)- Thus "solutions" of this 
form do not close, and form a coupled hierarchy. Indeed, 
it is equally valid to view H4t)|l as expressing the fourth 
cumulant feeding term in terms of lower cumulants. This 
and other similar relations are used later in the paper to 
simplify some expressions. 

Let us write the TBL scaling form in a more compact 
way in terms of the reseated characteristic function de- 
fined in (|21|l . It takes the form: 



Vi{u) 



p>l ai 



E 



T, 



v{u), 



(47) 

where we have defined the boundary-layer replica vector, 
and scaling function: 



u = exu/fi 

OO 

-(-^)-E^ E 



k—2 ai---ak 



(48) 
(49) 



Note that all dependence in T/, which has been made 
apparent in (|47|l . factors in front of all cumulants and 
that the TBL characteristic function v(u) approaches a 
fixed point form at large I. Below we will discuss the 
physical meaning of the set of (/-dependent) constants. 



/2p = ^f/.U0,---,0)/(xe)^^ 



(50) 



with 

1 



/2 = — i?"(0)/(xe)^ (remember that UaUf, = 
^{ua — UbY up to gauge terms). 
Inserting the I dependent form (|47|l into H23|) . taking 
into account all I dependent rescaling factors, and dis- 
carding terms small as Ti <^ e^, one gets the TBL-FRG 
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equations which must be satisfied by the constants f2p 
and the function v at their fixed point. The Tr In term 
in 1)23(1 takes the form: 



Tr In 

Mab ^ 



(5'^'' - aMab - dadbv{u) 

'-dadb[^f2p ^ Ua 
p> 1 ai , . . ,a2p 



(51) 
(52) 



with a — (ex)^/Ti, and da — d/dua here. 

However we see that because the monomials associated 
to each f2p is such that one rephca appears only once, the 
matrix Mab is independent of replica a and b. Thus in 
the expansion of ((51(1 to any order in M , any power of 
M or any product of the form: 

J2 Mabdbdcv{u) = Mab ^ dbdcv{u) = (53) 



vanishes when traced, because of translational invariance 
((T^ . Thus M can be set to zero (for n = 0) in |(SIJ and 



we finally obtain the non-trivial one loop Wilson fixed 
point TBL equation in a very compact form: 



p>l ai,..,a2p 



■^ia2p-^Trln 



rafc 



dadbv{u) 



= 0, 

(54) 

where all I dependence has disappeared, the trace be- 
ing over replica indices. The first term is, as discussed 
above, the only one surviving from the rescaling part and 
contains the scaling eigenvalues: 



X2 =e-2C = 0(e), 
X2p^d-2p{e-C)^0{l) 



p> 1. 



(55) 
(56) 



((54(1 is a highly compact form which, expanded in the 
number of replica sums, yields coupled boundary layer 
equations for r and s^-^\ and all higher cumulants. It 
is useful for later purpose to also give the schematic 
structure of these TBL equation for the coupled function 



oik) 



(ui,---,Ufc): 



for any k > 2, where, as in 1(411) all coefficients and de- 
tailed structure has been suppressed (they are identical 
to the one in 1(41(1 and given in Appendix IXjl . Note that 
the coefficients fk appear only in the first term and have 
no feedback in the others. Similarly, the naively first 
non-gauge term in s*^^p' proportional to mi • • • U2p (sub- 
dominant to f2p) drops out of the equation completely. 
Thus this term can be taken as zero, and s^^p^ thus is 
chosen to start as {t^P+^. 



E. physical consequences, correlations, droplets 

Let us now pause and analyze what has been achieved. 
We have derived, within the one loop Wilson FRG, cou- 
pled flow equations for all cumulants of the disorder. 
We have found a scaling form for its solution at large 
I (i.e. small T/) which exhibits a thermal boundary layer, 
consistent for all cumulants. It is expressed in terms of 
an infinite set of unknown parameters f2p and functions 
g(p) ^ These encode for the probability distribution of the 
coarse grained disorder, and obey the reduced fixed point 
equations (|5ll . [K7|l . 

The first observation is that there is no remaining 
small parameter in Eqs. ((54I57() . It is clear in the Tr In 
term. For the other terms it also holds if conven- 
tional (zero temperature) scaling holds. Indeed although 



/2 ^ 1/e, the product 2:2/2 = 1, and in the conventional 

zero T scaling S[l'^\iO, • • • , 0) - e^P, thus X2pf2p 0(1). 
The problem has thus become fully non-perturbative and 
one is thus in the uncomfortable situation where a full 
solution of this fixed point equation |(^ , El is required 
simultaneously for all cumulants of the disorder distri- 
bution, an infinite number of functions. Although some 
formal algebraic relations between cumulants can be ob- 
tained in the TBL, as discussed above for s^^\ no general 
solution of the full hierarchy appears possible. Thus at 
this stage we have no choice but to assume the existence 
of a global, well behaved solution to the TBL-FRG equa- 
tion and explore some of the consequences of this thermal 
boundary layer ansatz (TBLA). As we show below, it im- 
plies that a smooth zero temperature limit exists for dis- 
placement correlations determined by the f2k coefficients, 
i.e. a well defined limit exists as the bare temperature 
T — > O"*""^ The above definition of the f2k further implies 
that conventional scaling holds for all T = correlations, 
i.e. u^'^T^o ~ for > 1, as we now show. 



1. zero temperature limit: evading dimensional reduction 

Let us explore this further, one of the aim being to 
obtain zero temperature correlation functions (as well as 
T > ones). The coefficients f2p ((50() have an interest- 
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ing physical meaning, linked to the so called "dimensional 
reduction" (DR) property. They are the rescaled cumu- 
lants of the so called random force, precisely: 

:.(2p) 



F(0, ri)..F(0, rfe) = S',,'\{0, 0)S'{r,,- ■ ■ , r^m 

The model Q), as other disordered systems such as ran- 
dom field spin models, has the amazing DR exact prop- 
erty that if one computes any observable (e.g. correla- 
tion functions) at T = in perturbation theory in ana- 
lytic disorder cumulants, the result is equal to all orders 
in perturbation theory to the result in the simple model 
with only a random force, introduced by Larkin to model 
physics at short scales: 



rfl 



d^r 



\Vu[ 



u{r) •i^(0,r) 



which is simply, of course (in Fourier) 
uiq)^FiO,q)/q^ 



(59) 



(60) 



we will see that the f2p do indeed determine the corre- 
lations at T = (connected with respect to disorder), 
simply as expected from (16011 : 

u(gi)..M(g2p) V=o = (2TrfS'^{Y^ g^)CT=o((Ji, • • • , q2p) 

i 

CT=o{qi, ■ ■ • ,92p) = qi'-q2pS[Tii{0, • • • ,0), (61) 

except that from the TBLA these cumulants are I de- 
pendent and flow under RG, yielding a non trivial result, 
different from the '"naive DR estimate" which assumes 
them constant and is incorrect. Thus the structure l|47|l 
has nice properties: it does obey correctly the DR prop- 
erty, since setting naively T = by simply removing the 
term proportional to T; in (|47|l does indeed leave us only 
with a random force. At the same time, though, the TBL 
provides a novel and non-trivial mechanism to evade the 
"naive DR estimate" since it gives a non-trivial renormal- 
ization of the random force cumulants. It is particularly 
interesting since only two other mechanisms have been 
proposed previously to evade DR, namely (i) replica sym- 
metry breaking, only convincingly demonstrated in mean 
field models (ii) the cusp and non analyticity of the dis- 
order cumulants in the FRG, demonstrated in the statics 
only to one loop (although up to two loops at depinning). 
Here we have constructed another one, as apparent from 
the structure H47|l . 

Specifically the zero temperature correlations can be 
computed by restoring the I dependence, using from 122() 
that ^If.l/O, • • • , 0) = A;,-2PAf-2pe+2pC^(2p)^^^(0^ . . . ^ 0) 

and the TBLA property ^ that ^{^^^(O, • • • ,0) con- 
verges to its fixed point value e^^X^^f2p- We can then 
simply estimate the correlation when e.g. all momenta 
have similar magnitude gf 
one finds from H61|l : 



q^. Setting I = ln(g/Ao), 



Ct=o('?i 



>92p)|g2 



^p/2pe ^q 



2p„-2p(d-l+C) 



with bp = A^^'^^x^P. Taking into account the Fourier 
transform, this is the form expected from the scaling 
i.e. real space correlations ...itr-2p'^ o ^ i2pC 
[L the common scale). The prefactors of the correlations 
however, and thus the full probability distribution, are 
determined by the constants /2p. It is thus likely (but at 
this stage not established) that they should be universal 
(for long range and periodic disorder, and up to a com- 
mon non universal scale for short range disorder) . Similar 
correlations were computed recently for depinning'^^ . It 
is worth noting that, as was the case for i?"(0), the terms 
involving the higher cumulants of the random force f2p 
do not feed back into the renormalization of the nonlin- 
ear part of pinning force (and thus do not appear in the 
trln), as was the case for i?"(0). This is also to be ex- 
pected on physical grounds, since they only result in the 
shift 159() and can thus be eliminated by a translation. 



2. non-zero T, droplets 

In addition to the terms containing /2p which describe 
the zero temperature limit, (|47|l also contains a term pro- 
portional to Ti which describes the first correction pro- 
portional to T in the correlation functions. Remarkably, 
its content is precisely the one of the droplet picture. Let 
us indicate here how it works, remaining very schematic 
for now. A more precise and detailed formulation will be 
given in Section llVI 

Let us consider an arbitrary correlation amongst repli- 
cas (connected with respect to disorder). Since same 
replicas correspond to same thermal average, it corre- 
sponds to some disorder average of some product of ther- 
mal averages (itself connected^), and we write schemat- 
ically: 

{u-^^l ..uiii), = {2^Y5\Y. qi)c^,l::Zi (63) 

i 

= {u..u)..{u..u). (64) 
Lowest order perturbation theory gives 

rp2p 



c: 



gi,--92p 2 2 

yi--y2p 



y(2p)(0)^ 



(65) 



where here V'^^^^O) denotes a 2p-th derivative at u = 
0, the exact replica index structure of which we ignore 
here. Substituting the rescaled form for the characteristic 
function yields 1^''^^(0) = A^Af+'P'^y/'P^ (0). Inserting 
now the TBLA (|tZ|) gives 



f^ai..a2p _ 1^ 2p 



l^d+2pC-2pe 



2 2 -U2p+T^v(^P\0)) 



91 -95 



ip 



(66) 



Thus we find that at T 0, as expected: 



(62) 



{u..u)..{u..u) u{qi)..u{q2p) \t=o (67) 
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obtained above H61I62() (since at T = the correlations 
become independent of the rephca index - in the absence 
of RSB). The contribution of thermal excitations, to low- 
est order takes the form 

c-ct=o = iTqy^p-\-^p^''-^+<h;v^^p\o),m 

with bp — A^^'^^x^P^^. This is exactly the scaling ex- 
pected from the droplet picture at low temperature. In 
real space: 



C-C' 



T=0 



(69) 



with C = (ur-^..) .■■{■■Ur2 ) • The complicated set of TBL 
fixed point coefficients thus encode for the distribution of 
low energy excitation. This equivalence will be further 
explored below. 

Another question is what fixes the coefficients f2p ? 
Let us examine the FRG equation ljA15|) for the fourth 
cumulant (see Appendix^ . We need to take four deriva- 
tives at zero. Since R"(u) starts at m^, S'J^q as u^, s[^q as 
u^, only one term in addition to the rescaling term con- 
tributes. The analysis is identical for any 2p-th cumulant 
and one finds 



disi'^^io, ..0) = {d- 2pi0 - c))5iT(o, ..0) 



(2p), 



.0), 



(70) 



a simple generalization of the relation (|37|l for the sec- 
ond cumulant. Thus the f2p are determined from the 
hierarchy, and from the function v itself. 

A further hypothesis is that these f2p match the be- 
havior of the outer solution, in particular. 



f2p = lim ]im S[f\{m, • • • , U2p) = 5irU0+, • • • , 0^ 



(71) 



Here we have introduced the notation 0+ for the small ar- 
gument limit of the outer solution. Given the existence of 
the TBL, this is not a trivial statement. It is a particular 
case of a more general question of how (and which) zero 
temperature quantities match to T ^ objects within 
the TBL. The appealing result of the above postulate is 
that the f2p coefficients are truly properties of the sys- 
tem at zero temperature, despite being defined from H5U|) 
deep inside the TBL region (u = rather than u = 0+). 
This is, however, only an hypothesis. The structure and 
consistency of matching is explored and verified in some 
detail in d = in Section Ivl 

If one attempts to generalize the TBLA to match fur- 
ther and further (nonanalytic) derivatives of the outer 
solution, it seems natural to expect that: 



T, 



v{u) 



T, 



V2{u) 



(72) 



Another consequence of the TBL ansatz is that the 
scaling of moments of the pinning force derivatives 

{S'i^^}.n2k{'^i • • • ,0) for fc = li2i- • ■) '^ot modified from 
the naive results of Refs. liollsil Thus, in particular, the 
interpretation of the boundary layer as representing ther- 
mal activation, and moreover the equality of the bar- 
rier exponent with 6 remains unchanged. It is somewhat 
less clear, but we think likely, that the general foTin of 
the velocity-force characteristio^2i^ (in the low velocity 
limit) should be unmodified, as this relies only on the 
scaling assumption. We leave this question open to fu- 
ture analysis (see also the parallel work on the dynamics 

However, although the scaling of moments is pre- 
served, the prefactors are certainly modified by inclu- 
sion of higher cumulants. Indeed, we have shown that 
the ^(3)(-o^o,u) term in ^ modifies R'^^\fS) (and all 
higher derivatives) if the full hierarchy is solved self- 
consistently. This value appears explicitly in calculations 
of susceptibility fluctuations^ and in the velocity-force 
characteristic^^, so these results may be modified by the 
non-Gaussian distribution encoded in the higher-replica 
terms. Note however that s in H4U|) starts as and thus 
it cannot change the result r""(0) = 1, thus H37(l remains 
true. It would obviously be desirable to solve the hier- 
archy encoded in (|54|l to obtain these and other physical 
quantities. This is an important open problem. 



III. EXACT RG AND GENERAL IDENTITIES 
IN ANY DIMENSION 

The Wilson one loop approach was helpful in explor- 
ing the new physics of the thermal boundary layer. Since 
the non perturbative nature of the problem has been un- 
veiled, the validity of the approximations implicit in the 
one loop Wilson scheme, which is known to work in per- 
turbative situations, are here difficult to assess. In par- 
ticular the assumption of uncorrelated space points does 
not remain correct, and a more controlled and general 
technique is required to perform a detailed analysis. We 
therefore turn now to a more powerful Exact RG (ERG) 
method. Within the ERG, we will still not obtain a com- 
plete solution of the TBL hierarchy, but we will obtain 
some concrete results. For instance, in d = 0, up to a 
change in some coefficients (cp defined above) the one 
loop FRG equations are exact. This is very interesting 
both as a new way to study d — Q but even more so as 
a way to confirm the structure found here in an exactly 
solvable case, since in some cases exact solutions exist in 
d = 0. This paves the way for future ERG studies in 
d > 0, as discussed in Section Ivll 

In this section we review the ERG method, and obtain 
useful exact identities valid in any dimension. We will 
turn to its specific application in d = in the following 
section. As a general method, the ERG was introduced 
long ago^SiS and much studied since, but not so often 
for perturbative calculations, and rarely in the context of 
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disordered systems. In recent studies however a system- 
atic multilocal expansion scheme was developed^Si^ and 
applied to perturbative study of pinned systems and the 
FRG. We wiU find useful the version developed directly 
on the effective action functional*^®, and will borrow from 
it. In parallel, we will also derive simpler, but powerful, 
exact relations, directly for correlation functions using 
from scale invariance and STS. These will teach us much 
about the physics associated with the thermal boundary 
layer. 



effective action 



where the "proper vertices" T^'^i^^^lp^^" are the 2p-th 
functional derivative of T[u] in m = and are the sums 
of all IPI graphs with 2p external legs (note here 2p does 
not refer to the number of replicas, but to the number 
of external legs). A different expansion is the multilocal 
expansion: 



TM - To + 2^ / (q^ + m^)u-^u% - VM, (78) 



V[u]= IV{ur)+ I V'-^Hur,Ur',r-r') + 



(79) 



We are interested in calculating correlation functions 
in the replica theory, typically: 



' m^), where we have decom- 



{<<<ui) 



(73) 



as we recall correlations with an odd number of u fields 
vanish. A useful object in doing so is the effective action 
functional r[u]. We recall its definition as a Legendre 
transform: 



dh^W[h] 
W[h] = \nZ[h] 



W[h] 



(74) 

(75) 
(76) 



where Z[h] was defined in ©. It thus has a physical 
interpretation in terms of a free energy with imposed 
expectation values for the fields (e.g. probability distri- 
bution of the order parameter). 

Let us consider the perturbation theory of the repli- 
cated problem defined by Z[0] in (jS)), the free propagator 
being: 



T 



T^ab 



(77) 



and the interaction vertices come from the disorder cu- 
mulants. We have added a small mass which confines the 
manifold in a finite region in u, and provides an infrared 
cutoff. Other IR cutoff, such as finite size L, or momen- 
tum IR cutoff A' < A can also be considered (the UV 
cutoff is still denoted A). 

We recall that W[h] is the generating function of con- 
nected correlations, noted {■■)c, which are sums of con- 
nected graphs. r[u] is the generating function of one 
particle irreducible (IPI) graphs (which cannot be dis- 
connected by cutting one propagator) and thus contain 
loops. The interest of the effective action functional is 
that since it resums the loops, correlations are simply 
obtained from it as sums of tree diagrams. 

There are two useful expansions of the functional T [u] , 
the polynomial expansion, schematically; 



T[u] = r[o] 



2 ^ 

ab 



y{2)ab a b 



1 

4!^ 



r(4), 



with To = -inL'^X^ln(g2 

posed the interaction part, noted — V[w], into a local part 
defined by the function V{u) of a replica vector, a bilo- 
cal part V^^'>{ur,Ur',r — r') defined as a function of two 
replica vectors and one space argument, with by defini- 
tion (u, V, r) = 0, trilocal etc.. Each of these terms 
can then also be expanded in sums over different num- 
bers of replicas, as was done in 1^. Because of STS, 
V, .. contain only sums over k > 2 replicas (the 
single replica part is not corrected) and thus represent 
the "renormalized disorder" {V encodes all cumulants of 
the "renormalized random potential" at the same space 
points, y*-^^ at two different space points, etc.). 

From this Section on, we will make a convenient change 
of notation. We will denote the original action (i.e. the 
bare model) studied here as: 



m2)M>1, 



VnlUr 



(80) 



where the index denotes the bare model. The func- 
tion V{u) thus defines the "renormalized disorder" and 
is expanded as in Q in terms of the "renormalized cu- 
mulants" R,S^''^ while Vo{u) represent the bare disorder, 

(k) 

and bare cumulants Ro,Sq . A similar expansion for the 
bilocal terms, and higher, will be considered below. 



B. exact RG 

The aim is to compute T[u], e.g. the renormalized 
cumulants, as a function of the IR cutoff since some kind 
of fixed point behaviour with universality and scaling is 
expected in that limit. Remarkably, the effective action 
satisfies the exact functional equation: 

J X „ 



(81) 



Gxa,yb = (-V^ -I- m )5ab5xy - ^T""^ 



6u%6u^y 



where the inversion is both in replica and in internal 
space. This equation is supplemented by an initial con- 
dition, namely in the absence of any fiuctuation (with 
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all loops suppressed) the action and eflFective action are 
equal. For instance, if one takes m = oo, or A' = A, in 
that case: 



T[u]^S[u] , V[u] = / Vo[u. 



(82) 



i.e. the renormalized disorder is equal to the bare one 
{V = Vb, = 0, ..). Thus solving the Exact RG 

equation H81|l from the initial condition l|82(l allows, in 
principle to compute T[u\ in the small m limit. 



C. correlations from the effective action 

Once the effective action functional T[u] is known, any 
connected correlation function can be computed, since 
it is the sum of all connected tree diagrams made with 
the proper vertices. This means that it can be computed 
from expanding 

3 — r[u]\ tree 



{uiU2..U2p)c = {uiU2--U2pe '"')^ 



(83) 



(indices here mean both space and replicas) and using 
Wick's theorem, but (as indicated by the "tree" super- 
script) keeping only the connected tree diagrams. The 
two point function is thus: 

r<5afc + T^^^r(2)(g), (85) 



T- 



(g2 



where r'^^)(g) is the off-diagonal disorder at finite mo- 
mentum. At zero momentum it is exactly given by the 
local part V of the renormalized disorder 



r(2)(g = o) = -i?"(o) 



(86) 



as one easily sees that only two legs entering the S^*^-* 
vertex leaves k — 2 free replica sums, i.e. a factor n'^"^. 
However V'^^^q = 0)— r^^^(O) at finite momentum is given 
by the two replica part of the bilocal part V^^^ and higher. 
These can be computed order by order in the local part 
(useful if this one is of order e). In perturbation theory 
it contains all IPI graphs with a flowing momentum. It 
is thus in general simpler to compute only correlations 
at g = in presence of a mass since then one needs only 
the effective action at g = 0, i.e the local part of the 
renormalized disorder. 

The higher correlations can be computed in a similar 
way. For instance, schematically one has for a four point 
function: 

{U1U2U3U4) — {{U1U2) {ui,U4) + 2 perm) 

{uiUy) {u2Uy) {U'iUv) {uiUy) 



xF 



l',2',3',4' 

(4) 

l',2',3',4'- 



(87) 



More generally in expanding (|83|) there is thus a (F^^^) ^ 
on each line of these tree diagrams and the the vertices 
can be only F^^rf^),... 



Analyzing (and even more so, solving) the ERG equa- 
tion is not a trivial task since one has to deal with a 
functional. One way to analyze it in d > is the mul- 
tilocal expansion scheme (|79f) . which we will explore in 
Section IVH We first study it in d = since there it re- 
duces to a flow equation for a function. Interestingly this 
flow equation is very similar to the Wilson FRG equation 
studied in the previous Section. We emphasize though, 
that the relation between the cumulants i?, S'^^^ (i.e. the 
local part of the effective action) and the correlation func- 
tions at zero momentum are identical in any dimension. 
It is thus particularly interesting to study the model with 
a mass, i.e. a confining well for the interface, and study 
the correlations of its center of mass Ucdm since many of 
the d = formulae then apply. These exhibit scaling with 
m and amplitudes which have good limits as d — > (this 
is particularly interesting for the random field interface). 



D. identities from exact RG and STS 

Before we do so we first establish directly on the corre- 
lation functions two sets of simple and general identities. 

The first is an exact RG identity. The exact RG can 
indeed also be formulated directly on correlation func- 
tions. Usually it leads to complicated hierarchies whose 
meaning is unclear. Here it yields interesting constraints 
on the structure of the low energy excitations. Let us 
consider the average of any observable 0[u] which does 
not contain the mass m explicitly: 



(OM) 



where S[u] defined in (|8U|) has a simple explicit m depen- 
dence. Taking a derivative yields: 



<0[u]){ulul)) 

f 



m 



(89) 
(90) 



where the last equality holds for n = 0. This identity can 
be used, e.g. to generate relations between correlation 
with 2p and 2p + 2 u fields. 

The second is the consequence of STS. Let us consider 
again H88(l and in the numerator let us shift integration 
by — > u° -t- , this yields: 

(OH) = {0[u + v]e-^L-^^''^"-^i^J2.<) 
xe"^"/x('^^-)'+"'"' (91) 
To linear order this gives: 

T{J2du,0[u])=m'{0{u]J2ui) 
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-{0[u]J2wlul) 



(92) 



When applied to odd powers of u it generates identities 
between even ones. 

The consequences of these identities will be explored 
in the following sections. 



IV. ZERO DIMENSION: ERG, 
CORRELATIONS, DROPLET PROBABILITIES 
AND THERMAL BOUNDARY LAYER 

In this section, we consider the special case of zero (in- 
ternal) dimension of the manifold (d = 0). For this simple 
limit, many exact results are known by other methods, 
providing useful tests of our very general TBL Ansatz. 
Moreover, since the phenomenological droplet theory is 
known to be correct in d = 0, this limit provides a means 
of further understanding how droplet physics emerges 
from the TBLA. In particular, we will consider the con- 
nection between the TBLA and droplets in detail for low 
moment correlators, including tests of universal ampli- 
tude ratios with exactly solvable limits. Beyond the phys- 
ical connection of the TBLA to droplets, we will also use 
the d = model to understand the matching problem, 
more specifically to study the emergence of a well-defined 
but non-analytic zero temperature RG /3-function beyond 
one loop. The form of this T = /3-function was postu- 
lated in R,ef.l47t but can be confirmed from first principles 
only by a matching calculation such as considered here 
for d = 0. 

A manifold with zero internal dimensions corresponds 
to studying a point particle in N external dimensions. 
Specifically, specializing here to the case = 1, we con- 
sider the partition function of a particle in one (external) 
dimension: 



Zw — 



due' 



[^m^u^ + W{u)]/T 



(93) 



The particle is confined by a quadratic well and feels a 
random potential W{u), often chosen gaussian. For a 
glass phase to exist (as to ^ 0), the correlations should 
be long range: 



{W{ui) - W{U2W = 

= 2(i?o(0)-i?o(w)) 



K{u) 



(94) 



with a > 0. This can be seen as the d — Q limit of an in- 
terface model with long range disorder and a similar mass 
term. The following scaling relations are then obtained 
by replacing L 1/m: 



u ^ m , W ' 
6 ^ 2a/(2-a), 
C = 2/(2-a), 



(95) 
(96) 
(97) 



using the STS relation 61 = 2(C - 1). For a > (6* > 0) 
one expects a single ground state to dominate and the 



droplet picture to be exact. In the marginal case 9 = 
there can be a transition and a glass phase with a Gibbs 
measure dominated by a few sites Ua, the physics of the 
Derrida's random energy model and RSB. This was found 
for gaussian disorder with logarithmic correlations a = 0, 
K{u) ^ Inu, but also holds for uncorrelated disorder, 
provided the distribution has exponential tails. 

The case where W{u) is a Brownian process, the so 
called "toy model" with a = 1/2, has been much studied. 
It corresponds to the d = limit of an interface in random 
field disorder, with: 



C = e/3 = 4/3 , 61== 2/3. 



(98) 



In this case, powerful real space RG methods, or path 
integral techniques, allow one to obtain the full solution 
for the zero temperature glass fixed point and compute 
virtually any universal observable in statics (ground state 
and low temperature behaviour) and dynamics (Sinai 
model)^ii^2iS . We will use these exact results for 
a = 1/2 to test our methods. 



A. exact FRG in d = 

In rf = since there is no space, T[u] is purely local, 
V[u] = V{u) and thus the ERG equation becomes 

T 

-mdmV{u) = n- Tr[5cd j9„^(9„^ V"(w)]~^ 

TO^ 

(99) 

It simply describes the flow of a function of the replica 
vector, the characteristic function V{u) of the renormal- 
ized disorder given by © in terms of its cumulants (the 
trace and inversion refer to replica indices). One notes 
that it is almost identical to the Wilson FRG equation 
(|24l 1 ■ up to the coefficients Cp = 1 here. Here, however, 
in d = 0, it is exact. There is thus no need for new cal- 
culations to obtain the FRG equations for all cumulants, 
which are the same as obtained in Appendix^ Here we 
define the rescaled temperature and cumulants as: 



T = 2Tto'', 

R{u) = ^m-'^^R{um.'^), 



(100) 
(101) 



5«(wi, ■■■,uk) = (uiTO^, • • • , UfcTO^). 

(102) 

Now V{u) = V{um'^) instead of while V{u) is still 
given by (|49() in terms of the rescaled cumulants. Then 
we find that the ERG flow equations in d = for R, 5*^^' 
are exactly the same as (|25l , I26|l setting 7 = 3/4 and 
di —mdrn (and for the fourth cumulant S^'^'^ given in 
Appendix with, in addition 7' = 1/2). The schematic 
form of the FRG equations for all cumulants (|41|l . and 
the TBL form H57|) are also valid here. The detailed co- 
efficients can be retrieved from Appendix IXI 
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Following the analysis of Section^ we also consider 
the TBL form for the rescaled solution of the ERG equa- 
tion: 



p>l ai,..,a2p 



■Ua., + ^viu) + OiTn 



(103) 



with the same definition H50|l for the f2p and for the TBL 
cumulants. The only difference is that here we define 
{ex^ = = -(2 - 9)R*"{0). We expect a fixed point 
value for i?"(0) but here, in d = 0, e = 4 is not a small 
parameter any more and there is no longer an obvious 
relation between ^*"(0) and the R*"'{0+) in the T = 
non analytic theory. In fact in l|25|) we expect the term 
proportional to S^^^ to be important everywhere (while 
it could be neglected O(e^) outside of the TBL). The 
exact TBL equations in d = read: 



p>l ai,..,a2p 



sraf) 



-dadbv{u) 



^ 0, 
(104) 

with 0:2/2 = 1, X2p = —'2p{9 — C) for p > I. 

So we now explore the consequences of the ERG equa- 
tion on the correlation functions in d = 



B. comparison with droplet predictions 

In order to compare with the droplet theory we first 
discuss briefly how to calculate physical quantities in a 
low temperature expansion. One assumes that the be- 
haviour up to order T can be obtained by considering no 
more than two quasidegenerate minima in the rescaled 



We first introduce, in terms of the rescaled variables, the 
probability distribution P{xi) of the position xi of the 
absolute minimum, normalized as J dxiP{xi) = 1. The 
probability density (conditional on the global minimum 
being xi) that the secondary minimum is located at po- 
sition X2 with rescaled energy above the ground state 
e — Em^ , is denoted P{xi, X2,£). Thus one has 

y"dx2^ deP(xi,X2,e) = P(a;i). (105) 



We define the "droplet probability density" 
D{xi,X2) = P{xi,X2, e = 0). 



(106) 



Note the symmetry D{xi,X2) — D{x2,xi). This is the 
probability density that the absolute minimum is (al- 
most) degenerate, and does not obey a normalization 
condition. In the droplet picture it is finite and non- 
zeroS Let us illustrate the way to construct the low 
temperature expansion (see e.g.'*- for systematics in the 
toy model). Let us consider a generic disorder average 
of product of thermal averages of any observables of the 
rescaled position. One starts with restricting to two wells 
(more generally N): 

— — r + 00 !■ 

W_{Oi{x)) = J dE J dxidx2Pixi,X2,e) 
X UipO.ixi) + (1 - p)0.(a;2)) (107) 

(108) 



1 



P 



1 + e-'^-y^ 



landscape separated by a distance x = um 



0(1). 



p and \ — p being the (Gibbs) thermal occupation prob- 
abilities of the two wells. Then one adds and substracts 
the T = part, and introduces w = em~^ /T: 



X{{0,{x))= dxiP{xi)\{0,{xi) 

i i 

p + cc p ^ 

+Tm^ J dw J dxidx2P(xi, X2,Tm^w){Y\_ 
1 

P = ■ 

1 + e-^ 

We have assumed that the restriction to two wells is valid, 
that the thermal width of the packet in each well is small 
compared to the distance between the wells. This width 



(109) 

>,ixi) -f (1 - p)0,(x2)) -llO,{xi)) + OiT^m^'), (110) 

i 

(111) 

I 

results in subleading corrections^. 

For small m one gets the lowest order droplet result: 



{0{x)) - / dxiP(xi)Y[0,{xi) 



(112) 



r + oo r 

+Tm' j dw j dxidx2D{xi,X2)([\{pO,{x^) + (l-p)(),{x2))~\[0,{xi)) + 0(T' 



17 

(113) 



provided the integral converges, but it usually does since 
the one well part has been substracted. Generalization 
to N wells is studied in^. 

A first interesting constraint on these droplet functions 
is obtained from STS. The identity from Section [ffll spe- 
cialized to d = gives in replicas 

T{Y^du^O[u])^m'{0[u]Y^uf) (114) 

/ 

where 0{u) is any odd function. For the particular case 
where 0{ua) is a function of a single replica, one directly 
obtains (for n = 0) 

T{duO[u]) = m'[{0[u]u) - {0[u]){u)). (115) 

To 0{T), the expectation value on the left hand side can 
be evaluated directly at T = 0, while the right hand side 
is determined from averaging with respect to D{xi,X2). 
This gives in rescaled variables 

J dxiO'{xi)P{xi) = (116) 

2p{l ~ P) J dxidx2 [xi — X2)D{xi,X2)0{xi), 
I 



for any odd 0{xi). Thus one finds the implication of STS 
within the droplet picture: 

— P'{xi)= J dX2{xi — X2)D{xi,X2)- (117) 

As discussed in the Appendix, this is the most general re- 
lation between droplet probabilities that one can extract 
from STS. Upon multiplication by odd functions (e.g. 
powers) of a;i and integration it generates an infinite se- 
ries of identities between zero temperature observables 
and finite temperatures ones. 

One may also obtain an RG equation for the droplet 
probabilities. One considers the energy difference 

2 

E^^-{ul-ul) + V{u2)~Viu,) (118) 

between the two minima. Upon an infinitesimal variation 
of m, the total derivative comes only from the explicit de- 
pedence on m, since the contribution from the implicit 
variation of ui,U2 with m vanishes using the minimiza- 
tion condition. Thus, assuming no creation or annihi- 
lation of additional minima, one obtains taking into ac- 
count rescaling 



-mdmP{xi,X2,€)= [2C + e + Cxrd,^ + eede - {xl - xl)d,]Pixi,X2,e). (119) 

I 



for e > 0.^^ Integrating this relation over e and X2, and 
assuming mdmP{xi) = gives 

-C(P(xi)-f a;iP'(xi)) = y dx2{x\-x\)P){xx,X2). 

(120) 

This relation p2()|l obtained within the droplet picture 
from the above RG equation, can also be obtained by 
a similar method to that used above for STS, from the 
more general ERG relation (not itself assuming droplets) 

-ma™(OM) = ^ Y.^O\u\uUi), (121) 
/ 

a particular case of the ERG equations, written directly 
on the correlation functions given in Section IlIII It also 
generates another infinite set of identities between zero 
temperature observables and finite temperature ones, 
also discussed in the Appendix. 

We now compare the droplet predictions for the corre- 
lation functions to that of the ERG. We begin with those 



of the droplet model. 

1. correlation functions from the droplet theory 

At zero temperature we are concerned with the sam- 
ple to sample fluctuations of the location of the global 
minimum. The simplest such quantity is the variance: 

^T=o = m-^H4)p, (122) 
where we defined the average 

{0{xi))p = j dxiP{xi)0{xi). (123) 

The non-Gaussian nature of the distribution P{xi) is 
probed by the four-point correlator 

'^-i^'' ^m-^^{{x\)p~i{xl)l). (124) 

Other correlation functions probe the thermal fluctu- 
ations of the particle. In particular one can deflne the 
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sample dependent susceptibility Xs = {u^) — (u)^ and its 
moments. Although it is negligible in a typical disorder 
configuration, it becomes large when there are two quasi- 
degenerate minima with Xs = m~^^p(l — p){xi — ■ 
Its disorder average is the second cumulant of thermal 
fluctuations 



T 



{u^) — (w)^ ~ —m^m I dxidx2 [xi — X2)^D{xi, X2) 



(125) 



We used (jTU^ and p(l - p) = i (see Appendix IdJ. The 
symbol ~ (here and throughout this section) indicates 
that the first term on the right hand side of the first 
line is intended to give correctly the amplitude of the T- 
linear term. In fact, as is well known, this correlator is 
fully constrained by STS to give exactly the free propa- 
gator result. Within the droplet calculation this can be 
obtained from (|117|l upon multiplication by xi and inte- 
gration. For the present model (confined by a harmonic 
well) it is exact at any T (i.e. all higher order corrections 
cancel) . 

At quartic order one may define several thermal corre- 
lation functions. One of them is the second moment of 



)2)2 ^ ^m«-4C J ^^^^^^ (^^ _ 2.2)4^(^^ ,^^) 



(see Appendix IdI for details). There is one other combi- 
nation 



T 



{{u — («))■*) ~ ^'^^ J dxidx2 {Xi — X2YD{xi,X2) 

from the droplet calculation. These two combinations 
are however, related by STS which implies that: 



{{u-{u)Y)^?,{{u^)-{uyY 



(126) 



which again is an exact relation valid for all T (see Ap- 
pendix ^ and"*^ ) . Although it cannot be derived from 
(|117|l . the droplet theory is clearly compatible with it. It 
in fact entails a propertyS^ of the various moments of p. 

Remarkably the RG, combined with STS, allows to re- 
late these quartic T > observables to quadratic ones at 
T = 0. Here we can see it from the droplet theory. Indeed 
multiplying pi7|l by Sxf , (|120|l by (!>x\ and subtracting 
we obtain: 



^ J dxidx2{xi - X2)'^D{xi,X2) 



2- 



(127) 



A fixed point has been assumed. Similar relations exist 
at the fixed point between 2n and 2n — 2 observables for 
any n > 2. 

Finally, it is interesting to look for a combination which 
should have contributions only from three wells. Within 
the droplet theory the coefficient of T should vanish and 



thus this correlation can be used to test droplet theory. 
The simplest one occurs at sixth order and reads: 



6 



(128) 



= 2{U^){U-^){U) - (m2)3 + (m4)((u2) - (u)2) _ (y3)2 

this is the lowest order necessary to see 3 well physics, 
since the combination {(ua — UbY(ua — Uc){ub — Uc)) is 
zero by symmetrization (here a, 6, c denote any triplet 
of distinct replicas). Thus the prediction of the droplet 
theory is that: 



(129) 



correlation functions from the ERG 



In this section we give the relation between a set of 
correlation functions of u and the m-dependent cumu- 
lants i?,5'('=), using the properties of the effective action 
summarized above in Section lTTll Details can be found in 
the Appendix 10 We stress, as discussed above, that the 
relations given here between the R^S'^^^ and the correla- 
tions are exact, not only in d = 0, but in any dimension 
(if one considers zero momentum u — > Uq=Oj center 
of mass). The rescaling and, of course, the fixed point 
values, does however depend on d (see Sections ^ and 
IVI|) . Next, we insert the expected asymptotic flow of the 
cumulant for m — *■ using the TBL ansatz and discuss 
consequences. The resulting relations are then also ex- 
pected to be "exact" but with the following disclaimers 
(i) they are subjected to the TBL assumptions and (ii) 
there are accurate up to subleading corrections in m. For 
them we will use the symbol ~ and will, here, reserve the 
symbol = for relations which are exact by construction 



1. two point function 



We start with the two point function: 



1 



1 



uV) = ^,6-— i?"(o) 



(130) 



which contains a zero temperature part, and a thermal 
connected part which form is dictated by STS. Higher 
correlations will have similar properties and we will ex- 
plicitly write them up to six-th order. To classify them 
it is convenient, from now on in this Section, to use the 
convention that different replica indices mean non equal 
replicas. Then one has 



{uI) = {u^)=t\-—^R"{0) 

77?/ 777^ 
1 



{UaUb) = (7i)2 = jR"{0) 



(131) 
(132) 



The disorder cancels between the two lines and the sec- 
ond cumulant of thermal fluctuations is thus exactly 
given by (|125|l above. 
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The TELA gives, using Eq. fTUT|l . 



(133) 

This gives the second moment of the position of the ab- 
solute minimum. 



2. four point function 



The calculation of the four point function requires only 
one r^'*) vertex. Since S'^^^ starts as only R and 5*^'*^ 
can give a contribution. The calculation is performed in 
Appendix lUl There are five possible replica monomials. 
Their connected expectation values read: 



III 

{ulul)^^^R^'\0) + ^^ 
/ 2 \ Pi 

[UaUbUcjc = — g- 



{UaUbUcUd)c = 



Fa 



(134) 
(135) 
(136) 
(137) 
(138) 



we denote the fourth cumulant of the random force -F4 = 
S'J^{j^(0, 0, 0, 0). To obtain their full expectation values 
(O) = {0)c + {0)disc one must add their disconnected 
parts, given in ljC2|l in terms of the two points averages 
(|132|l . and which contain only T and i?"(0). 

One sees that these five connected correlations depend 
only on two independent quantities ^4 and and R""{Q). 
The first one measures the fourth order cumulant of the 
displacement at zero temperature. Since in that limit all 
replica averages in H134I) ( and in ^(221 ) coincide 



3u2 



.om-4C-l^(tn(0, 0,0,0) (139) 
lb 



in terms of the rescaled quantities, assuming a fixed 
point. Comparing with H124|) and generalizing, one finds 
that the full distribution P{x) of the position of the ab- 
solute minimum, is retrieved from the ERG fixed point 
from its cumulants: 



2p\c _ J_ q(2p) 



,0) 



(140) 



At T > from H134() there is only one independent ad- 
ditional observable (as discussed in the Appendix ICl this 
can be equivalently seen from the STS-ERG relations). 
One finds: 



1 



(141) 



ri?(4)(0) 



The fourth cumulant cancels in this combination and 
i?"(0) also cancels in the disconnected part (see Ap- 
pendixQ. In terms of the rescaled disorder the variance 
of the susceptibility fluctuations becomes: 



y2 



Xs - Xs 



R^-^HO) 



^7^^'^(0). 

4 ^ ' 



(142) 



At this stage assuming a fixed point for M^'>{0) would 
imply conventional thermal scaling for the fiuctuations 
Sxs = 0{T). The TELA allows to evade that, as R'^^\Q) 
blows up as TO ^ with ^('*)(0) ~ x^'r""{Q)/f with f = 
2TmP , which leads instead to much larger fluctuations: 



Xl-xf^ \{2 - e)(-^*"(0))(Tm«)TO-4f,(143) 

using the exact result r""(0) — 1. The scaling coincide 
with the result above, and the amplitude too, given the 
relation (|127|l which exists at the fixed point 



Xt ~ Xs 



T- 



(144) 



It should be stressed that this relation originates from 
the exact ERG identity: 

-mdmR"{Q) = (2 - e)R"{Q) + fiR""{0) (145) 

assuming that R"{0) reaches a limit. Eoth droplets 
and the TELA are consistent with this exact relation. 
Note that similar identities for the 2p-th cumulant where 
shown from the ERG iCni) and result in relations similar 
to H145|) for higher order observables (see Appendix Q. 

To check further consistency between the TEL and 
droplets, we now compute the three well quantity D3 
defined in (|128|l above. The calculation of the six point 
fmictions is performed in the Appendix lUl We find: 

D3 = ^TO"6Cj.3(g^(4)(Q) ^ 2f ^(6)(0) 



+sgi{0,0,0) + 3R('\0f 



(146) 



in rescaled variables. We note that from the TELA 
as are all last three terms, so naively 
D3 ~ T. However this is not so, as we find that these 
three terms actually cancel to leading order when using 
the ERG equation. Indeed using the fourth derivative in 
zero of the FRG equation (|25|l for R, it simplifies to: 

D3 = ^TO-6Cf3(2i?(4)(0) _„,9^i?(4)(0)) (147) 

~ ^(4 _ 0^){-R*"{O))m-^'^f'^ (148) 
32 

which is indeed 0{T^). It is likely that obtaining the 
correct prefactor here requires considering the next order 
correction to the TEL. Eut clearly the TEL solution of 
the ERG knows that D3 is a three well quantity. Finally, 
the third cumulant of the susceptibility is also consid- 
ered in the Appendix and shown to relate to the third 
cumulant TEL function ^222(0, 0,0). 
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D. comparison with the solution of the toy model 

It is interesting that one case of the model studied here, 
the so-called toy model, has been exactly solved. In that 
case: 



-4/3 



V{u) - ? 
C = 4/3 



,1/2 



,4/3 

= 2/3 



(149) 
(150) 
(151) 



Correlations and probability distributions are known, or 
accessible, and one can infer from them some character- 
istics of the fixed point, e.g. the (universal) numerical 
value of ^*"(0) (see below) etc.. 

We can in particular verify that the solution of the toy 
model obtained in Ref.^'^ obeys the exact relations from 
STS and ERG obtained here. Let us start with the simple 
relation between fourth and second moments. With the 
choice K{u) = 2a\u\ at large u in (|94|l . with (7 = 1, one 
has at T = in the limit to ^ 0: 



^"(0) 8/3 



C2 



+ 00 



du Ai'{iu) 



2^/3 /•+°° du -iu 
~3 



27r Ai{iu) 



(152) 
(153) 

1.05423856519 (154) 



This implies: 

-^*"(0) = 4c2 
On the other hand one haa^: 



(u-(u))4=3((m2)_(u)2)2 



/i4 



27/3 



+ 00 



du 



27r Ai{iuy Ai{iu) 



Ai'jiu) 

2 ^ ) 



(155) 

(156) 
(157) 

(158) 



The above relation (|142|l between this quantity and 
^(*)(0) imphes that one must have: 



(159) 



where T = 2Tm^. This is clearly consistent with the 
TBL ansatz. Thus for the above prefactor to be correct 
one needs the non trivial relation: 



M4 



4C2 



(160) 



which we have confirmed numerically. 

One can go further, since the full distribution P{xi) 
and the droplet distribution D{xi,X2) were obtained 
for the toy model. First one notes that if P{xi) and 
D{xi,X2) are solution of H117|) and l|120l) then: 

P{xi) ^ pP{pxi) (161) 
D{xi,X2) ^ p'^D{pxi,px2) (162) 



are also solutions (they correspond to different choices 
for (T and the rescaling factor of u). One member of this 
family (corresponding to the choice p = 1) is: 



Xl) 



P{xi) = 9{xi)9{~xi) 

D{xi,X2) — D{XI,X2)0{X2 

+D{x2,xi)9{xi - X2) 

D{xi,xi +y) = 2g{-xi)d{y)g{xi 



with the specific forms: 

'•+°° dX 
— ( 
, 2tt 



d{y) 



-i\x 



Ai{i\) 



(163) 
(164) 
(165) 
-y) (166) 



(167) 



„ -e-° , ; = y 6-1"=!^ (168) 
27r Ai{i\) ^ ^ ' 

^ ' s— 1 



in terms of the zeroes as of the Airy function. Comparing 
with"*" one further sees that the choice considered above 
corresponds to p = 2-^/'^. 

We check in Appendix IdI that these functions indeed 
satisfy the relations (|117|l and H120|) . An outstanding 
question is whether these relations could allow some an- 
alytical progress in the other cases. 

To conclude this Section we have established a precise 
mapping in d = between droplet probabilities and some 
thermal boundary layer quantities. Questions which de- 
serve further investigation are: (i) how much constraints 
the FRG-STS equations give on the set of all thermal 
boundary layer functions (associated to all disorder cu- 
mulants) (ii) whether D{xi, X2) may contain all the infor- 
mation encoded in this set (iii) what part of these func- 
tions is universal. Answers to these questions in d = 
would certainly be very helpful to understand the general 
case. Before discussing d > 0, we turn to the question 
of how the r = information, namely P(xi), can be 
retrieved from the TBL study. This necessitates consid- 
ering the so-called matching problem. 



V. MATCHING 

In this section, we use the ERG method for d = de- 
veloped in the previous section to consider the matching 
problem in some detail. As discussed in the previous sec- 
tion, most physical quantities are formally determined by 
derivatives of coupling functions in the effective action 
evaluated at zero argument. As such, they are defined 
deep within the TBL. Physically, however, one can divide 
these quantities into those describing zero temperature 
(sample-averaged ground state) properties and alterna- 
tively thermal (T > 0) fluctuations. While it is natural 
that the latter are determined by the TBL, which, as we 
have discussed, is an encapsulation of droplet physics, it 
is surprising for the former zero temperature quantities 
to require a TBL analysis. Thus one is led to expect 
some sort of matching of certain derivatives of coupling 
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functions at strictly zero argument (and T ^ 0+) to cor- 
responding limits as Uab ^ O"*" of T = functions. The 
subtlety in the latter definition is the order in which the 
limit of coinciding points should be taken. It will natu- 
rally emerge from a careful study of the matching. 

Unlike for R{u), all higher cumulants are described 
by functions of more than two separate points. As the 
outer solution for such cumulants is non-analytic when- 
ever any two points are brought nearby, each such non- 
analyticity must be resolved by some boundary layer 
structure. Therefore for each cumulant we will intro- 
duce a set of Partial Boundary Layers (PBLs), one for 
each distinct way of bringing subset (s) of points together. 
We will adopt a transparent notation to specify the dis- 
tinct PBLs. Each PBL consists of a set of q groups, all 
rii points within group i having been brought together. 
Thus for a given k-th cumulant there is one PBL for each 
integer partition of k, = ^- We label the corre- 

sponding PBL functions by the integers Ui in decreasing 
order. For instance s^'^^ is the full TBL introduced previ- 
ously and s'-^^^ is the only PBL associated with the third 
cumulant, and we will denote simply the outer solution 

We will illustrate this procedure explicitly for the third 
and fourth cumulants, making a full consistency check 
for the third and a partial check on the fourth. In the 
end this construction enables and justifies a expansion in 
powers of R of the fixed point equation for the function 
R(u), i.e. the beta function of the model. For comparison 
with the e expansion, we will keep e = 4 — d = 4 factors 
unsimplified. 



A. analysis of third cumulant 

The necessity of understanding matching in some de- 
tail is evident from inspection of the equation H25|) for 
R{u). We have already noted the feedback of the third 
cumulant in this equation, but have not emphasized that 

even for u of 0(1), this feedback term S'Jjq(0, 0, u) is re- 
quired for two arguments coincident with one another. 
Since the outer solution is non-analytic whenever pairs 
of arguments coincide, we should expect some kind of 
boundary layer softening of this singularity to render 
5J^^q(0, 0, u) well defined. Because, however, u remains of 
0(1), this requires a new Partial Boundary Layer (PBL) 
treatment. The tantalizing question is how such a PBL 
analysis can recover a sensible zero temperature limit, 
in which naively SIiq{0,0,u) above is replaced by some 
quantity defined entirely outside all BLs. 

We investigate the PBL for the third cumulant through 
its ERG equation it^ . The PBL will be taken to cor- 
respond to Ui2 = exui2/Ti = 0(1), wi3 — 0(1), and 
Ti <C as usual. We would like to constrain the PBL be- 
havior by assuming that 5ji|)(0,0, u) is 0(r°) (and 0{e^) 
hopefully in the e expansion) , and moreover is quadratic 
in u at small argument. It clearly should be, so as to 



correct the fixed point value of i?"(0) compared to a low- 
est order truncation. A naive scaling ansatz would then 
be S^^\u^'^^) ^ f^Xs'-'^^K''j-i2,ui3). However, it is readily 
seen that this is inconsistent in the above equation, be- 
cause then the thermal term would be much larger than 
all others (i.e. 0{Ti) as compared to 0{Tf )). Instead, 
the only way to consistently obtain the proper behavior 
of S{ii^{Q,Q,u) is to postulate a purely quadratic in ui2 
term with the aforementioned T^x scaling, a fully com- 
plex PBL term appearing only at 0{Tix), i.e. 

5(3)(mi,U2,U3) -T'fx"?20(wi3)+rfs(21)({ii2,Mi3) 
+0{f^), Ui2=0{\) , Wi3=0(l). (169) 

where (p and S are order 0(1) and even functions. The 
function (j) so defined is particularly significant because it 
gives the feedback of the third cumulant into the second: 

S^^i{Q,Q,u)^~2xmu). (170) 

Note that while does not depend on the choice of pa- 
rameterizing using ui2, Mi3 rather than ui2, U32, the func- 
tion s^^^^(ui2, U13) depends on this choice (as well as V'(u) 
defined below). 

In term of the function (j){u) the outer equation for the 
second cumulant thus reads: 

diR{u) = Q ^ {e - 'iC + Cudu)R{u) (171) 

+ \R"{uf - R"{0)R"{u) - 2x30(w). 

Expanded to O(u^) it yields the exact relation: 

R"'(O+)'-x'(l + 2x0"(O+)). (172) 

Here and below we will omit the * to denote fixed point 
quantities. 

If one considers the large {ti2 limit, one may expect 
to match to the outer solution. In particular, taking the 
small M12 limit of the outer solution, we expect (by a 
choice of gauge, see Appendix IE|| : 

^(^'("123) X'ul2cl^{u^3) + x'WuNiuis) + 0{ut,) 
ui2 = 0(l) , ui3 = 0(l). (173) 

Matching the two forms is an extremely strong con- 
straint, and completely determines 4>{u) from the outer 
solution (outside all BLs) - for this reason we have used 
the same function (f>{u) in Eqs. (|169I173|I . This fact 
explains the conceptual dilemma mentioned above: al- 
though (piu) is a quantity defined inside the PBL, it 
is wholly determined by the zero temperature solution. 
Thus the zero temperature limit of Eq. H25|l is well- 
defined. As a check on this matching requirement, we 
will verify that the equation for (t>{u) obtained by anal- 
ysis entirely with the PBL is equivalent to the equation 
for 0(m) obtained from the small 1*12 limit of the outer 
solution. 
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Note that all terms in the series in Eq. H173|) are pow- 
ers of |mi2|, since any odd power of U12 can be converted 
to a higher even power by symmetrization. Thinking 
further about this, there should be a hierarchy of such 
terms, since the cj) term can only capture the 0(u^2) co- 
efficient of the small U12 limit of the outer solution, the 
g(2i) term only the 0(1^121"^) etc. In general the coeffi- 
cient of 0(|ui2|") in the outer solution is matched by the 
large U12 ( ^ 11^2 ) of a generalization of S at order T" in 
the PBL expansion l|169(l . In particular one must have: 

S*^^'(W12,M13) ~ |U12|^'/'("13), W12 ^ 00 (174) 

In both regimes (inner PBL and outer small U12 limits). 



the higher and higher terms are progressively smaller, 
so we may hope to truncate this successfully. Unlike at 
0(u^2) for 4){u), at 0(7/^2), matching to the outer solu- 
tion (i.e. ?A(ui3)) only determines the large U12 limiting 
behavior of s'^^^ (-ui2, M13), which has the full PBL com- 
plexity otherwise. 



We first consider the PBL constraints implied by 
the S*^^^ equation. One finds, applying the ansatz of 
Eq. (|169|l . that all terms are 0{T^) or smaller (we do not, 
as appropriate to d = 0, assume e <C 1 for the moment). 
Keeping only these terms, and for simplicity dropping 
the feeding term from the fourth cumulant, one has 



= 9zS(3)(ki23) = TM2X[{'^^ - 2 - 4C)<^(U13) + CW130'(^13) 
+ 0"(U13)R"(U13) + 30'(7/i3)R"'(ui3)] 
+f2R"(ui3)[r"(«i2) +7(r"(iti2))' +4ii,(0-0,^ii2)] 

+TiX^[[s2q\ui2,U13) - S20^\0,U13)) r" {U12) + S^o^^(ui2,Ul3)] 

+rz'x'[- ^i?2^^''^'("l3) - ^420^(^12, 0,1^13) + 4l00(^*112,"l3)] + 0{ff). (175) 



R"(^13) 

X 



To obtain this equation we have performed a expansion 
in small U12 discarding some "gauge terms" as detailed 
in the Appendix. Note that the dangerous term, the con- 

tribution from to T] 5*200 vanishes being a pure gauge. 
For completeness we have included the feedback from the 
PBL31 and the PBL22 of the fourth cumulant (last line) 
which will be discussed in the next Section, but is not 
crucial for the following discussion. In the third line of 
Eq. 11751 we have grouped a contribution arising from 
the fi[R"]2, R"^no, and [R"]^ terms in Eq. This 
is simplified by eliminating s\^q{ui,ui,U2) using the full 
BL equation for the second cumulant: 

■5no("i; "1: "2) = ]^ul^ - ]^r"{ui2f - r"{ui2) (176) 



The resulting equation can be studied at large W12, us- 
ing the matching constraints Eq. (|174|l and r"{ui2) ~ 
r^|ui2| and the corresponding limit of the fourth cumu- 
lant quantities (discussed below): 



symi2sfio|,({tii2,ui3) ^ -2u\2g^''^^\ui3.) (177) 

Symi2420 («12,0,M13) ~ 2?i22g(22)(y^3) (lyg) 

Collecting all 0{u\2) terms (the dominant ones in this 
limit), one finds 



(2e - 2 - AO^{u) + Cu<l>'{u) + — ^ [2+^^- 2^^'^)' ^ 
+(/."(u)R"(u) + 30'(u)R"'(w) + &xrZi^{u) (179) 

This equation allows to compute 4'{u) in an expansion in powers of R as we will explicitly show below. 



R"'{u) 



Conceptually, one should regard 4>{u) and ili{u) as de- termined by the underlying zero temperature fixed point. 
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since although they obtain in the PBL, they are entirely 
constrained by matching. Thus (j){u) should be obtain- 

I 



where the last term arises from the PBL211 of the fourth 
cumulant. Eq. (|179|l should match the small Ui2 limit of 
the outer equation. Inserting the expansion of Eq. (|173|1 . 

I 



Note that these equations are similar but not identical, 
and in the derivation of Eq. H181|) we used neither the 
thermal terms in Eq. (|26|l nor the BL equation for r{ui2). 
Nevertheless, they are equivalent as can be seen by us- 
ing (|172|) in Eq. (|181|l and the matching requirement 
R"'{0+) ^ xrZ: one directly obtains Eq. ((T7^ . We have 
also used (|i^'45l li^'46|l relating the tilde quantities to un- 

I 



able by considering the outer equation for 5*^^^ : 



(180) 



together with fourth cumulant matching conditions (see 
below) and keeping terms to 0(1/^2)- One finds 



(181) 



tilde ones in the fourth cumulant matching. Note that 
4>{u) starts as (see below). The corresponding term 
is allowed in the PBL, while in the full BL it is 
not (the third cumulant starts as vP). 

The first true PBL quantity which is not fully deter- 
mined by matching is s^^^^ (ui2, U13). Using Eq. (|181ll to 
simplify Eq. H175|) . one finds 



^ (2e - 2 - 6C + Cu^a„J5(3)(ul23) + 3sym{R" {ui2)i-2x^Hui3) - ~S^r^^(uii-i)) 

+7R"(U12)R"(7.13)' - |R"(^12)R"("23)R"(«13)) " 3x Vm^'^") (u^^ , ?/l3) 

I 



(2e-2-4C)0(w) + C"0'(")+7 



R"(w) 



X 



R"'(0^ 



X 



R"'(m) 



+ (0"(ii) - (/)"(0))R"(u) + 30'(u)R"'(w) + 6i?"'(0+)V(u) 
-X(V;(31),(^)_^ 25(31) (^)_^~(22)(^))^0 

I 



= (4o'^(^i2,"i3) - 4o'^(0'"i3)) ^"(^12) + 4o'^(^ii2,«i3) - ^rZul^-^iux^) (182) 

R (^13)/ 1^/ n2 / lll\2~2 \ 
-5 {1-70{t (W12) -(t-oo) "12) 

-^420 ("12>0,M13) +U?25^^^Hwi3) + sfioQ (uil2 , U13) + 2Ui2ff ^^^^ ("w) 

I 

supplemented bv the boundary condition f ll74() . We note for s'^q\u,u). One finds 
that this equation can be formally but explicitly solved 

I 

sif\u,u) - sg'\0,u) = —l-^^{a{u) + biu) + 6u\ZHu) - (7 - i)^(r"(u)2 - (r-)^^^) (183) 
+ 14^0^(^,0,^.) - n'g^'^H^) ~ s^Ml{0,0,~u,u) - 2^25(31) (^)) (184) 



where a and b are arbitrary gauge functions. This solubil- ity is similar to the algebraic solution for r"{u) possible 
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in the full BL for the second cumulant. 

We have thus solved the matching problem up to the 
third cumulant. It is interesting to organise the resulting 
implicit solution for 4>{u) as an expansion in powers of R. 
To lowest order 0{R^), Eq. H181I) becomes: 



{2e - 2 - ACmu) + Cu(j)' {u) 



+7 



R"{u) 



R"'(0^ 



R"'{u) 



(185) 

0. (186) 
I 



Solving this linear equation and inserting feedback H170() 
from S in the R equation we obtain the " beta- function" 
to the same order, i.e. the condition for the fixed point 
in the asymptotic limit Ti — > 0: 



diR{u) = = (e - 4C + Cudu)R{u) 



+ ^R"iuf ~ R"{0)R"{u) + 
2e-2-4C 

^- c 
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fx^{R"{ux) - R"{0)){R"'{ux)^ ~ i?"'(0+)2) + 0(i?4) A > 
T /r ^x^iR"iux) - i?"(0))(i?"'(ua;)2 - R'"{0+f) + 0{R^) A < 



(187) 



(188) 



Both formula match up to a constant when A = 0. Set- 
ting e = 4, 7 = 3/4, A = 2 — 3a = ^^^^p-, the above 
result gives the exact beta function in d = to order 
0{R^). Higher powers can be obtained by studying H179|) 
or (|181|l to next order in the expansion in power of R, 
and considering the feedback from the fourth cumulant 
which is of the same order. Note that at this order, ob- 
tained here rigorously, this result coincide with the naive 
zero temperature procedure. This procedure consists in 
first solving the S equation (|26|l for arbitrary by drop- 
ping terms containing temperature and neglecting the 
RS, and fourth cumulant term, then evaluating the feed- 
back term from S in the R equation (|25|l which at this 
order has an unambiguous limit (see discussion near Eq. 
I31|l . It is unlikely at this stage that the absence of am- 
biguity will exist at higher order. Note furthermore that 
in the limit e — > and C — > the convolution operator 
in lfTK7|l becomes local, e.g. ^x^f{x) /(1)/A, and 
we recover the result of Ref4i, (Eq. 3.43) with X = 27. 
Thus the Wilson value 7 = 1 seems to reproduce, up to a 
factor of 2, the result X — 1 obtained in the e-expansion 
to two loop in Refi^. This is not so bad, since this d ~ Q 
calculation does not involve any spatial integration. It 
raises the hope that a lot of the structure of the pu- 
tative /3-function in a true e-expansion may already be 
contained in this d = model. 

To close the study of the third cumulant we examine 
the matching of the full BL to the partial one. This is 
complicated by the fact that in the full BL regime the 
rescaling term is negligible, while in fact it is not in the 
partial BL and the outer region. Thus as one leaves the 
full BL region, presumably the rescaling term becoming 
important tends to suppress s. We find that to match we 
must have, coming from the full BL: 

S^^\ui23) --uia-oo CiXW?3"l2 



+XUl3p{ui2) + X|ui3|/5(ui2) + 0(1) 

"123 = 0(1). (189) 

while the small argument behaviour from the partial BL 
side should be: 

5^^^(1212, W13) ~«i3^0 Wi3/3(ui2) -I- |mi3|/0(ui2) -t- 0(ui3) 

(/>(ui3) ~«i3-.o ciul^ + 0{ul^) 

U12=0{1) , U13 = 0(1), 

(190) 

with ci = i^"(0). 

B. analysis of fourth cumulant 

We now turn to the structure of matching for the 
fourth cumulant. This is of interest for a variety of phys- 
ical quantities. First, the fourth cumulant of u at zero 
temperature, (|139|l can be rewritten as: 

^- 3(^)2 = ^m-4CxV4 (191) 

and is the lowest order quantity expressing the devia- 
tion of the distribution of u from a Gaussian. Second, 
the beta function to next order (i?'*) requires the fourth 
cumulant feedback into the third. More generally, the 
fourth cumulant is illustrative of the nested structure of 
higher cumulant PBLs. Unlike for lower cumulants there 
is more than one path of successive matchings connecting 
the outer solution to the full TBL. 

Using the notations defined at the beginning 
of this Section, we need to discuss the PBLs : 

s(4)^s(31)^s(22)^g(211)^g(llll) ^ ^(4) fo^ ^j^g f^^^ ^^^-^^^ 

partial BL, and finally, the outer solution. Each of these 



has its zero T piece and their finite T parts. Let us non-trivial order in temperature, in each case: 
enumerate the various boundary layers. One has to first 
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5W(«i234) = hf^UlU2U3Ui + ff X'^ S^^\ui23i) , ^1234 = 0(1) (Ml BL) 

5^(^1234) = T^X Symi23(Ml2ui3)^^^'Hw34) + T^sOl) (iii23 ; "34), Ul23, "34 = 0(1) (PBL 31) 

^(^H"1234) = T'4s(22)(ui2;W34;Ul3), "12, "34, 1^13 - 0(1) (PBL 22) 

5(*H"1234) = f;VG220(211)(y,3,u,4) + f fxsf^") (^12; "13, "u), "12, "13, "14, "34 - 0(1) (PBL 211) 



where ^^^n)^ s^^^i), and s^'^^'> are order 0(1) func- 

tions. The first terms in the fuU BL, PBL31 and PBL211 
are zero temperature components and are anomalously 
large (have a lower power of temperature) compared to 
their associated s functions. These terms are necessary 
for a zero temperature limit to exist. Their form is con- 
strained by the requirements of permutation symmetry 
in each group of boundary layer variables. They have 
the additional crucial property that each contributes to 
its PBL equation at the same order in temperature as the 
corresponding s term, despite the lower power of temper- 
ature in its definition. This is because when two deriva- 
tives act on u BL variables in these terms it produces 
only gauge. The (/)*^2ii) function is a simple generalisa- 
tion of the (j) function introduced in PBL21 for the third 
cumulant. In the PBL31 we have omitted the quadratic 
term {u\2 + u'l'^+U2z)4>^^^^ {'^'ii) since the latter two terms 
in parenthesis are pure gauge and the first can be rewrit- 
ten up to gauge in the form above using that M13 T; 
(see Appendix rE|l . Finally, in PBL22 symmetry in both 
1 2 and 3 4 forbids zero temperature quadratic and 
cubic terms, while the naive quartic one has been sub- 
sumed into 5^22)^ Additional more detailed properties of 
these functions are summarized in the Appendix IFl 



, 092) 

nil 

i 

211 
/ \ 

31 22 

FIG. 1: Nested partial boundary layer structure for the fourth 
cumulant. The limit of two points in a given PBL coinciding 
takes one along the direction of the an arrow in the diagram. 



Particular pairs (inner and outer) of these PBLs are 
connected by taking a single argument (pair of points) 
in the outer PBL small (points close) and reciprocally 
in the inner PBL large (points far). This is illustrated 
in Fig. n The conditions for such asymptotic matching 
sequentially going inwards from the outer solution are 
considered in detail in Appendix IfI This determines the 
form of various small and large argument limits of the 
PBL functions defined above: 
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s(^")(t2i2;Wi3,Ui4) = |ui2|^V'^^"H"l3i'"l4) Ul2 ^ OO 

S(2")(U12;W13,U14) = "13P^^^H^12,"14) + |wi3|o-^^^n^l2,"l4) Ul3 ^ 

s(2")({tl2;Ul3,Ul4) = |W34|0''^^^''(^*12,U13) "34^0 

0(211) (^^3, = Ui4g(22)(^^3) ^3^^0 

S*^^^H^'12,"34;W13) = U34"123*^^''(W13) + l"34|0-(22n^'l2;Ml3) 1X34 ^ OO, {ti2fixed 

S*^^^H^'12,"34,W13) M?2^34 1^220^ (0, 0, 0) +Wi3p(ui2,%4) + |m13|o-(u12,W34) U13 

^('^H«34) = U34^('^^'(0) + iwi4SgnK4)^('^^"(0+) U34->0 

s(^^H^'123;W34) = ■"l3"l23*^^H^t34) + Mi3Sgn({ti3)M^2^^^^Hw34) + Ml3Wl2P'^^H'"34) + |wi3|w?20-'^^H^t34) "13 OO 

s(^^H'"l23,'"34) = i5-('^)(t2i23)M34 + P'-'^''('«123)|'«34| M34 ^ 0, {ti23fixed 

S^''H"1234) = 5XW34Sgn(u34)uf2V'*^^^''(0+) + U34X'5'''^n^'l23) + |m34|P*''H"123) M34 -> OO, {ti23fixed 

s(''Hui234) = Ul3y5(wi2,U34) + |w13|0-(ui2,M34) U13 OO, -1*12, ■lt34fixed 



r 



(193) 



The matching analysis results in some additional condi- 
tions not already implicit above. These are 

V,(31)(u) = ^(31)(^)_^_^(211)^Q^^^^ ^^g^^ 
5^''H«34) = ff(31)(„3,)_ 10(211) (0,^,34), (195) 

The various zero temperature terms can be determined 
sequentially by matching from the outer solution. In par- 
ticular, the small argument behavior of 5(''^ yields (/)^^^^\ 
which in turn yields ip'-^^\ g^^^'^ and g^^^\ which finally 
determine f^. This Larkin term thus matches all the way 
through to give 

/4 = 2^(3i)'(O) = -20(^i")(O,O), (196) 
h - 420' ((^,0,0) = 4g(22)(0) = -20(^1^(0,01197) 

As discussed in Appendix El these derivatives are un- 
ambiguous, and determined by a definite procedure dis- 
cussed therein. 



Unlike these zero temperature quantities, the PBL s 
functions themselves are not determined by matching 
alone. Instead, they obey PBL equations derived from 
the basic ERG fourth cumulant equation (given in (jA15|) ') 
and the PBL forms in ((T^ . The PBL211, PBL22, and 
PBL31 equations are given and analyzed in Appendix IFl 



C. expansion in powers of R 

The PBL analysis of the fourth cumulant enables the 
calculation of the beta function to 0{R^). To do so, 
we require 4>{u) which feeds in p71|l . To obtain it, we 
must solve H181|l perturbatively in R (equivalently x)- 
We write: 

0(u) =0o(u) -I- X01 («) + ■■• (198) 

The zeroth order solution 0o [u] obeys (I185|l , the solution 
of which is implicit in p87|l . The next order satisfies 



{d-W + 2001 (w) + (^ucj^^iu) = - ^'^{u))R"{u) ~ ?,cp'^{u)R"'{u) - 6R"'(0+)^o(w)] + m{u), (199) 

I 

where ?7o(u) = ■0g^i^'(M) -I- 2g'"^^\u) + gf^^\u) represents pendix and keep only the feeding {0{R'^)) term: 
the feeding from the fourth cumulant. 

{d - 3d + 3C)Mu) + Cui'oiu) (200) 



To determine i/'o (u) we use the expansion of the outer pi 



equation for S^^^ to cubic non-analytic order (see Ap- ~ 7X 



^R"'{0+) \R"'{u)f + R"(u)R"'(0+)R""(0^ 
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Note that the r.h.s starts as \u\ at smah u. The function 
r]Q{u) is determined from the appropriate expansions of 

I 



the 0^^^^) equation as detailed in AppendixIFl It satisfies 



^{R"'{uf ~ R"'{0+ff + 12 R"{u) {R"'{ufR""{u) - R"'{Q+f R""{0+)) + R" {uf R""{uf\ .(201) 



r 



It is straightforward to solve the linear equations 
(|200I201() for ipo.rja, and substitute them into (|199|l which 
in turn is easily solved. For arbitrary C, the resulting beta 
function to this order now involves double convolutions, 
generalizing the single one obtained in f H187|) . Thus, for 



J 



simplicity of presentation we give the result explicitly for 
the particular case C = 0. This choice is of formal in- 
terest for comparison with e-expansion calculations. One 
obtains then 



diR{u) = = (e - 4C + Cudu)R{u) 



1 ^ 



+-R"{uf - R"{Q)R"{u) + 2-i-is{R"{u) - R"[Q))[R"'{uf - R"'{Q^ f) 
+77^ [R"(u)[20(R"'(u)2 R""{u) - R"'{0+f R""(0+)) + 2R" " {u){R"' {uf - R"'{Q+f) 
+6 {R"'{uf - R"'{0+ff + 4 R"{uf {R""{uf + R"'{u) R^'H"))] 



(202) 
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7{R"'{uf - R"'{Q+ff + 2AR"(u) [R"' {uf R"" {u) - R"'{Q+ f i?""(0+)) + 2R"{uf R"" {uf 

I 



where here 75 — 1/(2 — 2e) and 7q = 1/(4 — 3e) (respec- 
tively 73,0 and 74^0 in the notations introduced below). 

This result is also of interest if one performs the 
i?-expansion in another consistent way, closer to the 
e expansion, by expanding simultaneously in Q taking 
C — 0{R). In that case, and to this order, there 
is only one additional 0{R^) term appearing in (|202ll . 
2x^C(20o('") + W(/)o(-u)) with 0o given by (|185|l setting 
C = 0. 

A remarkable property of the above /3-function is that, 
apart from an immaterial constant, its Taylor expansion 
in |m| does not contain a 0(|m|) but starts as u^ + \u\^ + ... 
It thus satisfies the physical requirement of the absence 
of " supercusp" (as discussed in^L^ in the dynamical ver- 
sion this can be traced to a requirement about potential- 
ity of the model). It implies some non trivial identities 
between the zero temperature functions introduced here. 
Examinations of the various terms on the right hand side 
of (|181|) shows that a sufficient condition for the coeffi- 
cient of lul to vanish is: 



7y'(0+) = 



(204) 



<^"(0+)-H2^/.'(0+) =0 



(203) 



One easily checks that these identities are indeed verified 
to lowest order in R, i.e. by the functions </>o, V'o and 
770 ( note that 770 is a linear combination of functions 
which do start as |u| so cancellations occur). A quite non 
trivial property is that there hold actually to all orders. 
From the exact equation (|171|l . it is again a sufficient 
condition for the absence of supercusp. Although we hold 
no general proof we have checked that it holds on some 
higher order terms. It seems to work term by term in 
the FRG equation with cancellations reminiscent of those 
occurring in properties such as dimensional reduction. 

We give here the correction to the second derivative, 
the u^ term in the /3-function, for any C,. One notes indeed 
that one can obtain relations between derivatives valid 
for arbitrary C,. We define ^m,n = — l/a^m,n where the 
Xm.n — d — mO -\-nC, = —(2m — 4 — {m — l)e + (2m — n)C,) 
are the naive bare eigenvalues (rescaling factor) of a 71-th 
derivative of a m-th cumulant at zero. These generalize 
the one defined above in (|104|l . Xm.m = Xm- Computing 
0"(O) with the above method, one finds: 



J 



= (e - 2C)R"(0) -I- R"'{Q+f + 8773,4R"'(0+)'R""(0+) + 4873,4(73,57 + 74,27') R"'(0)' R""'(0) 
+1273,4((273,4 + 473,5)7 + 1374,27')R"'(0+)' R""(0+)' + 0(i?') 



(205) 
(206) 



28 



this can be equivalently obtained by taking two derivative 
at zero of the beta function for arbitrary not shown 
here. 

We now compute the Larkin term to leading order in 
R. It can be expressed by taking the derivative at zero of 
the equation for t/j^^^) which, to lowest order reads (see 
Appendix : 

{d-ie + 304^'^ (u) + Cmi^^^'^'(w) 

= 37'x"4R"(u)R"'(w)3 (207) 
Using Eq. (|196(l yields, for arbitrary ( one finds: 

/4 = ^^-4^^7'x-'R"'(0+)^ + O(i?^) (208) 
that gives: 

K4 = ^^3^ ^^'~l'j^\ ' + 0{R) (209) 

(m2)2 e-29 ' y ' y ' 

given that the exact value in c? = for the case ( = 4/3, 
9 = 2/3 is^ K = —0.16173.. it clearly appears that the 
R expansion should also be performed as an expansion 
in small e and C. For the (relevant) random field case 
C = e/3 one finds: 

'^4 = -^e' + 0(63) (210) 

thus values of order the exact result in d = are already 
reached for e « 1. 

To obtain a reasonable estimate directly in c? = 
(within a fixed dimension expansion) one should clearly 
pursue the expansion to large order, i.e compute the 
fiiR) "function" up to values of R (i.e. its anomalous 
derivatives) where it changes sign. R is presumably not 
small in that regime. In fact we now the exact value of 
^"(0) = 4C2 « 4.2 at the fixed point. 

Let us now reflect on what has been accomplished. We 
have shown explicitly on the third cumulant, and to some 
degree of detail on the fourth one that in d = match- 
ing through the various PBL was possible from the outer 
solution all the way down to the full BL where the true 
derivatives are defined. Since there is little doubt that 
this extends to any cumulant in d = a consequence 
is that one should be able to compute directly, e.g. the 
beta function, from the outer equation, performing a non- 
analytic expansion of R(u) and taking the proper deriva- 
tives at the limit of coinciding points, in a manner free 
of ambiguities. 

This is what we have performed, using mathematica, 
up to four loops (order R^). The three loop result (order 
R^) reproduces exactly (|202|l above and the four loop 
result is detailed in the Appendix [HI 

Let us summarize the procedure here. We use 
schematic notations, parenthesis means a symmetrized 
expression, square bracket means taking the limit of two 
coinciding points and symmetrizing the result. We in- 
dicate schematically the derivatives. Here subscripts do 



not indicate derivatives but the order in the expansion. 
The iteration reads as follows. The lowest (two loop, R^) 
order reads as: 

= isiiR"R"R") (211) 
/3o = [4'^"] (212) 

the three loop order reads: 

^(4) = jQj'{R"R"R"R") (213) 
5f)=75((i?"5f") + [4'^"]) (214) 
Pi = [S['^"] (215) 

We have explicitly checked that the evaluation of the lim- 
its at coinciding arguments [..] never contain sign ambi- 
guities. This is the case also for the direct evaluation 
of 

4 

lim (n5nJ5^''(wi234) = -67Q7'^"'(0+)' 

(216) 

Upon taking the four derivatives, a choice of ordering 
of the arguments can be made. The limit can then be 
taken, and, since the expression is permutation invariant 
in 1, 2, 3, 4, is found to be independent of the choice. Thus 
the function can be defined by continuity at zero. That 
this is also equal to the true derivative at zero, which 
requires some knowledge of the full TEL is of course the 
non trivial property. We have demonstrated in some de- 
tail here how it arises. 
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VI. DISCUSSION AND PROSPECTS 

We have shown that the naive thermal boundary layer 
structure found in previous one loop studies of the FRG 
is really only the tip of the iceberg. In fact the TBL 
is, as we have shown, a truly non-perturbative feature 
of the field theory which extends to all cumulants. For 
the special case d = and TV = 1 we have explored the 
structure in some detail and unveiled the connection of 
the TBL to droplet physics. Also for this case, we have 
also been able to describe in detail the matching between 
the TBL and the zero temperature outer solution. This 
matching solves the ambiguity problem and explains how 
correlation functions calculated in the field theory have 
their expected zero temperature limits. It enables a sys- 
tematic evaluation of the T = beta function in powers 
of the renormalized disorder R{u), given explicitly here 
up to four loops. 

The detailed understanding of the thermal boundary 
layer described here in rf = is clearly a necessary first 
step towards progress in higher dimensions. The ultimate 
goal of such a program is a first principles derivation of 



the e-expansion, of matching, and a test of the droplet 
picture (or its extensions) in d > 0. The ERG described 
in Section IIIII appears as a promising technique to ap- 
proach the problem. We now indicate some first steps in 
this program. 

We employ the multilocal expansion of the full effective 
action functional discussed already in (|79|l : 

oo „ 

in terms of translationally-invariant functions ("j-local 

terms") of j Ui and variables, with Vi^^\u) = V{u) 
in previous notations (and here I — ln(l/m)). The 
ERG equation H81|l can itself be expanded in localitjiSiS 
yielding a set of FRG flow equations for the functions 
Vi and more conveniently for the rescaled functions 

Vi^\ui- ■ ■ Uj;ri...j) — m^'^Vi'\uirn^ ■ ■ ■Ujrn^]Xi,,j) 
with X = mr. By similar considerations as in Section 
we obtain a generalization of the TBLA 



p>l ai,..a2p 



with u = exu/Ti and T = m^T. The fj,2p{xi---j) are the 
multilocal extension of the random force terms and are 
proportional to the j point, fc-th cumulant of the ran- 
dom force (with /i,2 = -R"{0)/{xe)^). The TBLA is 
applied by first separating the ERG equations into cou- 
pled sets by a cumulant expansion (i.e. by number of 
replicas), and then applying the ansatz. Each q*^ cu- 
mulant equation is satisfied by the ansatz to the leading 
0{Ti) order for u ^ 0(1). As in rf = and in the one 
loop Wilson calculation, in this balance the random force 
terms contribute only through rescaling terms while the 
11^^^ functions survive in the remaining terms, the flow 
{—mdm) term being subdominant. Thus one obtains a 
set of coupled flxed point equations for a set of appropri- 
ately defined j-local fc*'' cumulant boundary layer func- 
tions s^^'^'K These equations are regular, thereby solving 
the 1/Ti divergence problem mentionned in the Introduc- 
tion. 

The apparently self-consistency of this generalized 
TBLA is encouraging. Furthermore, it embodies the 
lowest-order droplet scaling as we now describe. In- 
deed, since the ERG gives the exact effective action r[u], 
the above scaling form is a very strong statement with 
implications for all correlations. As in d = 0, these 
are obtained exactly by their tree level calculation from 
T[u\. For instance, for the center of mass fluctuations 



(uo = L j^Ur), 

{{uD - {uoVr (219) 

where the sum is over the number Uy of vertices, 
9?"'''^"* (m)]"'' I indicates schematically a combi- 
nation of derivatives of the local replica TBL function 
v^^''{u), with 2(n-|-ni)'^ total derivatives spread over rt*'' 
factors. The number rii is the number of internal lines, 
and for tree diagrams = riy — 1. Remarkably this pre- 
cisely reproduces the droplet behavior of thermal fluctu- 
ations, (|2J), with an amplitude determined from the local 
TBL functions. Of course, the full TBLA for the effective 
action determines all correlation functions at low temper- 
ature, including those at non-zero wavevector, which will 
have scaling dependent upon q in the manner expected 
from droplets. 

While these successes are encouraging, the further 
steps of matching and connecting to the outer solu- 
tion (hopefully describable in an e-expansion) remain 
a formidable task ahead. In Appendix we sketch 
some of the problems encountered in approaching these 
goals. Another intriguing question is the extension of the 
present methods to > 1 in d = 0. There too, a pre- 
liminary study indicates that the iV = 1 mechanism for 
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the resolution of zero temperature ambiguities is in some 
way modified. A better understanding of the mechanism 
of the (super-cusp) cancellations, which appear to occur 
term by term for iV = 1, is necessary. Of course, de- 
spite the progress made in this paper, we should not be 
surprised that some formidable challenges remain in the 
glass problem. We look forward to successful application 
and extension of the ideas elaborated herein. 



Acknowledgments 

L.B. was supported by NSF grant DMR-9985255, and 
the Packard foundation. Both L.B. and P.L.D. were sup- 
ported by the NSF-CNRS program through NSF grant 
INT-0089835, and CNRS Projet 10674. 



I 

APPENDIX A: GENERAL STRUCTURE OF THE FRG EQUATION 



In this Appendix we study the equation (|24(l . in terms of the rescaled function V{u) defined in l|21|l . namely: 

diViu) = OViu) + J2 CpffTr{{V"Y) (Al) 
p>i 

where O is the rescaling operator, for an arbitrary set of coefficients Cp. We expand the trace and project onto k 
replica terms. We want to evaluate: 

Tr{V"Y (A2) 

we decompose the matrix V" — A + B as: 

Vl = Wa + (1 - (A3) 

We can use the matrix identity: 

p 

Tr{A + BY = Tr{AP) + Yl Yl TriA"" BA""^ BA"'\.BA"'-) (A4) 

P 

= Tr{AP)+J2 J2 {my + l)Tr{BA'''^BA"'\.BA"'^) (A5) 

^"■^ mi,..mi,>0;i'+^^ rrti—p 

as can be easily checked by multiplying with and summing over p. Thus we obtain the convenient representation: 
TriVr = (A6) 

a 

'' = 1 mi,..m„>0;ii+^ mi=p ai,--a„ 

Using the definition of V in terms of the cumulants, we can now obtain the full form of the flow of the rescaled q 
replica term. 

To express it we define the following sets of constraints, CQ{q) as: 



l<m<q~l , {0 < ,^pi^q-l-m (A7) 

i=l 

the set Ci{q) as: 

V 

l<w<9 ,{0 < < 1 - v}i=i.,..v ,l<P = 'y + E"^*-^ (A8) 

and C2{q,v,m, {mi},p) as: 

V 

{PiO,Pili ■■Ptmi}i=l,..v ,Pi]>0, '^iPiO +Pil + ■■ +PinH) = q - P (A9) 
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The FRG equations for the rescaled cumulants (in the Wilson approach) read: 

:^{d-qe + Cu^^U^)S'>{Uc,) (AlO) 
m,{Pi}i=l,..meC'o(<;) 

" 1 

+ Yl Cpg!(m„ + 1)[]. 



v,{m,},peCi{q) {pio,P^l,■■P^m^}i=l....eC2iq,v,{mi},p) PiO ' (1 + Kl ) ! ■ ■ (1 + ) ! 

Sym[Sl^QP^^{Uy,Ul,Ua^o) ^ n '5'2(Mh.o'("l>"ai.i)^ 

( m2 2 \ / '"'^ \ 

JJ^ '^200..0^ ("2, Ua2fc2 ) I •■'^llU.'o ('"'I'-l' ^"i-o) I '5'20o'!'.0 " (^"J """^fcu ) )] 

fc2 = l / \fc„ = l / 

where the first sym is symmetrization with respect to the union of sets ui, ..Um, , ••"q^, with card{uai) ^ pi + 1 
and the second sym is symmetrization with respect to the union of sets ui, ..m^,, itaj^, ..Ua^^^^ , ..Mq^o, ..Mq^,^^ , with 
card{ua^o) = Pm and card{ua-^^ ) = Pik, + 1 for 1 < fcj < m^. 

Let us now derive the flow equations for the second, third and fourth cumulant. The definitions of the rescaled 
quantities are the Wilson ones, Eq. H22|) . only at the end we give the equations for the ERG defnitions pU2|) . 

(All) 



Ci, C2 


are: 










m — 




^ 








V = 


l,mi 


= 0,p = 


l,Pio 


= 1 




V = 


1, mi 


= l,p = 


2,Pio 


= Pii 


= 


V = 


2, mi 


= m,2 = 


0,J5 = 


2,Pio 


= P20 



which yield the following terms, in the same order: 

diR{u) = (e - 4C + (:udu)R{u) + 2cifiR"{u) + 2ci5^i|,(0, 0, u) - AC2R" iO)R" (u) + 2c2R"{uf (A12) 
For g = 3, one gets: 

(A13) 



TO = 


-- 1,P1 


= 1 








TO = 


= 2,pi 


= P2 =0 








V = 


l,mi 


= 0,p = 1,]510 


= 2 






V = 


l,mi 


= 1,P = 2,pio 


= l,Pii =0 






V — 


l,mi 


= 1,P = 2,pio 


= 0,pii = 1 






V — 


1, mi 


= 2,p = 3,pio 


= Pii = P12 = 







V ~ 


2, mi 


= m2 — 0,]3 = 


2,pio +]320 = 


1 




V — 


2, mi 


= 1,TO2 = 0,J5 


= 3,pio = pii 


= J520 


= 


V — 


2, mi 


= 0,TO2 = 1,J3 


= 3,pio = P20 


= J521 





V — 


3, mi 


= m2 = TO3 = 


0,p = 3,pio = 


P20 = 


P30 



These yields, in the order in which they appear (denoting u^- = Ui — Uj): 

diS^^Hui,U2,U3) = (-2 + 2e - 6C + CuAn)S^^Hui,U2,U3) + Sym[3cif;S';^o|)(ui,U2,U3) + 6c2fi^"(ui2)-R"(Ml3) 
+ 3ci^Ji|)o(m1,M1,U2,U3) + 12c2^Ji^(ui,Ul,U2)-R"(ui3) - 6c2^"(0)5;^o^(ui,M2,W3) 
-18c3i?"(0)i?"(ui2)i?"(ui3) - 12C2^{?|,(U1, U2, W3)/?"(mi2) + 6C3R" {U2l) R" (u^s) R" M 

+ l2c3R"{u2l)R"iui2)R"iu23) - 6c3R"{u3i)R"{ui2)R"iu23)] (A14) 

This is rearranged by redefining R"(u) = R"{u) — i?"(0), using the gauge invariance (terms which depend on less 
than q replicas can be erased in a g-rep term) as well as identities from translational invariance (STS) such as 
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sym52oo('^i! W2, Ma) = — 2sym5i^o(ui, ^2, W3). Finally setting Cp = l/(2p) one obtains the one loop Wilson FRG 
equations for q = 2, 3. 

The same is applied to g = 4. One checks that the change R"{u) — > R"{0) + R"{u) produces terms proportional to 
^"(0) which cancel or are only gauge terms. It thus simplifies into: 

dlS^^\ui234) = {d-4e + C«ia„J5W(«i234) + Sym[25^?|,oo(«11234) + 2T,5^^^o(«1234) (A15) 
+3f55^g|,(ui23)R"(wi4) + 47flR"(«i2)R"(ui3)R"(ui4) 

+6R"(ui2) (5^0(^1134) - S[%,im234)) + 6{S[%{un2)SifoM + S[%{m23)Si%{u,2A)) 

+ 127(5^?^(un2)R"(«13)R"(Ml4) + S^?^(«123)R"(W14)(-2R"(«12) + R"(«24)) + ^^^00("123)R"(W14)') 

+67'(R"(wi2)'R"(til3)(2R"(tXl4) + R"{u2i)) - 2R"(tX3l)R"(til4)R"(«12)R"(«23) + ^ R" (^41 ) R" («12) R" (^^23) R" (^34))] 

with 7 = 7' = 1 for Wilson and 7 = 3/4, 7' = 1/2 for ERG in d = 0. 
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APPENDIX B: STRUCTURE OF CUMULANTS 



We discuss the structure for n = 0. The disorder cumulant S^'^^ui, ..Uk) is a sum of (symmetrized) monomials of 
even degree N, u\^u^ ..u^^!' with 1 < Pfe < Pk-i ^ Pi and X^iLiPi ~ ^- Symmetrization is imphcit below so we do 
not always write the sym symbols. The number of possible such terms of same degree N is A{N^ k). These however 
must come in a smaller number of independent linear combinations (of same N) satisfying translational invariance. 
We denote it B{N, k). Shifting Ui — > + A each monomial yields: 



k 

,,pi,,P2 „Pfc _^ ,,pi,,P2 yPfc ^_ ^ Pi 

u 



^ES<"?-"r+o(A2) (Bi) 



=1 



imposing that combinations give zero - up to terms which can be gauge with respect to fc-th cumulant, allows to find 
them. It is sufficient to impose the condition on the 0(A) term, i.e. focus on infinitesimal translations. Also each 
replica must appear at least once (otherwise it is a gauge term). Let us determine the lowest possible value of N for 
a given k and the structure of the lowest terms. We must distinguish k even or odd. 

For k even the lowest possible term is simple and unique and has N = k (i.e. A{k, k) — B{N, k) ~ 1): 

5W(^i, -^k) = 5(^^.1(0, ..0)u^U2....Uk (B2) 



which is by itself symmetric and translationally invariant since the translation (|B1|I produces only k — 1-rep terms 
(pure gauge). One can determine easily the next term N = + 2. Let us do it first on the fourth cumulant. There 
are two choices, which give under translation (terms proportional to A): 

— > 3m^M2U3U4 (B3) 
U1U2U3U4 — > AU1U2U3U4 (B4) 

using symmetrization, so the translationally invariant combination is sym[-|uf U2U3U4 — julu^u^u^]. The result for 
the fc-th cumulant is the obvious generalization: 

S^''\ui, ..Uk) = ^fcS'3?.i(0, ..0) sym{^ulu2U3..Uk ~ ^ululu3..Uk) (B5) 

. There is thus a single term for = fc + 2, i.e. A{k + 2,k) — 2, B{k + 2, fc) = 1. In particular one can parameterize 
the fourth cumulant up to 0{u^) as: 

S''-'*'' (U1234) = (?4UlU2U3U4 — -(76^1^21*3^4(2^^ + 2^2 + 2^3 + 2^4 — U1U2 — U1W3 — U1U4 — U2U3 — U2U4 — U3U4) 

with 94 = S'nii(0, 0, 0, 0) and = -5'22n(0, 0, 0, 0). 

For odd k it is more subtle. The naive choice TV = fc + 1 yields again a single term: 

5('=)(ui, ..Uk) = ^4?^(0, ..0) symulu2....uk (B6) 

but, it is in fact forbidden. Indeed, the translation produces XuiU2....Uk, a variation which cannot be cancelled, since 
there is only one term. One then has to try = fc + 3. To find the possible combinations let us discuss first the three 
replica term fc = 3. There are three monomials at order u^: 

u\u2U3 — > Au\u2U'i (B7) 

M1M2U3 3u^U2U3 + 2u\u2U3 (B8) 

u\u\u^ — !■ Qu\u^uj, (B9) 

since there are two constraints, there is only one possible combination uf M2U3 — 2u\u2U3-\-u\u2U^. Again, this extends 
easily to any fc odd, the lowest A^ is thus A^ = fc + 3 with a single term: 

S^^\ui, ..Uk) = ^^S'lJ^.i(0, ..0) sym(uiM2U3 - 'i.ululu-i + ululul)u4..Uk (BIO) 

It is interesting to note that for the third cumulant complex analysis can be used. Indeed S'*-'^^(ui23) can be expressed 
as a function oi z = ui + ju2 + j^u^ and its complex conjugate z {j is the third root of unity satisfying = 1). 
That automatically enforces translational invariance. Since the permutation group is generated by z ^ jz and 2; — > z 
it is easy to express symmetrization. One drawback of this representation however is that gauge freedom is not 
immediately apparent, so we have not explored further this method. 
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APPENDIX C: CORRELATION FUNCTIONS 



We give here some details of the calculation of correlations from the effective action. We compute all four and six 
point correlations in terms of the disorder cumulants. In each case we also give the STS and ERG identities which 
relate these correlations. 

The calculation was outlined in Section UlTI Correlations are obtained from tree diagrams using the exact two point 
function and the proper vertices r(''), rC^), .. with 4, 6 external legs etc. These are then separated in number of 
replicas. The exact two point function is given in (|13()|l . 



1. four point functions 

The calculation of the four point function requires only one F^^^ vertex. Since S'^'^^ starts as only R and 5*^^^ 
can give a contribution. One simplification in the calculation is that dressing the external lines with a R"{0) gives 
zero, so any propagator from an external point to a vertex can be replaced by the free propagator TSat/'m'^- This is 
because if one contracts first these terms one has, for instance: 

i?"(0)^9„,^^(Me,M3,U,,) =0 (CI) 

/ egh 

from translational invariance. Next, a limited range of q-replica interactions can contribute to a given average. At 
linear order in T one finds that for a given correlation with p distinct replicas, interactions with q < p cannot 
contribute (since one can equivalently use "excluded vertices" where replica are constrained to be distinct). For a 
given correlation which is a monomial of order to, interactions with q > m cannot contribute (using un-excluded 
vertices, there remain free replica sums). This is why an observable such as (uaUbUcUd) with four distinct replicas 
gives zero when contracted to linear order with R and S^^\ and cannot involve fifth and higher cumulants either. 

The result for the connected averages (graphs where all proper vertices are connected) for the five possible monomials 
at order u** is given in the text in Ij 134(1 . In each (except the last) a term proportional to S^^^{0, 0, 0) + ^^^^(O, 0, 0) + 
5^^^(0,0,0) naively arises but vanishes by translational invariance. Here we give the disconnected parts: 

iOd^sc = 3{ul)^ (C2) 

{ulub)disc = S{ul){UaUb) (C3) 

{ulul)d^sc = {ulf + 2{UaUbf (C4) 

{ulubUc)dtsc = {ul){UaUb) + 2{UaUb)'^ (C5) 

{UaUbUcUd)disc = 3{UaUb)^ (C6) 

These can be used to check that the part involving R"{0) cancels in the calculation of the second moment of the 
sample dependent susceptibility x = ("^) ~ i'^)'^- the disconnected part: 



{ulul ~ 2ulubUc + UaUbUcUd) disc {{ul) - {UaUb)f = (C7) 

TO 



Let us now give the STS and ERG identities from H114|) and (|121|l . We find five relations between the two two-point 
functions and the five four-point functions. The ERG relations read: 

2C,T{uaUb) ~ ~Tmdm{uaUb) = 2m^{{u\ub) - (ulubUc)) 

2CT{ul) c -Tmdrniul) = m\{ui) - (uluD) (C8) 

where we recall a, 6, c are all distinct replicas, and one must distinguish several cases in the sum over / in (|121|l . 
Choosing respectively O = m^, O = u'^Ub and O = UaUbUc, one finds three STS relations: 

3l^{ul) = {ut)-{ulub) (C9) 

TO^ 

^{2{uaUb) + {uD) = {ulub) + {uluD ~ 2{ulubUc) (CIO) 

TO^ 

T 

3—{UaUb) = 3{{ulubUc) - {UaUbUcUd)) 
TO^ 

These five equations are not independent, otherwise all four point correlations would be determined from the two 
point ones. The difference of the two ERG equations is equivalent to the difference between the two first STS ones. 



35 



Thus only four are independent and we have five fourth order correlations. Generally at order n one will have iV„ 
71-th order correlations, Nn-2 RG equations and N^-i STS equations. -/V„ can be obtained from the boson partition 
sum nr=i = E„ (with iV4 - 5, A^6 = H, = 22, TViO = 42, ..). 

Note that the relation H12()|) is obtained by forming the linear combination of the three STS relations (with coeffi- 
cients 1, —3, 2 for the three successive lines, respectively). 

If we now compute all above averages using the droplet theory, one finds that the above five equations reduce to 
only two: 

2{xl)p = {x\)d ~ {xIx2)d (Cll) 
2C{xl)p - {x\)d - {xIxDd (C12) 

where (0(a;i, a;2))_D = ^ dxidx20{xi,X2)D(xi,X2)- The three STS equations yield a single equivalent one and the 
two RG equations a single one. Since there are three fourth order unknown we are left with e.g {x\) d undetermined. 

2. six point functions 

The calculation of the six point connected function requires either one T^^'^ vertex or two F'-'*^ vertices. It can 
require a priori up to the sixth cumulant. However, terms of order {S^^^Y ^^d higher contain too many replica sums 
to match the incoming ones. In both vertices S"*^^^ can be dropped as it starts as u^. The F^^^ vertices involve four 
derivatives in zero and thus can only involve R or S'^'^\ Thus to obtain e.g. (u^)c, we need to contract with 

1 „ ■^-^ 1 „(?,•) -^-^ 1 „f4i ■^-^ 1 



e/ e/g e/g/i efghij 

+ E 2(2^(4^)^^/^^^ + E. 2(2r^K24r^) (^^^)^^/(^^'^)S^^°'°^ ^^14) 

where in this equation only derivatives are indicated by superscript and schematic notations are used. We find the 
following result for the eleven distinct connected sixth point functions: 



{^uDc = 15q4T - 1596^^ - 15T^R^^\0) - T'^R^'^^O) + mg (C15) 

{ulub)a = 10^4^ - IbqeT^ - 10T3i?(4)(o) _ T^r((^){0) + me (C16) 

{u^uDc - 7q4T - 3qeT^ + 2sqT^ + hT^R'^'^^Q) + T^R^^Ho) + \2T^R^'^\of + (C17) 

{uluDc = QqiT + 'iqeT^ - isQT^ - 6T^R^^\Q) - T^R^^\Q) - 18T^ R^'^'> (Q)^ + me (C18) 

{utubUc)c = 694^ - 9qeT^ + seT^ + 6T^R^^\Q)^ + me (C19) 

(ululu,), = Aq^T - isgT^ - T3i?W(o) _ ST^ijW (0)2 + (C20) 

{ulubUcUd)c = SqiT ~ 3qeT^ + me (C21) 

{ululuDc = 3g4T + 3qeT^ + sgT^ + ST'^R^*^ (0) + eT^i?^ (o)2 + (C22) 

(ululucUd)c = 2q4T + qeT^ + me (C23) 

(ulubUcUdUe) c = qiT + me (C24) 

{UaUbUcUdUeUf)c = TUe (C25) 

Here for simplicity the mass m has been set to unity, but is easy to restore. We have denoted sg = *S'222(0, 0, 0), 

<Z4 = ^Si (0, 0, 0, 0) , 96 = ^^2n (0, 0,0,0) and me = (0, 0,0,0,0,0). 

We now list the ERG and STS identities to sixth order. Choosing the observable O in H114() successively as 

O = it^,0 = u'^Ub,0 = u^u^,0 = u^UbUc,0 = u^uluc,0 = u'^UbUcUd, O = UaUbUcUdU^ one finds seven STS identities: 

^H<) = {<) - (ulub) (C26) 

-^(Mulub) + {u*)) = {ulub) + {utuD - 2{utubu,} (C27) 

^(3{ulul) + 2{ulub)) - {utuD + {uluD - 2{ululu,) (C28) 
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T 

^{3{ulubUc) + "^iulub)) = {u'^UbUc) + 2{ululuc) - 3{ulubUcUd) (C29) 

^mulubUc) + {uluD) = 2{ululuc) + {vlvlc^) - 2,{ululu,ud) (C30) 
T 

-^{2{UaUbUcUd) + 3(ulubUc)) = {u^UbUcUd) + 3(w^M^"cMd) - 4(u^U6UcUdWe) (C31) 

T 

— 5{UaUbUcUd) = 5{ulubUcUdUe) - 5{UaUbUcUdUeUf) (C32) 

There are also five ERG relations: 

-mdr^Tiul) = m\{ul) - {uluD) (C33) 

-md^T{ulub) = m'iiulub) + - 2{ululu,)) (C34) 

-mdmT{ulubUc) = m^{{u^^UbUc) + 2{ululuc) - 3(u^ufucMd)) (C35) 

-mdmT{ulul) = m\2{utul) - 2{ululul)) (C36) 

-mdmT{UaUbUcUd) = m^iAiulubUcUd) - 4:{ulubUcUdUe)) (C37) 

If one now estimates the above correlations using the droplet theory, these twelve relations (seven STS and five 
ERG), to lowest order in Tm^, reduce again to two ! They read: 

5Tm''{xt)p - {xDd - {xIx2)d (C38) 
Tm\4C - mdm{xt)p = {xI)d - {x^xDo (C39) 

for STS and ERG respectively 

Finally, let us indicate the calculation of D3. One has: 

£»3 = {2a^ulc - a^ulc^ - a%c + a%l - a^uDc (C40) 

= 6^i?(^)(0) + 2^i?(e)(0) + 2se^ + 12^R(^\of 
m^" m^^ m^^ m^* 

with 56 = '^'222(01 0, 0). Upon rescaling this yields the expression in the text. 



Finally, the third cumulant of the susceptibility reads: 



{xs-x;r = {{u')-{u)^r (C41) 

= liiUa - UbfiUc - Ud)Hue " Uff) (C42) 
= {ululul - Sululucd + SulubUcUdUe - UaUbUcUdUeU f) c (C43) 

= S6^+6^i?W(0f (C44) 

It is of order Tin" as here no cancellation occurs. It directly gives a physical interpretation to the sixth cumulant 
TBL function (since we have already related R^^\0) to the variance of Xs)- 

APPENDIX D: DETAILS ON DROPLET CALCULATIONS 

1. some thermal moments 

Let us compute some thermal moments predicted by the two well droplets. We use the rescaled variables u = xm~^. 
Denoting y = xi — X2, one has: 

{{x - {x)D = (p(l - p)™ + (-1)^(1 - (Dl) 
for any power m. One can define the general correlations: 



Dmum,..m,=jl{{x-{X))^^) (D2) 



i=l 



with mi > 7712 > ..niq > 2 and ^ - is even. Defining: 

/OQ 
dxiD{xi,xi + y) 
-oo 

One then obtains: 

D — C T ymi+..+mg p)/ \ 



dzfj iz + i-zr-) 

J, z li (1 + 



(1 + z)- 

The symmetry z ^ 1/z can be used. One finds for example for m even: 



c - 

m 



1 4-'"V^r[l + m] 



2m + 4m2 2r[|+m] 

2 3r[2 + m]r[l + 2m] 



m.m.m 



6m + 27m2 + 27m3 r[3 + 3r 

There are three such quantities for the sixth cumulants: 

.1 



D, = {{x-{x)Y)=T^- J y'Diy) 

1 



D,, = {ix-{x)ry^T- I y^D{y) 



— OO 
OO 



1 

D222 = {{x-{x)W=T-J j'D{y) 



2. Check that the solution of toy model satisfies the ERG-STS equations 

We wiU need: 



27r Ai{iX) 
xgix)=l ^^e-— 



27r ^i(a)2 



— OO 



27r Ai(iA)2 



obtained taking successive derivatives. The STS equation reads, plugging in 



r+oo 

-P'(xi) = -g'{xi)g{-xi) + g'{-xi)g{xi) 2 / + y) + g{xi)g{~ 







It can be satisfied if: 
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+ 00 

dyyd{y)g{xi + y) = g'{xi) + h{xl)g{xi) (D18) 



One has: 



— OO 



+c>o 



27r ^i(iA2) J^.oo 2i7r Ai(zi) (zi - iA2)2 

C?A2^_,;,,^^ 1 ,Ai' 



(-r)'(*A2) (D21) 



27r Ai{i\2) Ai 



3 

I . 1 



_oo 27r Mi(iA2) Ai(iA2) 



It works, with h{xi) = x^. Similarly the FRG equation reads: 

/• + 00 

C(P(xi) + xiP'ixi)) = 2 / dyd{y)[{2xiy + y'')g{-xi)g{xi + y) + {-2xiy + y'')g{xi)g{-xi + y)] (D24) 







(D25) 



= 2xiP'{xi) + 2 / dyd{y)y^[g{-xi)g{xi + y) + g{xi)g{-x^ + y)] 
Jo 

if we use the STS equation. Thus we need to check that: 

r+oo 

CP(xi) - (2 - O^iP'i^i) - 2 / + y) + g{xi)g{-x^ + y)] (D26) 







Now we have again: 



+ 00 r+oo AJ(^\.\ r+oa 1 1-+°° 







, 27r Ai(iA2) 7_joo 2i7r yli(zi) (-zi + ^2)^ 

" ^e-^-^— 1— (4^)"(a2) (D29) 
27r Ai{i\2) Ai' ^ ' ^ ' 

dA2 1 _ . »A2A»'(»A2) 

_ 27r'' ^A»(*A2) A*(*A2)2 +^A*(*A2)4^ ^ ^ 

= 5(2^1) - -^(3(2^1) + xig'{xi)) - ^{g{xi) + xlg{xi)) = ^g{xi) - ^xig'{xi) - ^xlg{xi) (D31) 

which we have also checked numerically. One has thus: 

1'+°° 4 2 2 

2 / dyd{y)y'^[g{xi)g{-xi - y) + g{-xi)g{xi - y)] = -9{xi)g{-Xi) - -Xig' {xi)g{~Xi) + -Xig' {-Xi)g[xi) 







= (:P{xi)-{2-0xiP'{x^) (D32) 

and the ERG equation is also satisfied. 

Let us close by noting that under the rescaling l|162l) the function defined in \V)'i\ transforms as: 

D{y) = p^b{py) (D33) 

with: 

D{y) = 2d{y) J dx,9{-x,)g{x, + y) ^ 2d{y) j l^^r^ (D34) 
One can check that the cumulants computed iniii are recovered with the choice p = 2^^/'^. 
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APPENDIX E: SYMMETRIZED TAYLOR EXPANSION OF CUMULANTS 

The matching calculations in the main text extensively involve Taylor expansions of cumulants. This procedure 
exhibits some subtleties since the cumulants themselves are defined only up to the gauge transformation: 

S'-''\ui,---,Uk) S'-''\ui, - ■ ■ ,Uk) + synikifiui, - ■ ■ ,Uk-i)). (El) 

with arbitrary /. Thus in general since the specific function S^*^-* is not defined without specifying a gauge the results 
of the Taylor expansion may themselves not be gauge invariant. 

In particular when two arguments ui and U2 of S^''\ui, ■ ■ ■ , ut) are brought close together it is desirable to avoid 
odd terms in ui — U2- One may for instance expand around the midpoint (ui + U2)/2, or perform the half sum of two 
expansions around each point respectively using the permutation symmetry of S'^^^ . For instance: 

S^^\ui,U2,U3) = g(3)(- "l +"2 ^ ^"^ ^"3) + ^("1 - U2)^ S'2oo(ui,Ui,U3) - S[%{ui, Ui, U3) + 0{{ui - U2)^) 
S^^Hui,U2, us) - l [^(^'("l, Wl, "3) + ^^^^^2, U2, U3) \ - - U2)^S[%iui,Ui,U3) + 0((ui - ^2)^). (E2) 



2 



2 



Note that the zero-th order term in the first case is not gauge-invariant, while it is in the second. We will thus use the 
second method to perform short distance expansions. Note that taking the derivative 818283 on both sides recovers 
a gauge invariant quantities identical in both cases. 

Let us illustrate the procedure to higher orders. Let us consider for simplicity an analytic and symmetric function 
f{x,y) = f{y,x). The resuh is: 

f[ui3,U23\ = /[U13,W13] - ^W?2/ll["l3,Wl3] + ^^12 (2/l3 ["13, Wls] + 3/22 ["13, ""w]) + OlwD (E3) 

This is performed iteratively. When the lowest order term is even one extracts it from the sum and keeps iterating the 
rest. When it is odd, one first symmetrizes the full expression in 1, 2 then writes W23 — W13 — W12 and finally Taylor 
expand again in U12. What is actually done is to use the freedom: 

/[mi3, W23] ^ /["13, ""23] + 5[ui3, U23] (E4) 

where g is antisymmetric in its two arguments. Then g is determined order by order so that all odd terms are made 
to vanish. One can formalize this procedure as follows. If one assumes the result to be: 

f[x,x-y]+ g[x, x-y]^Yl y^"(l>2n{x) (E5) 

n 

then one has 

fix, y) = f{y, x) = i ^(x - y)2"(02„(a:) + <^2„(2/)) (E6) 

71 

g{x,y) =-g{y,x) = ^^{X - yf"'{(p2n{x) " (p2n{y)) 



One can show that any symmetric and globally even function can indeed be decomposed as ljE6p . 

Using the gauge freedom, one can also make to vanish any particular Taylor coefficient. Note the following property 
of the Taylor expansion upon a gauge transformation. Consider: 

S^^Huus) 5(3) (,.123) + h{u23) (E7) 

Consider the expansion: 

S^^\U123) = J2<2MU13) (E8) 
n 

writing h{u2z) = h{uis ~ ^12) and expanding in U12 one sees that the gauge tranformation becomes: 

M^is) ^ + (-l)"A^'^"Hwi3) (E9) 
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Note that each term in that expansion is not gauge invariant, it is only when the series is resummed to all orders that 
gauge invariance becomes manifest. This helps to clarify the proper definition of the function (j){u) in the PBL21. 
First let us note that the feedback from third to second cumulant involves the combination 2(j)2{u) + (t>'i{u) which is 
indeed invariant upon the above gauge transformation, as physically expected. Although one could have in general a 
non vanishing 0i (u) it can be set to zero by the above gauge transformation. The small argument expansion and the 
definition of (j){u) given in the text thus corresponds to a choice of gauge, and thus with this definition is gauge 
invariant (up to a constant). 

Let us illustrate the symmetrized expansion in the calculation of the feeding terms for the third cumulant functions 
(j>{u) and i^{u). Consider the feeding term of the outer S equation. It can be rewritten, up to gauge: 

^R"(ui2)(R"(7ii3) - R"(7i23))' + ^R"(«12)'(R"(^il3) + R"(U23)) + J(R"(ui3) + R"(u23))' (ElO) 

The first term gives |R"'(0+)R"'(Mi3)^|ui2p + 0{u\^), the second gives: 

'lR"'{Q+fuU2R"{u^:i) - R"'(ui3)wi2) + 7R"'(0+)R""(0+)R"(ui3)|?^i2|' + Oiul^) (Ell) 

The cubic analytic odd term is eliminated when performing the symmetrized expansion explained in the Appendix, 
which we also use to perform the expansion of the last term. The latter is analytic and gives — it^27R"('"i3)R"'('^i3)^ + 
0{u\2)- Putting everything together we find 

7[(R"'(0+)2 - R"'(ui3)')R"(«i3)]«L + 7[^R"'(0+)R"'(iii3)' + R"'(0+)R""(0+)R"(ui3)]|ui2|' + Oiut^) (E12) 
yielding the feeding terms given in the text. 

APPENDIX F: FOURTH CUMULANT MATCHING 

In this Appendix we give details about the matching procedure explained in the text. We first derive general relations 
between asymptotic behaviors of the various functions which parameterize successive partial boundary layers, based on 
consistency requirements that the TEL scaling holds. Then we derive exact fixed point equations for these functions, 
for which, in some cases formal solutions can be given. The general matching relations can then in turn be checked. 

1. general analysis of the matching 

The possible forms of the full BL, PBL211, PBL31, PBL22 and outer solution, and their functions, are defined 
in H193|) . Let us first give some useful properties of these functions. We can assume permutation symmetry in the 
u variables inside each TBL function, e.g. that s*-'^^-' ('U123; M34) is symmetric in {(123 independently of the variable 
W34 since it is true up to corrections that are higher order in Tj. The other symmetries are the following: t/;*^^^) is 
odd, s^^^^(ui2; M34; U13) = s^^^H^34; ^^12; ^13) — s^^^^('iti2; U34; — U13), (/(^^^^^ and s^^^^^ are symmetric in the two last 
arguments (exchange of 3 and 4). Note a quartic polynomial can be incorporated in s^^^^ and thus one need to 
consider only degree two and three polynomial. Degree one polynomial produce gauge terms. Other possible cubic 
combinations vanish, symi23(ui2W23U2i) = 0, symi23(u^2) = 0- quadratic combinations are equivalent up to gauge 
terms. Finally, note that only and tp^'^^^' are gauge invariant. 

One might alternatively choose to express the form in the PBL31 in terms of functions of M14 rather than 1*34, i.e. 

5(^^(^1234) =rfxSymi23(wi2W23)V''''^(wi4)+Tfs(31)(^^23;Wl^^ , "123,"34=0(1) (PBL31). (Fl) 

Clearly since 1x3 and ui are completely equivalent in the PBL31, ^i^'^^^ = -(/j^'^^), s^'^^^ — s^'^^^. However, the expression 
in terms of uu is more awkward when taking the necessary limit of large U13, M23, since this form singles out ui over 
W2, artificially breaking the symmetry mi U2, which is preserved by this limit. 

First consider matching of the outer solution with the PBL211. One can expand the outer solution at small U12. 
We can always choose this expansion to be even in U12 (see Appendix and the zero-th order term vanishes being 
pure gauge. Hence we expect: 

^'''(«1234)^xV2<^^'"n"l3,?^14) + x'kl2|'V^^'"^(«13,«14)+0(wt2) 

"12, Ml3, Wl4, U34 = 0(1). (F2) 
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As already discussed for the third cumulant, for the TBL scahng to hold, the first non-analyticity must arise at order 
|ui2|'^. This should be compared to the large 1112 limit of PBL211. This implies: 

s'^^"^(ui2;Ul3,Ml4) ~«i2->oo |m12|^i/'^^"H"13,'"14), (F3) 

which defines V^^"^- To match, one has then (^(^n) = 0(2ii) ^ ^(211) ^ ^(211)^ 

Next consider matching PBL211 to PBL31. From PBL211, we must take the limit of small M13: 

0(211) (y,3, = 0('"Ho,^/34) - «i3^^'''(^i34) + ul,~g^^^\us4) + ul-iSgn{u,3)h(^^\u3^) + O^g), (F4) 

s(2")(wi2;Ui3,Ui4) = s(2"^(wi2;0,U34) +Wl3/5(^^Hwi2,W34) + Iwisl^^^^^ ("12, U34) + 0(w?3)- (F5) 

We have omitted a naively possible |mi3| term in the expansion of cf)^^^^"^ (uis, uu) which again is forbidden by matching 
and TBL scaling. As discussed in Section M such a term would give rise to a supercusp. 

It is useful to note at this stage that (j}^'^^^\ui3 , un) has some gauge redundancy. In particular, one may tranform 
0^2^^^ (ui3, ^14) (?!'(2^^Hmi3, U14) + /(mi3) + /(mi4) changing S'^"') by gauge terms. The imphcations for the expansion 
coefficients in (|F4|) are straightforward to determine. In particular, the transformation by /(M13) leads to independent 
constant shifts in -0(^1), 5(^1), /i^'^i) and all higher functions. The transformation by /(W14) leads to 

0(211) (0,„34) ^ ^<^'''H0,U,i) + f{u34), (F6) 

^(3i)(^3^) ^ i;<^^'Hu3i)-f'M, (F7) 
g^^'Hu^d - 3<''H«34) + i/"(«34). (F8) 

Non-analytic coefficients such as h^^^^u^^) are invariant under the latter transformation. 

We should compare to the large -U13 limit of PBL31. This requires some care due to the cubic invariant. It is useful 
to rewrite 

symi23(Mi2W23) = ^("13 + "23) + ^"12 " Wl3U?2- (F9) 

When multiplied by ^^^^^u^^), the first two terms are gauge, and the second term is both antisymmetric (hence 
higher order in T) and does not grow at large ^13, hence is subdominant to the dominant final term ^ U13. Hence 

fh symi23 (ui2uL) V^''^ ("34) - -f^x^ul^uisi^^'''^ (W34) + Oift) (FIO) 
in this limit. Next consider the large ■U13 behaviour of s^'^i); 

s(3^'(m123;M34) = U?35(^^n"l2;W34) + W?3Sgn(Ml3)/i(^^'(wi2;"34) + Ul3P^^^Hwi2;W34) + |mi3 Ict'^^) (ui2 ; U34) + 0(1). 

It cannot grow faster because it must match PBL 211 which is no larger than T^. Comparing the two limits, we find: 

^(31) (y) = ^^^'\u) + ^g''\0,u), (FU) 
9^''\ui2,U3,) = uW''\us4)=ul2{9^''\u34)-l4f'\0,^3^)), (F12) 

/l'''^("l2,«34) = U?2/»^''H"34), (F13) 

P^''\ui2,u,^) - U?2/5'''^(«34), (F14) 
a^''\ui2,U3i) ^ ul^a^^'\u3i). (F15) 

To obtain the dependence on ^^^n) in the ffi-st two lines above, one must take some care with the 0{u1^) term in 
(lF4ll : 

w220(211) (0,^34) = ul2(b'^^^'\0,Uu~Ui3) 

= ul^cl^^^''^ (0, wi4) - uisul2^g''^ (0, Uli) + ^ul^^ul^^g''^ (0, ui4) + • • • (F16) 

The first term above is gauge and can be dropped. After dropping this, we can express the result back in terms 
of ui3 and W34 by writing U14 = M34 + M13 and expanding in M13. Because a gauge term has been discarded at the 
intermediate stage, one does not arrive back at the starting expression! Indeed, 

ii220(211) (0, U34) = -Wi3U?2</'0l"^ (0, "34) - ^"?3^?242"^ (0, "34) + ' ' ' (F17) 
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This leads directly to (|F11IF12|) . Such symmetrization procedure is also discussed in Appendix |E| It is useful for 
calculations to note the alternative expansion: 

^(2"n^l3,"l4) =0('"nO,^^14)-lil3V'^''H"l4)+wL(ff^''n"l4) + ^(''^'(«^^ (F18) 

which defines the partial derivatives. Finally note that t/j^-^^) (up to a constant) and .g*-^^^ are gauge invariant. 
Next we match PBL211 to PBL22. This requires the (different from above) small U34 limit of the PBL211: 

0(211) (y,3,U,4) = |u34|V;('')(ui3) + "^43^''H«13) + 0(|7.34n, (F19) 
s(211)(Sl2;"l3,Ui4) = |u34|a(22)(ui2,Ui3) + 0(1^341'). (F20) 

In the large £134 limit of PBL22, one has generally 

s''^^\ui2,U34:;Ui3) = M345(^^^(wi2;Ui3) + |w34|o-(^^^ (^12; U13), (F21) 

up to gauge terms independent of M34, and we used the fact that s^^^> is even in -it34. Matching gives then 

^(22) ^^(22) 

and 

ff(''HSl2,Ui3) = U?2ff^''H^13), (F22) 

Vi(22) = 0. (F23) 
Next we match PBL31 to BL4. The small U34 expansion of PBL31 can be written as: 

V'(3in"34) = ^i34V''''^'(0) + i^.l4Sgn(w34)^(''^"(0+) + .. (F24) 

s(31)(ui23,U34) =^(^n^l23)u34 + p(^n*il23)|"34| +0(^34) (F25) 

where ct^'') is odd and p'-^-' is even. We have used that is odd. On the other hand the large argument limit of PB4, 
taking point 4 away from 1, 2, 3 (i.e. large U34), reads: 

S*-'*n^il234) = ("34)^5^'^'' ("123) + (u34)^Sgn(w34)/l('^^(Ml23) + W340-('^Hm123 ) + 1^34 |p*-''n^*123 ) + CSt 

(F26) 

where the matching implies: 

h^'Hu,23)^xlul,^'^''^"iO+) (F28) 

5'^n^i23) = (F29) 

and 5('*''(mi23) = because ip^^^^^ is odd. Note that by adding gauge terms the quadratic functions h'^^^ and can 
be made symmetric in 123. Finally, the matching of the Larkin term in BL4 and the small argument behaviour of 
the PBL31 yields: 

U = 2V'(3i)'(0) (F30) 

(as we find below, these derivatives exist and are non-ambiguous). 

Finally, we match PBL22 to BL4. The small U13 expansion of PBL22 can be written as: 

s(^^''(ui2,U34,Ul3) = Wi2^'34^'522o''(0'0'0) + "l3/5(ui2 , U34) + | U13 |ct(w12 , ""34) (F31) 

On the other hand we must consider the limit of PB4 when taking both points 34 far from 12, i.e. the large U13 limit 
with U12 and U34 fixed: 

^(''^"1234) = '«13P("12,'«34) + |wi3k(ui2,U34) + CSt (F32) 

One sees however that one must have p = because p{—u, v) = p{u, w), p(u, —v) = p(u, v) thus p{—u, —v) = p{u, v), 
but p must be odd. Thus the matching implies p = and 

0'(W12, W34) = XCT({ti2, U34) (F33) 
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Finally we can combine the various matching stages to obtain expressions for the Larkin term f^. This is an 
important test since one may approach the full BL from either the PBL22 or PBL31. These two approaches give 

U = 2V'(^^)'(0), (F34) 
/4 = 4fo^(0,0,0)^4.g(22)(o)^ 

Consider first the approach from PBL211. We require V'*-'^^-' which is obtained using 0^^^^-' (iti3, U14) = (/)'^^^^^(iti3, U34 + 

i'^^'Kus,) = -<pg''\0,uu)-<t>T\0,^34)- (F36) 

Thus using (|FTT)| 

^(3l)(^,34) ^ (0,^.34). (F37) 

We also require g^^^-* (^13). From the PBL211, taking small U34, one writes 

'/''^"H"13,'«13 + ■"34) = 0'^"' (^13, W13) + 0Oi"Vl3, "13)^34 + ^^02"^ ("13, '«13)W34 H (F38) 

The first term is gauge and can be dropped. The second is usefully rewritten by symmetrizing in 34 (allowed in 
PBL211). 



J.(211)/ \ 1^^,(211)/ \ -(211)/ \\ , 1 ,(211). N 

9^ '[Ul3,Uu) = -(001 ("13,W13) - 001 (W14,U14))M34 + :7(/'o2 ("13, Ml3)W: 



2 

34 



-i0fi"^(ui3,wi3)"i4 + --- (F39) 



We therefore read off 



~9^''\u^,) = -i4fi)(j,i3,ui3). (F40) 

Note that first derivatives at the point (u, u) such as ^lo^^' ("13, "13) is changed upon adding a gauge term /, and can 
thus be set to zero upon choice of gauge. 

Finally (jF37p and ljF40|) can be applied to equations (|F34p and ljF35|l respectively to obtain the Larkin term. Both 
happily agree, and give 

U = -20fi")(O,O), (F4f) 

the derivative being unambiguous thanks to the TBLA. Note that (j)'^^-^^^ has unambiguous first derivatives at the 
points (m, 0) and (u, u). It has also unambiguous second crossed derivatives 11 at these points {u,u). There are 
ambiguities however in the second derivatives 20 at point (u, 0) (because of the h term) and one should be careful 
there. 

Let us note the useful properties: 

4ioo("ii2,"i3) - ~2ul^g^^'\u,3) (F42) 
4foH"i2,0,wi3) - 2u?2g(22)(^^^) 

^)(0,^) + l^^^''\u,u) - 4'o"^(0,") = -(^(^^^u) + 2.g(3i)(«) + ^(22)(^)) (F44) 

which is consistent with the equation from the PBL21 above since: 

_^ 2g(3i)(y) ^ ^(3i)'(y) + 2.g(3i)(u) (F45) 

~(22)(^)^^(22)(^) (F46) 



44 



2. Outer equation 

Having made precise the various limiting behaviors of the PEL functions we now derive the FRG equations that 
they obey. 

We start by writing the fourth cumulant equation ljA15|) in the outer region where all ~ 0(1). Neglecting 
feedback from fifth cumulant it reads: 

5,5(^H«i234) - - (d-4^^ + C^i^5„J5(4)(ul234)+ sym[6R''K2)(-2xV'"Hwi3,Mi4)-^|^^^ (F47) 

+6(-2x'(/.(^/l2)^^00(«134) + S[%{UU3)S[11{U,24)) 

+ 127(-2X^0(U12)R"(?/13)R"(U14) + ^i?|,(ui23)R"(«14)(-2R"(ui2) + R"(U24)) + ^^^0^("123)R"(U14)') 
+67'(R"("12)'R"("13)(2R"(«14) + R"(^i24)) - 2R"(u3l)R"(zil4)R"(^/l2)R"(«23) + ^ R" (^41 ) R" (^^12) R" (^^23) R" ("34))] 

We are now going to study the small ui2 limit of this equation. Before doing so, it is instructive to give the detailed 
result for the various small argument expansions defined in the previous Section of the i?^ feeding term. These will 
feed in the various RG equations for their corresponding functions. We denote these feeding terms by the subscript 
i?^ but one should not forget that they are only feeding terms, not the true function to this order. Symmetrized 
Taylor expansion to order yields: 



, =0^T^(ui3,^^i4) = sym34[3R''K3)R''K4)R('Ho+)'-2R''K3)R''K4)R(')K3)' (F48) 



R' 

-R"(U14)' R^'\uu)' - 4R"(7.i3) R"iu3i) R^^^H^U)' + 2R"(?/14) R"(U34) R(''("13)' - R"(U34)' R(')(ui3) 

-R"(zii3)' R^'Huis) R^'Hum) ~ R"(ui3) R"(wi4) R('H^^i3) R'-'Hum) 



+ 2R"(ui3) R"(U34) R^^Huis) R^'Hww) - ^R"(^i34)' R^'^(U13) R^'H^m)] 

We now expand this in U13 at fixed 1*34. One finds 

0k'"^(O+,"34) - -^R"K4)' R(^H«34)' 

i''iPiu3i) = 2 R"(«34) R('H"34)' + 5 R"(U34)' R('H^^34) R(')(u34) 



^i'i^("34) = -3R"(«34) R('^(^i34)' 

-g'^^PiuM) = -2 R(=^)(0+)' R(3)(u34)' - 6 R"(«34) R^'\0+f R^^Ho+) - 5 R"(«34) R(^n"34)' R(^H"34) 

-^[R"(U34)' R^'Hu^^f + R"(W34)' R('H«34) R^'Hu3i)] 

g'^\v,,) = - 2R(3)(0+)' R(3)(«34)' + \R^^\u,,t 

-6 R"(«34) R(^)(0+)' R(^)(0) + yR"(^34) R^'\u3,f R^'Hu,,)] 

h^^Piuu) = I R'^'\0+f R(3)(^34) - ^ R(3)(0+) R(3)(^34)^ 

We next expand instead the 0^^^^^ feeding term in M34 at fixed M13, discarding the zeroth order gauge piece (and 
symmetrizing this by averaging with the same expansion at fixed U14). This confirms that ^^^'^^ui^) = and gives 

~g<^^!\u,,) = -3 R(3)(0+)' R(^\u,,f + ^R(3'(^13)' + 6 R"(U13) R'^'("13)' R<''(^13) + R"(^13)' R(^'(«i3)'. 

The small argument behaviours are: 



4'"^(0+,^i34) = -^R"'(0+)V- yR"'(0+f R""(0+)|zi|3 + O(zi4) 



i>'iP{u) = 2R"\0+fu + 12R"'{0+fR""{0+)uhgn{u) + Oiu^) 

tP^^I\u) = -3R"'(0+)*w - ^R"'(0+)3R""(0+)u2sgn(u) + 0{u^) 
" 2 
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g<-^P{u) = -2R"'(0+)4 - 15R"'(0+)3R""(0+)|^i| + 0{u^) 
9^r'\^) = ^R"'(0+)' + ^R"'(0+)3r""(0+)|«| + 0{u') 



2 ' ' 2 

I / \ 



U^l\u) ^ -3R"'(0+)3R""(0+)w + 0{uhgn(u)) 
~g'^^\u) = -^R"'(0+)'* + 6R'"(0+)3r""(0+)|u| + 0{u^) 

Similarly, the Taylor expansion to order |ui2p yields the feeding term for V-''-^^^-' (""iSi ■"14)- It reads: 

V'i'^^lwia, W14) = 2sym[R"(^.i4) R"'(ui3)' R^'Ho+) + iR"(u34)' R(=')(0+) R(3)(wi3)2 

+R"(wi3) R(3no+)R('n«i3) R^'H^^m) - ^R"("34) R('no+) R(3Hwi3)R('n^ii4) + ^R"'(0+) R(^)(0+) R" {u,:,)R" {u,^)] 

with the following small argument behaviour: 



^i'"^ (1*13,^^14) = ^i'"^ (0,1^14) + W13(^R"'(0+)' R"'(^/14) - ^R"'(0+) R"'(U14)') 



+2|ui3|R"'(0+)'(R"'(wi4)' + 3R"(ui4)R""(0+)) (F49) 
We can now perform Taylor expansion to 0(1*^2) of tl^c complete outer equation ljF47p . One finds: 

= X^ul^[{d - 4^^ + 2C + C"i35„i3 + C"i4a„, J0(2") (mi3, U14) (F50) 
+6R"'[0+]V3(2"H«i3, W14) + 2sym34[ - {iP^^'{uu) + 2g^^^\u,,)) R" {u^^) + ^'^^^\u,,) R'^^\u^:i) - g'^^^\ui^)R" {u^i) 

+2R(3)(yi4)42n)(ui3,Ul4) + R"(^/14)42"'("13,"14) ~ R^'^ (^13) ^io"^ (""IS, "34) 

+ R"(wi3)</)fi"^(ui3,Ul4) - ^R"("34)0fi"^ ("13,^^14)] (F51) 

+2x^ui2sym34 [x'^(18'0(ui3)V'(^ti4) + 3(?!)'(ui3) (j)' {uu) - (f){ui4) 0"(ui3) - 0(^34) 0"(ui3)) 

-X'/'"(w14)'S'ioi(mi,U3,'U4) - 0" (0) S'^q^, (mi , U3 , M4) +X0'(ui4)5'}o2(ui,'"3,U4) (F52) 

/q\ 2 

2x0 (wi4)S'yi(wi,U3,U4)J (F53) 

+27«?2x'sym34[l2^(wi4)R"(^/i3)R"'(0+) + R"(^/i3) R"(^/i4) (</'"("i3) " 0"(O)) (F54) 

+ iR"(«14)' 0"(ui3) + R"(ui3) R"(U34) <^"(W13) - R"(W14) R"(w34) 4>" (uxz) 

+ ^R"("34)'0"K3) + 2 0'(ui4) R"(ui3) R(3)(«13) + 2 0'(7.i3) R"(«14) R^'H^ls) + 2 0'(«14) R"(?/l4) R('H"13) 
+4 0'(wi3) R"(«34) R^'^WIS) - 2 0'(wi4) R"(W34) R('H"13) + ^(^^m) R^'^("13)' + '/'(^^34) R('Hw13)')] 

+27u?2sym34 [iR"'(0+)24^|,(ui, U3, U4) + R^'^ (uu)' sf^^{ux,u^, u^) - R!' {u^^) R^^) (^14) S%{ux,u^, u^) 

-R"(U14) R(3)(wi3)^i?l(^^l,M3,?/4) + R"(u34) R(3Hwi3)^i?l(?^l,W3,«4) + R"(ui4) R^^^^m) 4ol ("l : «3, "4)] 

+7''"i20fl"''(Mi3,ui4) (F55) 

where (mi3, 1x14) as defined in ljF48n above and we have used extensively all the already derived relations between 
partial derivatives of (\P'^^'> . Note that we have performed an expansion symmetric in 1,2 as explained in the Appendix. 
Thus terms such as Mf2SgnMi2/i(Mi4)R"(Mi3) which do appear, e.g., in expanding (/!'^^^^^(u2i, U24)R"(u23) are in fact of 
higher order ~ |ui2p- 

Next, the resulting equation for (/)(^^^) ("13, M14) can be further expanded for small U13 (one uses '0(0) — 0, 0(0) 0), 
to obtain an equation for ■)/)(^^). In principle we should be able to obtain it by taking the first derivative with respect 
to U13 at fixed U14. It indeed works for the i?^ feeding term, but for the other terms it produces a non analytic piece 
juisl, i.e. the above equation expands as (schematically): 

= u\^(ur^U{u^4) + |Mi3|i?(wi4) + 0(u\^) (F56) 
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Let us give each piece separately. The H equation gives: 

= (d - 461 + 3C + Cu9u)V'^^^^ (u) - 6R"'(0+)7^io"^ (0, ") (F57) 
+ R"'(w)(-g(22)(o) _ 3g(22)(u) _ ^(3i)'(o) + 3^;^^^'^' (u)) + R" {u)ij^^^'>" {u) - 2R"'(0+)/i(3i) (u) 

+x2(-i50"(O)0'(u) - 9(j)"{u)(^'{u)) (F58) 

+7X^^(30' (u)R"'(0+)2 + 60"(u)R"(u)R"'(u) + 90'(u)R"'(u)2) (F59) 

-37'x"^R"(u)R"'(u)^ (F60) 

note that we write the term ^^^^\o,u) in both equations since it has an analytic piece and a non-analytic one 
(subscript NA) as one sees from (IF49|) above. The H equation gives: 



= R"(u)(5(22)'(o+) + 2.g(3i)'(0+) + ^(3i)"(o+)) _ R'" (0+)(6g(^^\u) + b^'-^'^'iu)) 
-6R"'(0+),/;fo")(0,z.)Ar^- 18x^(0" (0) + 2V;'(0+))^(7.) (F61) 
-67X~'r (0) + 2^'(0+)]R"(«)R"'(0+) (F62) 

This should vanish identically, in the same way that a u^jl^isl term does not appear in the third cumulant, which 
would produce 4'{u) ^ \u\ and the infamous supercusp! 

One easily checks that to lowest order in R the various combinations in this equation do actually vanish. First the 
combination 4>"(Q) + 2ip'{0'^) vanishes. Next to lowest order, one has (from above): 

-R"'(0+)(65(3i)(u) + 5^pi3^)'{u)) = 12R"'(0+)3(R"'(u)2 + 3R"(u)R""(0+)2) (F63) 
-6R"'(0+)i^io"^ (0, u)na = -12R"'(0+)3(R"'(u)2 + 3R"(u)R""(0+)2) (F64) 

so that the non-analytic piece of these two terms in the RS'^'^'^ term also cancel. Finally we note that the feeding 
terms to order 0{R'^) displayed above: 

ff(^.^"(0+) + 2g^^P'iO+) + ^P^^P"{0+) = (F65) 

also exactly cancel. 

When the feedback from fifth cumulant is also included, these cancellations should keep occurring order by order. 
We have checked on the equation for tp that 4>"{0) + 2tp'{0~^) also vanishes to next order from the RS term in the 
equation for S^. There seems to be additional exact relations to higher orders, we find: 

(/>"(0) -h2V''(0+) = (F66) 
24p'(0+) + UiP"{0+) + 40"' (0+) = (F67) 

this seems to work in some algebraic way, term by term, reminiscent of dimensional reduction type cancellations, 
according to a mechanism which remains to be understood to all orders. 



3. Equation for PBL 211 

We now proceed as in the study of the third cumulant matching in the text, where we derived two equations obeyed 
by (j){u) (then shown to be consistent). We have derived above the equation for 0*^^^^^ (ui3, M14) as a small argument 
limit of the outer solution. We now derive the equation for (/)(^^^) (mi3, U14) directly in the PBL211, by examining 
(jA15|l in the regime where U12 = 0{1), all other Uij being of order one. We use the various PBL functions defined in 
11931 Grouping all terms of lowest order 0[Tf ) we obtain: 

= x'wLM + 2C - 40 + C(ui3c'«i3 + "i45„,J]0(2ii) (ui3, U14) + x'4o"^ ("12, U13, uu) (F68) 

(F69) 

+ 2x'^i?2Sym34 [ - ^('')'(«14) R"(^^13) + ^^''^ («14) R('^ («13) 

+ 2 R(3)(ui4) (/.^'i")(«i3, «14) + R"(til4) (/.^'"^("IS, "14) - R^'^(U13) <^fo"H-"l3, ""34) 
+ R"(ui3)4l"^(«13,^^14) - ^R"(w34)0f"^(ui3,«14)] 
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+X^sym34[ - xr"{ui2) s^oo^^O- '"13, "w) +xr"{ui2) 4o"H^'i2, "w, wm) 

-R"(u34) 4o0^(0>'"12, -U13) + R"(W14) Siioo(^l5^1,'"2,'Wl3) + R"(wi4) Siioo(U2, t(2, Ul, Wls)] 



+X^sym34R"(ui4)s^o^^(ui2,Mi3) + ^r-"(ui2)S':^oo("i34) 



'20 ^."■12 5 "is; -r 2' ^"■12;<-'200V 

+2{(^2sym34 [x''(3(/>'(mi3) (I)'{ui4) - H^u) ^^"(wia) - ^^(^34) (/-"(uia)) 

(3) ^ 

-X0"('Wl4)<S'^Ol('"l''"3,'W4) + X?^'(wi4)5'i02('"l''"3,W4) - "^^^ ^"^^ '^T' - 2 X (/)' (uu) S^^i{ui,U3,U4)] 

+sym34 [ - 2x'* 4o^^ (0' ""13) 520^' (^'12, uu) + x^ 4o^^ (^12, ^^13) ^20^^ (^12> "14) + s^mi {^r-,^^-, U2) S^2ooiui,U3, u^)] 

+27-&?2XSym34 [R"(«13) R"(«14) <P"{U13) + ^R"(«14)' <P" M + R"(«13) R"(M34) -^"(wis) - R"(w14) R"(W34) cj>" (uis) 

+ ^R"(W34)' Cl>"{ui3) + 2 0'(W14) R"(W13) R^'H^IS) + 2(/.'(ui3) R"(7ii4) R^'Hwis) + 2 0'(«i4) R"(wi4) R^'H«13) 
+4(/.'(ui3) R"(W34) R^'Hwis) - 2 0'(^ii4) R"(U34) R^''(wi3) + (biu^^) R^^^Hu.sf + </'(«34) R^'H^w)')] 

+27Ui2X~^syni34[R(^)(Mi4)^ S'^o^i(wi, M3, M4) - R"(ui4) R^'^H^w) S'^ok^*!, ""3, M4) 

-R"(U14) R(3H"l3)5ni("l,"3,W4) + R"(U34) R(-''H«-13)5S("1,«3,^^4) + R"(ui4) R^'H^w) ^olWl, ^^3, W4)] 

+27sym34 [ - 2x' r"(ui2) R"(ui4) 4o'^(0> ^13) + 2 x' r"(ui2) R"(«i4) s^o'^("i2, "w) 

+ R"(«13) R"(«14) sfiUwi, Wl, «2) + (^12)' 4oo(«l 

+(2 7r"(wi2) + 37V'(ui2)2)R"(Mi3)R"(ui4) +7'x"'^i?2(0fl"^K3>"i4) " 3R"(«i3) R"K4) R^'\0+f) 

we recall 7 = 3/4 and 7' = 1/2 for ERG. 

Taking the large limit one checks that one recovers the previous small 11^2 limit of the outer equation for 5*. 
The non trivial terms which grow as give respectively: 

T2x3sym34r"(«i2)4oo'^(^i2,^^i3,«i4) - 6x^u?2R"'(0+)^('"H«i3, ^14) (F70) 

-Tf x'sym34R"(«34)4oo (0, U12, -U13) ^ -2x*«?2sym34R"(«34)5^^') (^13) (F71) 

2Tf x'sym34R"(tii4)4ioo("ii2, ^^13) ^ -4xV2Sym34R"(ui4)5^^'^(wi3) (F72) 

(F73) 

Using the full BL of R, 4n(wi22) = - ^^"(uis)^ - /'(mis) ~ ^uj^il - {rZ?) and that {r'^f = 1 + 2x</'"(0) 
yields: 

sym34S^'A ("122)5^^^0 (^^134) ^ -x'w?2'A"(0)4oo(^^i34) (F74) 

and similarly: 

2fjVm34[2x^r"(wi2) R"(W14) 4o'^(«12,«13) + R"(W13) R"{ui4)s'^^o{ui,Ui,U2) + ^r"(ui2)' S§1{U1,U3,U4)] 

~ «?2sym34[24x'V(^^i4)R"(«i3)R"'(0+) - 2x'0"(O)R"(ui3)R"(wi4) + R"'(0+)'4oo(^^i,«3,W4)] (F75) 
the non trivial term in the feeding term is: 

37'T'fr"(ui2)'R"(^ii3)R"(Mi4) ^ 37'w?2R"'(0+)'R"(wi3)R"(^^i4) (F76) 
Let us write the resulting equation obtained from PBL211 when the leading order equation is used. One gets: 



= X^-4o"^("l2,Wl3!Wl4) +X^'r"{ui2) {s^2tta\ui2,Ui3,Ui4) - 5200^^(0, Wl3,Ml4)) " 6x^Ui2?'J!!^ V'^^"H'"l3 , "w) 

X'sym34[- R"(u34) (4oo (0,wi2,-ui3) - 2ul29'^^^\ui3))] 
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+2x'syni34[R"(Mi4) (4ioo(^^ii2, "is) + 2ul^g^^'Hui3))] 

+47X^syni34[R"(ui4) [r"(ui2) (s^^^^ (t2i2, U13) - s^^^' (0, U13)) + ^s^o^^(t2i2, U13) - GrZuui^i^is)]] 
+ [2r"(ui2) + 4il,(«ii2) ^^?2+7(r-"(Si2)'-(r^)'^i?2)]4?o(«i,"3,K4) 

+ [27(4iU"ii2) J^^l2) + 2ir"{u,,) +3y(r"(«i2)2 - (r-)2^i?2)]R"(«i3)R"(^^i4) (F77) 

It can be simplified slightly by using the full BL of R. Note that it is a linear equation for s*^^^^\ the homogeneous 
part has the same structure as found in the third cumulant matching in the text, for s^^-'^' above and s^^\ Thus it 
can be similarly solved, formally, in terms of the (more complicated) homogeneous part. 



4. PBL 31 

Considering the ERG equation IjAlSp for 5'^^^ in the regime where ui, U2 and U3 are close together, and inserting 
the various PBL forms from p93|) as appropriate, we find at leading order 0{Tf): 
The equation is of order 0{Tj^). One finds: 

= -xsymi23[M?2"i3](t^ + 3C - 40 + C{uMdu,M'-''^Hu3i) + 2x'symi23[4ooo("i23, M34)] (F78) 

+ ^XWl3R"'('"34)So2o^(-iil2,0,M34) + 3x^r"('Ui2)(sfio|,(-iiii3, U34) - Siioo(ui23, U34)) (F79) 



+ 3xR"(u34)srooi("l23, M34)) + 2 R"("34)sri00(^1123) (F80) 

-5x'^i?2^"("l3)^(''''("34) - X«?2^13(^(''^"("34)R"(«34) + 3^^'''>' MR'" {U34)) (F81) 

+ 2^^'^'("34)U13?-"('"12) + ^R" ("34)5200(^123) + ^X^»'"('«12)S(20)(W13,M34) (F82) 

+ 9x3^22Mi30'(«34)0"(w34) + 3x'Ml3sfiU^*112)'/''(^^34) (F83) 

-^X^4o^^(0, "34)400 ("123) - 3x^4o^^("l2,U34)(slfo("l23) " sfi|,(uii3)) (F84) 

+7(3x'"l3r"(?il2) V(«34) + UiU2?i3(-12R"(li34)R"'(w34)0"(w34) ~ 18(R"'(u34))'0'(^^34)) (F85) 

-3xV" (t2i2)r"(t2i3)4o'Ho, 7/34) +3xV"({ti3)2420) ("12,^34) (F86) 

+3x^(2r"(ui2)r"(iti3) - r"(i2i3)r"(u23))S(2o) ("12' "34) (F87) 

+6r"(ui2)R"(«34)(4io(^iii3) - 4io(«i23))) + 37 sym[r" (^12)/' (^13)] R"(m34) (F88) 

+7'(9sym[(r"(ui2))'r"(ui3)]R"(w34) - 3r"(ui2)r"(ui3)r"(w23)R"(^^34) (F89) 

_^ t~ ~ „(~ ,, R"(^34)(R"^(^34))^ „ -2 , ^"{u3im"'{u3,))\ 

+sym(ui3Mi2r (M12)) 5 3sym[ui2M23] 5 ). (F90) 



Upon taking the limit 'U13 large the whole equation should grow linearly in ui^ and we should obtain an equation 
of the form: 

^*13[Ul2/("34) +g('"l2,U34)] + |W13|[W?2/(W34) +ff(wi2,M34)] =0 (F91) 

The parts / and / must match the small U13 limit of the outer equation discussed previously (and the corresponding 
non analytic piece should vanish). The parts g and g should match the small {ti3 limit of PBL211. All the obtained 
equations should be independent. We will not present these equations in detail here. 

One notes that again the equation for s*'"^^' is linear. It has a slightly different structure from the previous PBL ones, 
there being an additional first derivative term w.r.t. the TBL variables. This should not prevent it being formally 
solvable. The inhomogeneous piece can be simplified using the equations discussed above, the one obtained derived 
above for i/)^^^), and various other BL equations. We will not give further details here. 
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5. PBL 22 

Considering the ERG equation ljA15|l for S^^'^ in the regime where ui is close to U2 and U3 close to U4 and inserting 
the various PBL forms from H193() as appropriate, we find at leading order 0{Tj^): 

= r"{u34) (s[)20^(?ii2,'ii34,Ml3) " So20 (^12 1 0, U13 ) ) + r" {U12) (Saoo^ (?ii2 , ^34 , M13) " SjoQ (0; "34, "is)) (F92) 
+420 (^12, "34,^13) + S2OO ("12,"34,W13)] (F93) 

The result is rather simple as many terms do not contribute to this order. The rescaling term gives zero to order 
because there is no cubic polynomial allowed by symmetry. The feeding TRRR and RRRR term also yields zero 
to this order. For the TRRR term one checks that the r" {ui2)r" {U34) term cannot appear and for the RRRR term 
one checks that although it appears there are powerful cancellations between the various resulting terms. The other 
terms result from Taylor expansions in U12 and M34 not to high enough order to produce anything not antisymmetric 
(one needs at least U12U34 but this is higher order in T). 

We can solve this generally, remembering that the l.h.s. can be gauge, suppressing the dependence on U13 (which 
acts only as a parameter) . The left hand side can be taken as an arbitrary function of U12 plus the same function of 
M34 (since the right hand side is symmetric in the two variables). Doing this, one finds the general solution 

(22), ~ ~ ^ (22). ~ „ ^ 5(^12, U34) , a(tti2) + a(M34) , , 

S02 (m12,U34,W13 - ^2 Ul2,0,Ui3 = , + ^ ^ , (F94) 

l + r"(u34) l + r"(u34) 

where a{x) is arbitrary and g(x,y) — ~g{y,x) is antisymmetric but otherwise arbitrary. The most trivial but likely 
physical solution corresponds to a{x) — g{x,y) — 0: 



,(22),,, 



({il2, U34, U13) = '«12l'345'^^^("l3)- (F95) 



This matches properly provided that iT({ti2, U13) = 0, i.e. that s^^^ is sufficiently smooth. We have not checked 
whether the (here neglected) fifth cumulant will produce some additional non trivial feeding term to s*-^^^. This does 
not however modify the conclusion that s'^^) can be obtained formally in terms of higher PLB cumulant functions. 

APPENDIX G: BETA FUNCTION TO HIGHER ORDERS IN R 

The iteration to four loop order reads: 

Sf^ = -ip-i^{R"R"R"R"R") (Gl) 

S'f^ ^ -fQ{{R"st^") + 7(i?"i?"5f' ") + [S-^''"]) (G2) 

5f =75((i?"5f^") + [5l'^"]) (G3) 

h = [4'^"] (G4) 

The only term needed not yet determined is the R^ feeding term into the fifth cumulant. The general expression 
for the feeding term of the p-th cumulant 



c'^p\{TrW^- R'U-R'L,) (G5) 

aa\ ..ap 

Wab^SahJ2^ac-Kb (G6) 



where Rab = R{uab), the last term being the p + 1-th replica term which appear in the trace and must thus be 
subtracted. Here Cp = Cp = l/(2p) for Wilson and Cp = \j7P which comes from rescaling for the ERG. As can also be 
extracted from Appendix A, the feeding term into the fifth cumulant equation (not written) reads: 

1275sym[5R"(ui2)R"(ui3)R"(wi4)(R"(wi2) - R"(w23))(R"(^ii5) + R"("25)) (G7) 

+ (5R"(W14) - R"(W45))R"(«15)R"(U12)R"(^/23)R"(«34)] 

with 75 = 1 for Wilson and 75 = 5/16 for the ERG. 
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At each step of the calculation, upon applying the [..]i2 we check carefully that no ambiguity arises. In general, upon 
inserting the non analytic expansion for R"{ui — U2) in powers of |mi — M2I, a very large number of terms proportional 
to sgn(ui — U2) arise. They however amazingly cancel when the limit U2 — > Wi is taken. In fact this serves as a useful 
and non trivial check at each stage of the calculation. 

The final result for the four loop beta function diR = P{R) is: 

04 loop{R) = (G8) 
c^oR"'{ufR""{u) + cnR'"(0+)'R""(n) + (cn - c,o)R"'{0+fR"'{ufR""{0+) - (cio + cn)R"'{0+fR"'{ufR""{u) 
C2oR"(u)R"'(m)2R""(m)2 + c2iR"(w)R"'(0+)2r""(0+)2 + c22R"{u)R"'{0+fR""{0+)R""{u) 

-(C20 + C21 + C22)R"(w)R"'(0+)2R""(u)2 

+C3oR"(^^)R"'(«)'R""'(w) +C3iR"(«)R"'(0+)2R"'(m)R""'(w) - (C30 +C3i)R"(w)R"'(0 + )3R""'(0+) 

+C4oR"(w)2r""(u)3 + c5oR"(w)'R"'(w)R""(u)R""'(w) + ceoR"{.ufR"'{u)^R"""{u) + ceiR"{u)^R"'{0+fR"""{u) 
+croR"iufR""'{uf + csoR"iufR""iu)R"""iu) + c<,oR"iufR"'{u)R"""'{u) 

The various terms have been grouped together with their associated anomalous ("counter") terms. The coefficient aj 
are found to satisfy the constraint that the supcrcusp vanishes for each group i, i.e. the associated anomalous terms 
successfully substract the |u| piece of the main term (indexed j = 0). Terms with even powers of R"{u) do not need 
anomalous term, nevertheless one notes the existence of the cei term. The table of coefficients Cij is found to be: 



ClO 




372 75 7p 7q 7s + 354 7' 7q^ 75 + 248 7^ 7q 75^ + 156 7' jq 75^ + 138 7 75^ 


(G9) 


Cll 




-12 75 7p 7q 7s - 18 7' 7Q^ 7s - 8 7^ 7q 7s^ - 6 7' 7q 7^^ - 12 775^ 


(GIO) 


C20 




67275 7P7Q7S + 6647'7q^7s + 4887^ 7q7s^ + 3647' 7q7s^ + 346775^ 


(Gil) 


C21 




-624 75 7p 7q 7s - 588 7' jq'^'Js - 432 7^ 7(5 7s^ - 300 7' 7q 7s^ - 264 7 7s^ 


(G12) 


C22 




-24 7' 7Q^ 7s - 24 7' 7q 75^ - 32 7 75-'^ 


(G13) 


C30 




192 75 7P7Q7S + 220 7' 7q^ 7s + 192 7^ 7q 75^ + 148 7'7q7s^ + 166 7 7s^ 


(G14) 


C31 




-28 7' 7Q^ 7s - 32 7^ 7Q 75^ - 28 7' 7q 75^ - 46 7 75^ 


(G15) 


C40 




48 75 7p 7q 7s + 64 7' 7q^ 7s + 48 7^ 7q 7s^ + 52 7' 7q 75^ + 52 7 75^ 


(G16) 


C50 




96 75 7p 7q 7s + 184 7' 7q^ 7s + 144 7^ 7q 7s^ + 172 7' jq 7s^ + 224 7 75^ 


(G17) 


C60 




24 7' 7qSs + 32 7^ 70 7s^ + 24 7' 7q 7s^ + 50 7 7s^ 


(G18) 


C61 




-277s 


(G19) 


C70 




47'7q^7s + 47'7q7s^ + 12775^ 


(G20) 


C80 




47'7(3^7S +47'7Q7S^ + 1677s^ 


(G21) 


C90 




477S^ 


(G22) 



We have also computed the correction to the fourth cumulant, yielding: 

h = ^_^g^^^ X-'R"^(0+)'(7^ + 287'7qR""(0+) + 16(2757P + 7%s)R""(0+)) (G23) 
up to 0{R^) terms (strictly for C = 0). 

APPENDIX H: ABOVE FOUR DIMENSIONS 

In this section we consider the behavior for d > 4, and in particular the fate of the TBL in this case. Four 
dimensions plays to a certain extent the role of an upper critical dimension. In particular, for d > 4 and a weak and 
smooth pinning potential, the Larkin approximation (of Taylor expanding the pinning potential to obtain a random 
force) is self-consistent in the sense that the deformations {{u{r) — m(0))^) remain bounded and small. The self- 
consistency of this argument suggests that under these conditions there is an (almost) unique local minimum energy 
conformation of the elastic medium, and consequently no glassy behavior. (We include the proviso "almost" here to 
indicate that the self-consistent Larkin argument, which considers only the variance of the manifold displacement, 
may incorrectly miss rare samples/regions with anomalously strong pinning, arising from the tails of the disorder 
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distribution.) Nevertheless, clearly for strong enough pinning - or a weak but non-smooth pinning potential - there 
must be multiple metastable states for the manifold. 

To analyze the nature of these low-energy states and their consequences for the boundary layer, we consider both 
the FRG and a complementary 1/d expansion approach, both of which support the same qualitative behavior. For 
a weak, smooth, pinning potential, both approaches indicate that the effective action remains analytic at T = 0, 
and consequently that the TBL is absent. Since this implies the variance of the curvature of the effective potential, 
(F"(m))^ is bounded (and indeed small - see below), this indicates the absence of metastability. For stronger disorder, 
the effective potential does become non-analytic at T = 0, and there is thus a boundary layer for T > 0. Unlike for 
d < 4, however, the width of this boundary layer is 0{T) rather than 0{T/L^). Our interpretation of this result is 
that the metastable states for d > 4 occur at small scales, and hence between states with bounded energy differences 
from the ground state. Physically, even for strong disorder the manifold remains asymptotically flat on large scales 
(with a bounded displacement larger than the transverse correlation length of the pinning potential) , because it is too 
costly in elastic energy to distort the manifold over large distances. The metastable states thus correspond to small 
local distortions of the manifold atop this generally flat background. 



1. FRG for d > 4 



We consider the Wilson FRG. First, note that the distinction between strong and weak pinning can be observed 
from very simple considerations at T = 0. In particular, from l|25f) . and ignoring higher cumulants (e.g. 5^^^) for the 
moment, one finds by differentiating: 

azi?(4)(0) = ei?(4)(0) + 3 [i?(*)(0)j' + f',i?(*^)(0). (HI) 

For d > 4, e < 0, and for i?W(o) < |e|/3 (at low T) it flows to zero, as do all higher derivatives and cumulants. 
Taking ( — e/2(< 0), only the quadratic Gaussian Larkin force term survives under renormalization. This is the weak 
pinning state. From (|H1|) . however, for stronger disorder (and low T), one moves outside the domain of stability, and 
the disorder correlator becomes non-analytic (at T = 0) at some finite scale ^ = Zi as for d < 4. 

In this case for Z > /i we must return to general functional analysis. The detailed behavior of the FRG for scales 
of order li is quite complex, both inside and outside the boundary layer. However, considering the behavior outside 
the boundary layer (i.e. at T/ — 0), it is clear from H25|l that the full force-force correlator A;(m) does not have a 
stable fixed point, and that ultimately the flow of A; will be toward zero. Higher cumulants also flow toward zero, 
and moreover do so more rapidly with scale than does the second cumulant. Hence there will be another scale I2 > h, 
beyond which the entire function R{u) will become small, so that the 0{R^ ~ A^) terms can be neglected. Once this 
occurs, there is no longer any feedback to continue to narrow the boundary layer, which will have some width set by 
its narrowness at I = I2, presumably of 0(Ti^). In particular, in the linearized FRG equation for A;, the TjA" term 
can then be neglected (presuming as usual T/ <^ 1), and one has 

A,(u) - e~l'l'A(u) torl>l2. (H2) 

Here to avoid confusion we have not made any rescaling, i.e. used C — 0. In the same scheme it is appropriate to take 
6 ^ d — 2. This behavior obtains hence uniformly for all u, indicating that indeed there is no further sharpening of the 
boundary layer. The physical interpretation is that the important metastable states occur at finite scales L ^ A~^e'^, 
with consequent finite energy separations. 



2. Large d Limit 



Next we consider the large dimension limit d ^ 00. The FRG analysis sketched above suggests that the qualitative 
physics is independent of dimension for d > 4, so we may expect the large d approach to recover similar behavior. 
The results to leading order are identical to those for a fully-connected "mean field" model, so we will present the two 
in tandem. For the large dimension analysis, it is convenient to consider the model on a d-dimensional hypercubic 
lattice. 

pH[{u,}] = il^(u,-u,)2+^y,(u,), (H3) 
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where K = K/T and Vi{ui) = Vi{ui)/T. We consider in particular the free energy (effective action) ro[Mo] for the 
uniform (zero) mode of the displacement field Uq, in a given disorder realization, defined according to 

ro[Mo,A] = -lny"[d?/,]e-'^^[^"'>l+S.^("'~"°\ (H4) 
d),To[uo,X] = 0. (H5) 



Here A is a Lagrange multiplier used to define ro[uo] — Tq[uq, X(uq)] in the usual way via the condition in IHSp . 
One may also readily calculate a variety of other similar quantities at large d, e.g. the full effective action for all 
coordinates in a particular disorder realization, or the effective action for the replicated model. The present choice 
seems simplest and the most physically transparent. 

The large-d limit of ro[uo] can be calculated using the efficient method of Georges and Yedidia,^^^ which gives the 
result, to leading (zeroth) order in l/d, 

ToK] = ^(ff.(A) + A«o) - f (H6) 

where A(mo) is determined by the condition ^^igliX) + uq) = and 

(?,(x) = - In / dwe-^"'-^'(")+=^". (H7) 



is proportional to the free energy of the "toy model" with potential Vi{u). To leading order in 1/d the last term in 
(|Ij^6|) can be simplified to ^ 9iW9jW ~ -^""o' where N is the number of sites. With this replacement, the result 
is identical to the exact effective action for the fully connected mean-field model, 

ij i 

The effective action for this model is easily and exactly obtained by mean-field theory. One finds 

TfcAM = ro[uo] - E + ^^^*o - Nkul, (H9) 



with the same condition l|H5p determining A as above. Using the fully-connected form in ljH9|l . one can then generally 
show that the renormalized random force obeys 

^^^^=A-2.o, (HIO) 
NK 

where X — X/K. 

We now use this result to evaluate the two-point function of the "renormalized pinning potential" defined by 
ro[uo] (note that because we consider the uniform mode, there is no mean elastic contribution to the energy for a 
displacement Ug). For concreteness, we focus on the periodic CDW case, and take Vi{ui) — V{ui — Pi)/T, with /3i 
uniformly distributed between and 27r. 



Renormalized disorder correlator for weak pinning 



Consider first the case of a weak and smooth potential, (e.g. V{u) = a cosm for a <C K) with | V^(it)|, | V^'(m)| ■ ■ ■ ^ K. 
In this case, we expect from the FRG analysis, that the renormalized disorder will remain analytic at T = 0, and 
hence that there will be no TBL. Thus we concentrate on zero temperature, for which 



9.{X) 



Ku^ + V^{u) - Xu 



This can be evaluated perturbatively in Vi{u) to give 

9^{X) 



4K^ 



(Hll) 



(H12) 
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Then ljH10|) gives perturbatively 

~X^2uo + ^J2 i^/K) - ^ E i^VauoWiuo) + ^ E ^VauoW^'iuo) + OiVyk'). (H13) 



It is straightforward to obtain the force-force correlator by directly averaging using Eqs. (|H10IH13|) . One finds 



{T'aiuo) - nK))2 = NK^ [Ai(uo - u'o) + - u'^) + 0{V^/K^)\ , (HI4) 

with 

= >-jg[VCTI'[>-"(,a + ^.) + v"(/i-~.)l 

Jo 27r X-^ 

Provided the potential V{uo) is smooth (the existence of all derivatives is a sufficient condition), the force-force 
correlator in l|HI4l) is analytic at zero temperature to all orders in V/K. Since it is analytic at zero temperature, there 
is a smooth limit as T — > 0, without any boundary layer. This can also be explicitly verified by a direct calculation 
of (?i(A) perturbatively in V from IjHTp . This describes the weak-disorder behavior in large d. 



b. Renormalized disorder correlator for strong pinning 

For stronger disorder, or for a non-smooth pinning potential, the above analysis does not obtain. We now consider 
the behavior explicitly for such a case, taking for concreteness the Villain potential: 



^g-f(u-/3,-2.n)^^ (HIT) 



with P uniformly distributed between < (3 < lix. In the low temperature limit with a = a/ K fixed as T — > 0, the 
single site "free energy" is then given approximately by 

n 

9n{P, A) - ^(/3 + 2^n)2 - ^ (A + a(/3 + 27rn))' . (HI9) 
2 2(a + 2) 

Consider first the zero temperature limit, for which this simplifies to 

5,(A) =min5„(/3„A). (H20) 

n 

The minimum occurs for a given n when A„_i_„ < A < A„.„+i, where 

A„,n+i = 2/3 + {2n + I)27r (H2I) 
which is the point where <?„ = gn+i- The derivative of g has a jump at these values with: 

g'iiX, T = 0) = y — ^(A + a(A + 2TTn))9{2f3, + [2n - I)27r < A < 2/3^ + (2n + I)27r). (H22) 
^ — ' a + 2 

n 

Hence g[a,tT = Q consists of alternating linear segments with slope — I/(a + 2) , and jumps of magnitude g'(A„^„+i + 
0^) — g'{\n,n+i ~ 0^) 27ra/(a + 2). At low but non-zero temperature, the jumps in g[{X) are rounded over a region 
of width 5\ - R-^ = T/K. 

To obtain the second cumulant of the effective potential we use that T'q{uo) — A which gives: 



(r;,(uo) - T',{u'^)Y = N'K\\{uo) - XK) - 2{uo - u',))^, (H23) 
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FIG. 2: The function U{X) defined in tfie text for the A''-site d = oo Villain model for the simplest limit a — oo. At zero 
temperature, it consists of N "shocks" (step discontinuities) in each 47r interval of A, each of magnitude 2n/N. At T > 0, these 
steps are smoothed by an amount 5\ ~ — T/K. 



where as usual A is determined from 

uo = L/(A) = -l^g:(A.). (H24) 

i 

For simplicity (though it is not essential), we now study the limit a = oo where g[{X) is flat and jumps exactly by 2t: 
at each shock. Note that the infinite a limit amounts to making the model discrete, since each coordinate Ui is then 
constrained to sit exactly at one of the minima of the Villain potential, m — Pi + 2'Kni. This model is still non-trivial 
since K is still necessary to hold all particles near to one another (i.e. to keep Ui — ni> small). The function U{X) for 
this case is shown schematically in Fig. |5] Graphically the solution to ljH24|l is obtained by intersecting a horizontal 
line with the function Uo(A) for a given disorder realization. We consider first the behavior at T = and show that 
it indeed exhibits the expected absolute value cusp in the force-force correlation function. From the figure it is clear 
that (at T — 0) the solution to (|H24I) is always such that A is equal one of the shock values A = 2f3i + {2ni ~\- l)27r), 
with Ui an integer. Let us, for a fixed iV, order the shocks by increasing values of A, i.e. A^ < A^' for m < m' . 

To obtain the two point statistics of the A(uo) random function we intersect with two lines U{X) — uq and U{X') = Mq. 
For given values of uq and u'q and a given configuration of disorder, the solution is U (Am) — uq and U{Xm') ~ Uq with 
the above-mentioned ordering. Clearly from the figure one has, since U{0) = jj ^ - f3i, 

m-m' = n= _ , (H25) 

where [x\ indicates the greatest integer less than or equal to x. This involves the "center of mass" of the Pi variables. 
On the other hand we know that the A™ are of the form 2/3^ -I- {2ni + l)27r) for some Pi^rii but the differences 
Am — Am' = 2Pi — 2Pii depend only on the differences Pi — Pii . Thus we can use that the differences are statistically 
uncorrelated with the center of mass. This is true because the distribution of the Pi is uniform. Thus we may consider 
for fixed n = m — m' the distribution over the differences. 

This problem is formulated as follows. One chooses N independent variables 6i = 2Pi on the interval [0,47r] with 
uniform distribution. Since the quantities of interest depend only upon the differences, we may arbitrarily fix the first 
one at zero, and order them as f?o = < < 02 < ••• < ^N-i < 0n ~ 47r. The distribution oi ip = On ~ Oq is the 
object of interest. Switching to difference variables 5k — Ok — Ok-i, k — 1 . . . N , this probability reads 

„ N N n 

PnW =ZnJ l[dS^6{Y,S^ - 4^)<5(^<5, - V), (H26) 

i—l i—1 i=X 

where Zn is a normalization constant. The first delta function here expresses the constraint that the sum of the 
difference variables gives 0^ — 9q = An. Standard techniques (using a Fourier representation of the delta functions) 
give the normalized result 

«.(*) - (4^)1- r(„)r'("'-^) - 

Thus, finally, the probability distribution of Am ^ A,,,/ is obtained by the sum 

P{X{uo) ~ X{u',) = ^) =^P„(^)Prob{[ ^""~J^'^' ] - ^ n}. (H28) 



55 



Clearly the fluctuations of n are 0(1), while n itself is 0{N). Hence, in the large system limit ^ 1 the center of 
mass disorder is negligible and we can take 

X ^ n/N = {uq- u'o)/ (271). (H29) 

Then, using the Sterling formula, one has 

^ exp(27r7V(uo ln(A) + (2^ - uo) ln(-i^^^)) (H30) 
Iuq 2{2tt — uq) 

where the denominators have been chosen so that it is approximately normalized, and we have taken u'q — without 
loss of generality since P is translationally invariant. Due to N in the exponential, it is strongly peaked. We write 
'0 = 2uq + / to bring out the physical /, representing the fluctuation of the force difference per site (in units of K). 
Then 



-Pri(V') ~ exp 



2.iV(.oln(l + ^) + (2.-.o)ln(l-^^^) 



exp 



(H31) 
(H32) 



D{uo) = ■\uo{2tt - Uo), (H33) 

TT 



-2 Uo 



3 Uq{2t: — Uo) 
The force correlator is finally obtained as: 



(T'oiuo) - r^,K))2 = N'K'i\{uo) - \{u'o) - 2(uo - u'o)f = K'Nuo{27r - uo)/tt\ (H35) 

Indeed, as expected, this form has the linear cusp at small uq. Higher cumulants are readily extracted by further 
expansion around the saddle point / = 0, e.g. 



(r[,(uo) - r[,K))3 = -24iVi^^^(uo)I5(uo)^ (H36) 

This also exhibits the non- analytic behavior expected at T = 0. 

Turning finally to T > 0, we note that shock contributing to U{X) in Fig. |21is smoothed on the scale 5X ~ T/K. 
For large N ^ T/K each shock will overlap many others over this range. Therefore the T > solution is expected to 
be fully analytic in uq, and rounded over a scale ^ T/K. Thus one finds agreement with the naive FRG treatment 
for d > 4 of Sec. lH II the thermal boundary layer width does not diverge with scale (in this case N since we consider 
the center of mass) but remains finite. 

APPENDIX I: ISSUES IN EXTENSIONS TO > 1 OR d > 

1. more than one component 

In this Appendix we attempt the calculation of the two loop /3- function for d = and A^ > 1, by the direct method 
of expansion in powers of R of Section IV CI and show that it fails. 

In the N component model the zero temperature exact RG equations truncated to order B? read: 

-mdraR{u) = - K + Cu'du^)R{u) + ]^R'l^{uf + ^^.^^^,S{ul,U2Mu^=u,=u (II) 

-ma„S'(ui,U2,M3) = (-2 + 2e-6C + C<5„>)5(ui,U2,U3) (12) 
+7(3sym[R;;.(ui2)R;'fc(wi2)Rl',(^ii3)] - K,{u^2)R%{u23)R'L{u:ii)) (13) 

where summation over repeated indices is implicit, arguments are now A^ components vectors u\ i = 1, ..A^. Consid- 
ering the model with 0{N) symmetry, we define the functions: 

R{u) ^ Nh{u^/N) (14) 
h{x) = b{x^) (15) 



56 



for any N. To compute the beta function by the direct method requires the calculation of the feedback of S into R, 
denoted as: 

f32ioop = isi[Tr{R"R"R")"] (16) 

where we use the schematic notations defined in Section FV CI and the Tr denote the trace to be performed in 0{N) 
space (and Q is set to zero for now). 

Performing the necessary contractions, and assuming that h{x) has a well defined Taylor expansion near zero with 
h' [x) = we find the fixed point equation for h: 



1 T,u ^2 N-lh'{xf 



+4,(^/."(.)^(.'"(.)^-/.'"(o)^) + (^-^)^^'(^)' 3/.'(.)^/."(.) , i-h"i.r 



-''^^h"'m{N + i)^ + h"{x)]) 

where h'{x) = h'{x) — h"{0)x. The first line contains the rescaling term and the one loop 0{h^) terms, while the 
second and third line is the result of the calculation to 0{h^) of the above-mentioned feedback. Extension for finite 
C using convolutions as in (|187|l is straightforward. The non-anomalous terms, present in an analytic theory, have 
the same structure as obtained in previous similar field theory calculations. The anomalous terms all contain h"'{0) 
which plays the role of i?"'(0) for = 1. In fact, for = 1, h{x) identifies to R{x) and one recovers the beta function 
obtained in 1)187(1 . One nice feature of the calculation leading to the above is that when taking vectors ui and U2 
close one can check explicitly that terms containing (mi2 ■ iiis)^ vanish, thus the result does not depend on how the 
points are brought together, i.e. it is unambiguous. 

However, the above is not an acceptable /3-function. Indeed one finds that the behaviour at small argument of the 
r.h.s. of (jni) is: 

-l^^h'"{0fx + O{x') (IS) 

This singularity (supercusp) occurring in the feeding term at two loop is inconsistent with the working hypothesis 
(that the behavior be similar to the one loop result which starts as m^), necessary for this direct R expansion method 
to work. Thus, for A'^ > 1 one of the assumptions leading to this expansion must break down. 

One also notes that a two loop term survives for N = oo, which is surprising in view of the results found in (). It 

reads: ~^h"'{0)'^^^-^ — —jh"'{0)'^b'{x^). It is easy to trace it to the expansion of the leading order term: 

S = i^slNb'{^)b'i^)ib'i^) + b'C^)) + 0(1) (19) 

expanded to order u^j, keeping only the leading order in N . 

The fact that the most naive extension of the non-analytic expansion in power of R fails here signals that the 
matching problem must be reexamined for > 1. Work is in progress in that direction. 



2. higher dimensions 



It is instructive to consider how the naive one loop zero temperature FRG is encapsulated in the ERG. Since it 
requires keeping only terms up to second order in R{u), it is straightforward to write the corresponding multilocal 
terms used in this formalism. We require up to bilocal order for the second cumulant, which is a term in the effective 
action of the form 

^I'^-^f E^«'<'=^12). (110) 

Naively truncating to this order, the rescaled local and bilocal disorder correlators obey the equations 

diR{u) ^{d-29 + Cudu)R{u) (111) 
+fdg{x = 0)R"iu) + 2f f dg{x)Rii{u,u,x) + \{ f dg{x)g{x)){R"{uf ~ 2R"{0)R"iu)) 

Jx J X 
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and 

diR{ui,U2,x) {2d - 29 + (uidu, + Xidx,)R{ui,U2,x) + F{ui,U2,x) - 6''-{x) / F{ui,U2,x) (112) 

J X 

F{ui,U2,x) = fdg{x = 0){R2o{ui,U2,x) + ^02(^1, "2, a;)) + 2fdg{x)Rii{ui,U2,x) 

+ \dg{x)g{x)R"{u,)R"{u2). (113) 

We use the same rescaling factor 1/4 as in the text for second cumulant (and an extra mf" for the bilocal term). 
The projector in II12I extracts the purely local part. Here the Fourier transforms of g{x),dg{x) are g{p) ~ pq^fj 

9g{p) = (p2^i)2 ■ Note that if T is taken to zero in lllll) . it becomes exactly the naive one loop FRG equation of the 
Wilson approach for R{u), with dg{x)g{x) — l/(327r^) in d = 4 (yielding the usual universal number appearing at 
one loop). From R{u) one obtains, at least naively, the q — Q correlations to lowest order in e. 

Attempts to "improve" on this calculation however lead to difficulties. These occur already in considering the two 
point correlation function at non-zero momentum, which is determined from i?ii(0, 0, q). Naively solving (jI12p in the 
zero temperature limit gives the outer solution: 

Router{ui,U2,q) = (^(g) R" (ui ) R" (w2 ) , (114) 

which, we know from considerations in the present paper can at best hold only for Ui ~ 0(1). The momentum integral 
is found to be: 

/•OC 

(t){q)= dl'e'^'{h{qe'^')-h{0)), (115) 
Jo 

and h(q) = \ (^2+1)2 (fc+J)2+i- 

Taking two derivatives gives Rii{ui,U2, q) = (j>{q)R"' {ui)R"' {u2)- Unlike similar local (i?"'(u))^ terms, this does 
not have an unambiguous limit for small ui,U2. This is indeed a sign of difficulty for the matching. In particular, 
the naive Larkin term in the bilocal TBL, i?TSL(ui, U2, <?) ~ f2,2{q)uiU2, does not match the non-analytic outer 
solution at small argument Router{ui, W2, ?) ^ 0(9)(-R"'(O^))^|ui||m2|. Thus the matching problem is more complex in 
this case and requires further analysis. Without such a prescription, the zero temperature beta function becomes itself 
ambiguous at 0{R'^), since at this order there is a feeding term in the local equation of the form ^ R"(u)Rii(u, m, •), 
an expression which naively starts as |m| at small u, which seems problematic. Work is in progress to tackle these 
issues. , 
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d{—e)P{x2,xi,—e) being normalized to 2. The RG equa- 
tion 111911 is established, strictly, for P for any e, but 
clearly also applies to P for any e > 0. The symmetry 
of D{xi,X2) = P{xi,X2,e = 0"*") then follows if P is con- 
tinuous (which is implied by the RG equation for xi ^ x^)- 
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®^ the condition exactly for d = 4 becomes Tj ^ \^ and all 
other considerations are very similar 

®* Although as an extra precaution one may want to distin- 
guish T — > 0^ (studied here, liniT^oli'mi^oo) and T — 
(absolute minimum of the energy limi^aolirnT->o) , we will 
denote T = everywhere below for simplicity (also, at this 
stage we do not see any reason why the two should diff'er) . 

®^ in the simplest version of the droplet theory one assumes 
that p is uncorrelated from xi,X2- This restriction can be 
lifte d although it would cha nge the amplitudes 
that 6p2(l -p)2 = p{l-p) 



